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Abstract. 


I 


Modern  optimal  control  methods  are  used  to  develop  a 
multiple-lnput-multlple-output  controller.  Focus  is  made  on  a 
three-controller  configuration  exhibiting  high  controllability 
but  low  observability  in  the  first  controller,  a  median  amount 
of  each  in  the  second,  and  low  controllability  but  high 
observability  in  the  third.  These  characteristics  are  due  to 


V 

i 


the  technique  used  to  suppress  control  and  observation 
spillover  among  the  controllers.  A  control  model  for  large 
space  structures,  which  employs  full  state  feedback  using 


deterministic  observers,  is  developed  and  implemented  in  a 
computer  simulation.  The  technique  of  spillover  suppression 


and  the  conditions  assuring  stability  of  the  control  system 


are  developed  and  implemented  as  well. 

The  simulation  is  tailored  to  address  the  control  of  the 
Draper-2  large  flexible  space  structure  model.  The  model  has 


been  used  previously  for  optical  pointing  ( 1 ine-of -sight 
<LGS))  studies.  Here,  position  sensors  and  point  force 
actuators  are  used  to  effect  feedback  control  (regulation)  of 
the  damped  unforced  structural  vibrations.  The  simulation  can 
output  both  the  unsuppressed  and  suppressed  case  open- loop  and 


Hi 

& 


closed-loop  eigenvalues  and  the  LOS  time  response  for 
stability  and  performance  analysis. 

With  the  control  problem  formulated  for  modal  control. 


an  investigation  is  made  into  the  effects  on  time  response  of 


assigning  three  groups  of  four  modes  in  a  permutative  fashion 


1 


to  the  three  controllers.  A  fourth  residual  set  of  eight 
modes  is  carried  without  spillover  suppression  to  represent 
the  unntodelled  modes  of  a  real  structure.  The  groups  are 
assembled  based  upon  a  previous  investigation'' s  results  from 
applying  modal  cost  analysis  for  LOS  performance.  Simple  high 
frequency  truncation  is  also  used  for  comparison. 

Controllers  based  on  modal  cost  analysis  alone  are  found 
to  yield  marginal  stability  and  mediocre  LOS  performance  due 
to  little  insight  into  the  sensitivity  of  the  residual  modes 
to  spillover.  However,  specific  problem  modes  are  readily 
identified  by  examining  the  results  of  an  internal  balancing 
analysis  for  modal  sensitivity.  Simple  frequency  truncation 
is  found  to  give  the  best  time  response  here  when  the  modes 
most  contributory  to  LOS  are  assigned  to  the  controllers  with 
more  controllability.  However,  the  relatively  small  quantity 
of  modes  and  the  overwhelmingly  large  relative  contribution  of 
the  three  rigid  body  modes  included  may  obscure  some  conclu¬ 
sions.  Results  indicate  more  revealing  results  might  be  ob¬ 
tained  if  more  modes  are  added  to  the  model  and/or  if  some  of 
the  residual  modes  are  suppressed. 


MODAL  ASSIGNMENT  EFFECTS 
ON  DECENTRALIZED  CONTROL 
OF  A  LARGE  SPACE  STRUCTURE 

I .  Introduction 

Now  that  the  space  shuttle  Is  demonstrating  the  ability 
to  routinely  deliver  sizeable  payloads  to  space,  the  construe 
tion  of  large  apace  structures  for  commercial  manufacturing, 
scientific  research,  and  military  use  will  soon  be  reality. 
There  are  still  practical  economic  constraints,  however,  on 
the  amount  of  material  the  shuttle  can  deliver  for  construc¬ 
tion.  Efficient  structural  design  then  dictates  the  use  of 
thin,  lightweight  elements  in  trusslike  frameworks.  Although 
helping  eliminate  the  problem  of  delivering  materials,  this 
introduces  a  control  problem  due  to  the  undesirable  flexi¬ 
bility  inherent  in  such  structures.  Often  the  large  space 
structure  is  Intended  for  high  accuracy  pointing  of  antennas 
and  lasers.  Active  regulation  (control)  of  structural  vibra¬ 
tions  is  thus  required.  Unfortunately  such  structures  are 
lightly  damped  and  have  numerous  low-frequency  vibrational 
modes,  often  clumped  together,  within  the  bandwidth  of  prac¬ 
tical  control  systems.  Furthermore,  there  can  be  hundreds  of 
modes  to  be  concerned  with  even  after  techniques  (1,  2)  have 
been  applied  to  reduce  the  number  of  modes  from  infinity  for 
continuous  system,  to  thousands  for  a  finite  element  struc- 


tural  model,  to  those  the  designer  concludes  to  still  have 
significant  effect  on  the  stability  and  time  response  of  the 
control  system.  Modern  state-space  methods  of  control  theory 
lend  themselves  well,  however,  to  this  concept  of  a  discrete 
(finite)  structural  model.  The  more  modes  a  designer  can  keep 
in  the  model  of  the  system/space  structure,  the  less  inac¬ 
curate  it  will  be.  But  the  computational  burden  to  an  online 
computer,  functioning  as  the  controller,  grows  also. 

Therefore  the  designer  may  need  to  reduce  this  computa¬ 
tional  burden  by  employing  multiple  controllers.  Each  would 
bear  only  a  portion  of  the  computational  load,  and  ideally 
control  only  those  modes  assigned  to  it,  suppressing  the 
effects  of  other  modes.  In  other  words,  ideally  all  modes 
would  be  either  controlled  or  suppressed.  Practical  con¬ 
straints  can  still  prevent  there  being  enough  controllers 
and/or  few  enough  modes  to  permit  this,  however.  The  remain¬ 
ing  modes  are  called  residual  modes.  The  effect  of  these 
modes  on  the  system  is  called  observation  spillover  as  they 
contaminate  the  sensor  (observation)  data.  In  contrast,  the 
control  applied  to  the  controlled  modes  can  affect  the 
residual  modes.  This  is  called  control  spillover,  and  the 
residual  modes  may  drive  the  system  either  unstable  or  even 
more  stable.  Insofar  as  the  suppression  of  spillover  is 
concerned,  however.  Calico  and  Aldridge  (3)  successfully 
showed  that  a  transformation  matrix  could  be  generated  by 
singular  value  decomposition  of  the  control  and  observation 
matrices  which,  when  applied  to  the  feedback  gain  matrices. 


would  drive  spillover  terms  to  zero,  effecting  suppression. 

They  employed  a  mathematical  model  of  a  fictitious  nontrivial 
apace  structure  in  a  simulation  whose  computer  code  allowed 
for  the  inclusion  of  either  three  controllers  with  a  set  of 
residual  modes  (numbering  as  few  as  zero)  or  four  controllers 
with  no  residual  modes.  The  code  contained  an  observer  to 
provide  full  state  feedback.  Calico  and  Thyfault  (4)  used  a 
similar  version  with  altered  code  to  investigate  the  effect  of 
direct  output  feedback  in  place  of  the  observer.  In  neither 
case  was  there  time  to  carry  the  investigation  to  examining 
time  response,  which  could  have  been  quite  revealing  of  the 
overall  performance  effects  of  the  various  ways  in  which  modes 
are  assigned  to  controllers. 

—The  thrust  of  this  thesis  is  to  implement  time  response 
output  for  the  nontrivial  model  and  investigate  the  effect  on 
time  response  of  certain  modal  assignments  by  fixed  groups  to  any 
of  the  three  controllers  with  another  fixed  group  assigned  as 
residuals.  The  inclusion  of  residuals  provides  in  a  limited 
sense  a  truth  model  for  the  analysis  of  stability  and  per¬ 
formance.  There  may  be  significant  differences  in  time  response 
caused  by  the  fact  that  a  controller  is  known  to  provide  more 
relative  controllability  and  less  observability,  or  vice  versa, 
for  its  assigned  modes.  The  investigation  will  involve  initially 
the  fixing  of  certain  baseline  parameters  to  allow  a  parallel 
comparison  of  reasonable  results.  Then  line-of -sight  pointing 
and  defocus  performance  will  be  generated  for  comparison  among 
various  cases.  The  performance  is  actually  accomplished  by  using 


position  sensors  to  provide  modal  amplitudes  measurement  with  an 
observer  in  place  for  the  unmeasured  states,  and  using  point 
force  actuators  to  execute  control.  All  processes  are  assumed  to 
be  deterministic,  as  stochastic  estimation  and  control  is  not 
within  the  scope  of  this  investigation.  The  simulation  model 
is  a  version  of  the  so-called  ^*^raper-2'*  space  structure  model 
originated  at  the  Charles  Stark  Draper  Laboratories. 

The  next  section  -w^-tlr^descrlbej- the  selected  model  config¬ 
uration  and  discuss  its  finite  element  representation.  Then 
the  modal  control  and  matrix  transformation  methods  will  be 
explained.  After war d^~^he  implementation  of  the  simulation 
program  will  be  discussed.  -Finalli^  l^he  last  sections  will 
detail  the  investigation,  results,  conclusions  and  recom¬ 
mendations.  •  ,  ,  .  -  I  C.  ■- f-j'.  .  . 
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II .  Selected  Model  Configuration 


The  model  selected  is  basically  the  same  as  that  used  by 
Calico  at  al  in  their  investigations.  It  is  a  fictitious, 
nontrivial  space  structure  developed  by  Charles  Stark  Draper 
Laboratories  for  just  such  control  problem  studies.  Several 
versions  of  the  model,  incorporating  only  minor  changes,  have 
indeed  been  used  extensively  throughout  the  community.  The 
version  selected  here  is  that  used  by  Draper  Laboratories 
themselves  in  a  study  done  for  the  Defense  Advanced  Research 
Projects  Agency  (DoD)  (S) .  It  varies  slightly  from  that 
version  used  by  Calico  at  al  because  of  differing  masses  at 
some  locations.  It  was  selected,  however,  because  the  study 
using  it  Included  conversion  data  to  obtain  physical,  line-of- 
sight  output  directly  from  the  generated  modal  coordinates. 

The  model  provides  a  realistic  flexible  space  structure 
simulator  for  study,  and  includes  a  flexible  optical  support 
structure  in  a  trusslike  framework  and  an  equipment  section 
with  solar  panels.  The  optical  support  structure  consists  of 
three  trusses:  the  upper  mirror  truss  containing  a  convex 
primary  mirror  and  concave  tertiary  mirror,  the  lower  mirror 
truss  containing  the  flat  secondary  mirror  and  image  focussing 
plane,  and  the  metering  truss  to  maintain  mirror  separation. 
Figure  2-1  depicts  the  overall  structure,  and  Figure  2-2  the 
finite  element  representation.  Figure  2-3  gives  the  model's 
general  structural  dimensions.  Then  Figure  2-4  shows  the 
optical  path  of  a  ray  of  light  reflecting  from  the  mirrors. 
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fOUlfMiNT  SCCriON 


Appendix  A  gives  the  necessary  input  cards  for  a  NASTRAN 
run,  and  pertinent  data  for  the  Draper-2  model.  From  the  '‘C0NH2'' 
input  cards  it  can  been  seen  that  the  model  weighs  some  9304  kg, 
and  from  Fig  2-3,  is  24  m  high. 

Of  interest  in  this  investigation  is  the  calculation  of 
the  optical  line-of -sight  <L0S) .  For  Draper-2  the  optical 
line-of -sight  refers  to  how  the  ray  of  light  shown  in  Fig  2-4 
deviates  in  the  x-,  y-,  or  z-direction  from  a  nominal  "focal” 
point.  The  x-y  planar  deviations  can  be  (and  are  in  this 
investigation)  put  into  a  more  meaningful  radial  deviation, 
while  deviations  in  the  z-direction  can  be  understood  better 
as  a  measure  of  defocus  of  the  optical  ray.  The  motions  of 
the  structural  members,  of  course,  are  what  causes  these 
deviations,  and  these  motions  bear  kinematical  relationship  to 
each  other  since  the  members  are  connected.  Ref  (5>  provides 
the  basic  algorithm  for  calculating  these  LOS  deviations  or 
errors,  while  the  full  expansion  of  the  basic  algorithm  can  be 
found  in  Ref  (3) .  Fortunately  it  is  not  necessary  to  employ 
the  fullblown  set  of  equations  to  calculate  the  LOS  error  at 
some  time  t  in  the  time  history  of  the  motion  of  the  struc¬ 
ture.  When  modal  analysis  is  performed,  as  in  this  investiga¬ 
tion,  a  vector  of  modal  coordinates  at  some  time  t  can  be 
multiplied  by  some  output  matrix  which  effects  modal  super¬ 
position  and  yields  the  three  quantities  of  interest  in  an 
"LOS  vector."  This  matrix  is  referred  to  as  <t)L0S  «nd  is  of 
row  dimension  3  and  column  dimension  n,  where  n  is  the  number 
of  modes  in  the  model.  Table  2-1  provides  the  matrix  in  its 


transposed  form  for  convenient  reading.  The  column  labels 
LOSX  and  LOSY  refer  to  a  normalization  imposed  on  the  x-  and 
y-LOS  errors,  wherein  they  are  divided  by  a  factor  of  8.051. 


Table  2-1.  (5> 

Transpose  of  Matrix  <^L0S 

Mode  Free  LOSX  LOSY  Def  ocus 

1  0.0  0.2705E-04  0.3778E-03  0.9313E-09 


2 

0.0 

-.1580E-04 

3 

0.0 

0.4927E-03 

4 

0.0 

0.1166E-03 

5 

0.0 

-.3502E-04 

6 

0.0 

- .7835E-03 

7 

0.1455 

-.2641E-06 

8 

0.2832 

0.3489E-06 

9 

0.3173 

-.1505E-05 

10 

0 . 3329 

-.2293E-03 

11 

0.4432 

0.6041E-05 

12 

0.5779 

0.7359E-03 

13 

0.5814 

0.2330E-05 

14 

1.2238 

-.1027E-05 

15 

1.3002 

0.3759E-05 

18 

1 . 3475 

-.1589E-03 

17 

1.7209 

0.2738E-04 

18 

1.8187 

0.2590E-05 

19 

1.8187 

0.1151E-04 

20 

1.8892 

- .3597E-13 

21 

2.3635 

0.2237E-04 

22 

2 . 9895 

0. 1345E-03 

23 

3.1795 

- .4355E-05 

24 

3.3873 

0.5573E-06 

25 

5.1617 

- . 1814E-13 

26 

5.2603 

- .3395E-04 

27 

7.8769 

0. 1038E-14 

28 

8.1168 

-.4275E-03 

29 

8 . 3600 

0.2916E-03 

30 

8.5706 

0. 1868E-01 

31 

8.8135 

- .3822E-04 

32 

8.8135 

-.4494E-04 

33 

11.3462 

0.3762E-05 

34 

11.4978 

-.1225E-02 

35 

12.7258 

0. 1098E-01 

38 

13.5832 

0.2097E-04 

37 

13.7141 

-.7143E-03 

38 

14.1604 

0.3953E-02 

39 

15.6523 

-.1414E-01 

O . 8426E-06  - . 1863E-07 
-.3133E-04  -.9604E-09 
0.1155E-03  0.4366E-10 
-.8110E-03  -.2721E-08 
0.4676E-04  0.6985E-09 
0.3263E-03  -.4904E-07 
-.2141E-03  -.1099E-05 
-.3452E-06  -.5238E-05 
0.4407E-06  -.1467E-06 
0.1066E-05  0.1559E-04 
-.6768E-05  0.1906E-05 
0.8731E-03  -.4276E-06 
-.2362E-03  -.4831E-07 
-.1519E-03  -.7763E-05 
-.1413E-05  0.9364E-05 
-.2777E-02  -.4955E-04 
-.4976E-06  0.2634E-05 
-.2364E-05  -.5493E-06 
-.3518E-14  0.3030E-14 
0.1376E-01  0.9028E-05 
0.2724E-05  - . 1509E-04 
-.2602E-03  0.2500E-05 
-.4536E-03  -.4258E-06 
-.4398E-15  0.2601E-14 
-.7013E-06  0.4747E-05 
0.7264E-16  -.9614E-16 
-.1701E-02  -.1950E-03 
0.4775E-06  -.3638E-05 
-.1108E-04  0.1072E-02 
-.2919E-06  -.3755E-05 
0.7780E-07  0.9546E-06 
0.8995E-05  -.5240E-06 
0.1870E-02  0,3982E-03 
-.5662E-03  -.1593E-01 
0.7475E-03  -.6865E-04 
-.lOllE-01  0.1142E-02 
0 . 5935E-02  - . 1084E-02 
0.3213E-02  -.4966E-02 


Table 

2-1.  cont'd 

Transpose  of 

Matrix  (^LOS 

LOSX 

LOSY 

Def ocus 

40 

16.0724 

0.1907E-02 

0.1044E-01 

-.2294E-03 

41 

16.5248 

-.6780E-02 

-.2220E-02 

0.1413E-01 

42 

16.7453 

0.3a53E-02 

-.9540E-03 

-.1375E-01 

43 

17.1555 

0. 1132E-03 

-.8405E-02 

-.1739E-02 

44 

17.8283 

0.3244E-02 

-.1168E-02 

0.7174E-02 

45 

19.0713 

-.6439E-03 

-.1446E-01 

-.1545E-02 

46 

23.7716 

-.2266E-01 

0. 1669E-01 

0.5217E-01 

47 

24.4140 

0.9058E-02 

0.3711E-01 

- .2293E-01 

48 

25 . 9089 

0.2688E-03 

0.1205E-01 

-.3317E-02 

49 

26.3625 

-.1444E-01 

-. 1042E-01 

0.2382E-02 

50 

26.4292 

0.9416E-14 

-.6257E-14 

- . 1485E-14 

Note  that  the  first  six  rows  are  for  the  rigid  body  modes 
which  have  zero  frequency.  The  portion  of  the  table  used  in 
this  investigation  included  mode  numbers  4  through  23,  total¬ 
ling  twenty  modes. 

For  this  investigation*  a  set  of  twenty-one  collocated 
pairs  of  point  force  actuators  and  position  displacement 
sensors  were  employed  at  the  same  locations  in  the  structure 
as  in  Ref  (3) .  These  locations  and  the  sensor /actuator 
orientations  in  direction  cosines  are  given  in  Table  2-2. 

There  is  a  sensor  and  actuator  for  every  motion  which  affects 
the  LOS  deviations.  Notice  that  the  orientations  are  conven¬ 
iently  aligned  with  the  coordinate  axes  of  the  model. 

Strictly  speaking*  any  orientation  that  is  physically  possible 
could  be  used  as  long  as  proper  direction  cosines  are  used. 

With  the  structural  model  described*  the  system  mathema¬ 
tical  model  and  control  derivation  is  at  hand. 


Vi  J 


S' 


Table  2-2 

Draper-2 

Sensor/Actuator  Locations 

and  Direction  Cosines 

Pair 

Node 

X 

y 

z 

1 

9 

0 

1 

0 

2 

9 

0 

0 

1 

3 

10 

0 

0 

1 

4 

11 

1 

0 

0 

5 

11 

0 

1 

0 

6 

11 

0 

0 

1 

7 

12 

0 

0 

1 

S 

27 

1 

0 

0 

9 

27 

0 

1 

0 

10 

27 

0 

0 

1 

11 

28 

0 

0 

1 

12 

29 

0 

1 

0 

13 

29 

0 

0 

1 

14 

30 

0 

0 

1 

15 

32 

0 

0 

1 

16 

33 

0 

0 

1 

17 

34 

1 

0 

0 

18 

34 

0 

1 

0 

19 

34 

0 

0 

1 

20 

35 

0 

1 

0 

21 

35 

0 

0 

1 

<N . B .  These 

map  to 

only  translational 

components 

of  the  modal 

Matrix  described  in 

Section  III.) 
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The  motion  of  a  large  space  structure,  such  as  Draper-2, 
can  be  described  as  that  of  a  typical  multiple-degree-of - 
freedom  system  with  light  damping.  Here  the  system  will  be 
deterministic,  with  active  control  to  regulate  vibrations. 

The  equations  of  motion  can  be  expressed  in  matrix  form  as 

Mg  ♦  Eg  ♦  Kg  =  Du  <1) 

where 

H  -  nxn  mass  matrix 
E  -  nxn  damping  matrix 
K  -  nxn  stiffness  matrix 
D  -  nxm  actuator  mapping  matrix 
g  *  nxl  vector  of  physical  coordinates  q 
u  *  mxl  control  force  input  vector 
and  n  is  the  number  of  degrees  of  freedom  (corresponding  in 
turn  to  the  number  of  modes  of  vibration  in  the  system) ,  while 
m  is  the  number  of  actuators  employed  to  effect  control. 

The  matrices  M  and  K  may  be  obtained  by  finite  element 
analysis  of  the  structure  (for  example  by  NASTRAN) .  The 
control  designer  may  then  decide  how  many  modes  must  be  kept 
in  his  model,  and  how  many  and  what  type  actuators  there 
should  be  and  where  they  should  be  placed. 

Proceeding  now  toward  the  objective  of  determining  the 
control  u,  transformation  is  made  to  state  space  form. 


Solving  the  undamped,  f ree-vibration  equation  of  motion 


Mg  -  Kg  =  0 


yields  the  system's  natural  frequencies  k=l...n,  and  the 

nxn  modal  matrix  <{>.  The  columns  of  <I>  are  the  associated 
eigenvectors  determined  from  det  I  MC~w)<2j  -  K  I  =  0  .  The 

modal  matrix  is  such  that 


g  =  <t>!l 


where  2  nxl  vector  of  modal  coordinates.  The 

properties  of  <{>  are  such  that  when  Eq  (3)  is  directly 
substituted  into  Eq  (1)  which  is  then  premultiplied  by  the 


result  can  be  expressed  as 

\  1 .  r\  1 

'  L  \J'  L  \  ‘ 

The  second  order  matrix  equation  may  now  be  converted  to  a 


first  order  state  space  expression 


X  =  Ax  ♦  Bu 


where 


-  2nx2n  plant  matrix 

-  2nxm  control  mapping  or  input  matrix 


-  2nxl  state  vector 


u  -  mxl  control  input  vector 


C‘ 
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These  state  space  matrices  have  form 


0 

I 

\ 

\ 

-0)2 

-2  fu) 

[  \ 

subscript  k  implied) 


B  * 


O 


(6) 


X  s 


To  obtain  feedback  of  the  time  dependent  motion  of  the 
structure,  so  as  to  determine  the  control  vector  u,  an 
output  equation  is  introduced. 

Y  =  Cpg  +  Cvg  (”7) 


where 

Cp  -  npxn  displacement  matrix 
Cy  -  nyxn  velocity  matrix 

np  -  number  of  position  (displacement)  sensors 
ny  * 


number  of  position  rate  (velocity)  sensors 


In  state  space  form,  the  output  equation  becomes 

Y  =  Cx  (8> 

where 

C  *  C  Cp4)  I  Cv<I)  ]  (9) 

Control  Model 

The  complete  finite  element  model  of  the  structure  is 
described  now  by  the  2n-dimenaional  vector  x.  Recall  that  n 
represents  the  number  of  modes  in  the  model ,  already  reduced 
from  infinity  for  the  real  structure.  Many  of  the  high- 
frequency  modes  can  usually  be  discarded  as  of  no  consequence 
in  a  reduced  order  model.  As  previously  mentioned,  there  are 
techniques  available  to  aid  in  deciding  how  many  and  which 
modes  to  keep  to  maintain  some  desired  accuracy  in  reproducing 
the  dynamics  of  the  structure.  The  finalized  reduced  order 
model,  the  control  model,  may  still  contain  a  great  number  of 
modes,  enough  to  place  practical  constraint  on  an  online 
computer  attempting  to  do  realtime  computation  of  the  control 
u  as  well  as  model  the  plant  A  and  receive  feedback  of  the 
output  y.  The  computational  task  may  be  apportioned  among  N 
controllers  reducing  the  unwieldy  sizes  of  matrices  to  multi¬ 
ply.  Doing  this  also  means  more  modes  can  be  controlled 
overall  than  by  using  a  single  controller.  As  is  the  case  in 
this  investigation,  each  of  the  N  controllers  would  control  a 
subset  ni  of  the  n  modes,  where  n  now  represents  the  actual 


C  I  Cvi^^  3 


(13) 


I 

I 


I 


c 


# 


t- 


& 


Ci  = 


Note  that  the  modal  components  o£  xi  (and  xr>  are  simply 
arbitrary  subsets  or  groupings  of  the  modal  states  without 
regard  to  any  particular  ordering.  In  other  words,  the  modes 
can  be  assigned  in  arbitrary  groups  to  the  various  controllers 
or  to  the  residual  set  in  the  control  model . 

Interestingly,  when  only  position  sensors  are  used  and 
these  are  collocated  with  the  actuators  in  pairs,  as  in  this 
investigation,  the  right  partition  of  the  Ci  matrices  is  zero. 
Dropping  the  subscript  i  momentarily  it  is  seen  that 

ccp<ia  =  c<J>Td]T  a  cdT(I>]  (14> 

This  is  convenient,  because  now  obviously  Cp  =  0*^. 

These  matrix  forma  and  equations  can  be  implemented  in  a 
computer  simulation  for  any  structure.  Only  the  dimensioning 
of  the  matrices  would  vary  based  upon  the  desired  number  of 
modes  to  keep  in  the  model  and  the  applicable  number  of 
sensor 3/ actuators . 

Modal  Control 

As  previously  indicated  the  type  of  control  in  this 
investigation  is  full  state  feedback  control,  specifically 
applied  to  a  multiple-controller  configuration.  The  desired 
control  is  given  by: 

N 

y  =  E  Gi2S1 

i  =  l 
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where  the  Gi  are  the  control  gain  matrices  to  be  determined 


Feedback  in  the  form  of  Eq  (15)  assumes  complete  know¬ 
ledge  of  the  state  variables.  This  is  impossible  in  a  real 
structure.  Here,  the  system  model  representing  the  real  struc¬ 
ture  is,  however,  completely  observable.  This  means  all  the 
states  can  be  reconstructed  by  use  of  an  observer  to  “observe" 
the  unmeasured  states.  <6)  The  state  equations  become 

«i  =  Ai«i  ♦  Biu  ♦  Ki<y  -  ^i>  (16) 

9i  =  Ci«i  (17) 

The  jti's  are  the  reconstructed  states,  the  y's  the  estimated 
outputs,  and  the  Ki's  are  observer  gain  matrices. 

The  objective  in  selecting  the  Ki's  is  to  have  the 
reconstruction  error 

§1  =  5li  -  xi  (18) 

go  asymptotically  to  zero  in  steady  state.  An  idea  of  how 
this  can  be  done  can  be  obtained  by  examining  the  case  of  a 
single  controller.  Here,  y  would  just  be  yi.  Using  this  fact 
and  subtracting  Eq  (16)  from  Eq  (10),  then  substituting  Eq  (17) 
and  Eq  (8)  (subscripted),  the  observer-error  state  equation  is 

§1  =  (Ai  -  KiCi)ei  (19) 

The  Ki  are  selected  so  that  the  eigenvalues  of  (At  -  KiCi)  are 
stable.  Thus  the  steady  state  error  will  go  to  zero.  Notice 
in  Eq  (19)  that  there  is  no  input  term.  The  error  is  excited 
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& 


41 


only  by  the  initial  conditiona  for  the  plant,  and  is  thus  not 
affected  by  control  inputs.  This  is,  however,  only  the  case 
when  there  is  a  single  controller.  With  multiple  controllers, 
Eq  (16)  still  applies  and  there  will  be  spillover  due  to 
output  y  from  all  sensors. 

Proceeding,  with  all  the  states  now  available  in  the 
vector  Hi,  the  control  is 


ui  =  53 

i  =  l 


(20) 


Determination  of  the  control  gain  matrices  as  well  as  the 
observer  gain  matrices  will  now  be  discussed. 

For  the  development  of  these  gain  matrices  linear  optimal 
regulator  theory  is  applied.  In  both  cases  this  involves 
construction  of  a  quadratic  performance  index  J.  For  a  control 
gain  matrix  the  index  is 


J  2  1/2 


r- 


iTQiXi  ♦  yTRj^u)  dt 


(21) 


where 


Q  -  a  positive  semidefinite  square  weighting  matrix 
R  *  a  positive  definite  square  weighting  matrix 

J  represents  a  compromise  between  obtaining  minimum  deviation 
of  the  present  state  from  the  desired  final  value  and  reaching 
the  desired  final  value  with  minimum  control  input  energy. 

Thus  the  objective  is  to  minimize  J  subject  to  the  state 
Eq  (10).  The  optimal  solution  is  shown  in  (6)  to  be  given  by 
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-Ri-^BiTSi 


(22) 


9^ 


where  Si  ia  the  solution  to  the  steady  state  matrix  Riccati 
equation : 


SiAt  ♦  AiTSi  -  SiBiRi-lBtTSi  ♦  Qi  »  0 


(23) 


A  similar  approach  is  taken  to  finding  the  observer  gain 
matrices  Ki.  To  build  a  similar  performance  index  expression, 
an  equation  similar  to  Eq  (19)  can  be  written: 


wi  =  (Ai  -  KiCi)Twj^ 


(24) 


The  eigenvalues  for  (Ai  -  KiCi)'!*  are  the  same  as  for  those  in 
Eq  (19)  since  only  the  transpose  of  the  matrix  is  involved. 

Eq  (24)  can  be  rewritten 


Si  *  Ai^wi  -  CiTri 


(25) 


which  provides  the  form  of  an  input 


ri  =  KiTwj 


(26) 


Now,  a  similar  performance  index  can  be  written 

r 

J  *  1/2  I  (wiTQoiWi  *  riTPoiri)  dt 
Jo 

where  the  weighting  matrices  are  not  necessarily  the  same  as 


(27) 


those  for  the  controller  gain  matrices. 

When  Eq  (27)  is  minimized  subject  to  Eq  (25)  the  solution 
is  optimally 


Ki  =  *Roi-lCiPi 


(28) 
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where  Pi  is  the  solution  to  the  Riccati  equation: 


P 1 A 1  AiPi  ■"  P iCi'^RQl ^ ^Cj^Pi  ♦  Qoi  ”  ^  (29) 

Eq  (22)  and  (28)  provide  the  forms  for  the  control  and 
observer  gain  matrices.  The  designer  may  then  determine  the 
weightings  “Q*‘  and  "R“  as  desired.  It  may  be  desired  to 
improve  response  settling  time,  stabilize  a  control- 
destabilized  system,  or  build  uniformity  into  the  closed-loop 
system's  modal  damping  ratios. 

While  each  separate  controller  may  be  designed  such  that 
it  is  stable,  the  individual  controllers  are  still  coupled 
when  applied  simultaneously  as  in  Eq  (20).  In  other  words, 
there  is  spillover  among  the  controllers.  The  next  subsection 
will  describe  how  it  arises  and  indicate  what  is  to  be  done  to 
suppress  it  and  provide  stable  control  in  its  presence. 

Section  IV  will  detail  the  means  of  suppression. 

Multiple  Controllers 

Repeated  here  for  clarity  are  the  pertinent  equations  for 
multiple  controllers: 


xi  =  AiXi  Biy  (10) 

N 

y  =  ^  CiXi  Crxr  ( 12) 

i  =  l 

=  Ai«t  ♦  BiU  >  Ki(y  -  yi)  (16) 

9i  =  Ci«i  (17) 

ei  =  «i  -  Xi  (18) 


Subtracting  Eq  (16)  from  Eq  (10),  rearranging,  and  finally 
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substituting  Eqs  (12),  (17),  and  (18),  a  proper  expression  for 


the  error  states  associated  with  multiple  controllers  follows: 

r 

§i  “  -  KiCi)ei  ♦  ^  KiCicXk  i  =  1,...,N  (30) 

k=l  k  =  1 , . . . ,N,r 

k9«i 

Using  Eq  (15),  the  state  equations  for  controlled  states,  with 
error  states  included,  as  well  as  for  the  residual  states  are 

N 

Xi  =  (Ai  BiGi)xi  ♦  BiGigi  ♦  ^  Bt^Gijl)<  i  =  1 . N  (31) 

k=l 

kj<i 

and 

N 

Xr  “  A^X}^  ♦  B]pGi)li  (31a) 

i  =  l 


Notice  there  is  no  residual  error  state  equation  because  no 
attempt  is  made  to  observe  the  residual  states. 

An  augmented  state  vector  z  can  now  be  assembled  from 
the  controlled  state  vectors  xi  and  their  associated  error 
states  ei,  wherein  Eqs  (30)  and  (31)  are  combined  so  that 
the  state  and  error  vectors  for  individual  controllers  are 
kept  ad3acent: 

z  =  (  xiT,  eiT,  ...  xnT,  eN^»  (32) 


With  Eqs  (30)  and  (31)  so  combined,  a  single  closed-loop 
matrix  state  equation  can  be  assembled  in  terms  of  z 


BiGi 

. .  BiGx 

BiGi 

. .  BiGn 

BiGn 

0 

0 

• 

• 

Ai-KiCi 

a 

• 

..  KiCi 

a 

a 

0 

a 

a 

• •  KiCn 

a 

a 

0 

a 

a 

KiCr 

a 

z  = 

• 

• 

BiGi 

■ 

a 

BiGi 

a 

. . Ai ♦BiGi 

a 

a 

BiGi 

a 

a 

. .  BiGN 

0 

a 

a 

0 

z 

KiCi 

0 

.  .  0 

Ai-Kid 

..  KiCN 

0 

KiCr 

• 

• 

■ 

; 

: 

: 

:33> 

BnGi 

BnGi 

•  •  BNGi 

BfiGi 

.  .  An^BnGn 

BnGn 

0 

KnCi 

0 

..  KNGi 

0 

. .  0  An 

-KnCn 

KNCr 

BrGi 

BrGi 

. .  BpGi 

BrGi 

. .  BrGN 

BrGN 

Ar 

The  entire  control  model  is  now  contained  In  thla  single  first 
order  matrix  differential  equation.  The  large  matrix  above 
has  been  referred  to  as  the  system  ma^or  matrix.  In  its 
present  general  form  it  la  not  easy  to  see  if  such  a  system  is 
stable  or  not,  but  this  will  soon  be  examined. 

Now,  for  three  controllers  N  =  3  ,  and  the  augmented 

state  vector  appears 

z  =  (  xiT,  eiT,  >  62^.  egT,  (34) 


The  components  of  z  are  then  assembled  based  upon 


The  observer  applied  is  of  form 


iii  =  AiJti  *  Biu  +  Ki(Y  -  ?i>  i  =  1,2,3 


where 


^i  =  CiJli  i  =  1,2,3 
u  =  Giili  G2JI2  *  ^3^3 


(39) 


(40) 


Applying  properly  chosen  Ki  such  that  the  individual  ei  go  to 
zero  with  time,  the  expressions  for  the  ei's  are 


®1  =  <Ai-KiCi)ei  K1C2X2  ♦  K1C3X3  +  KiCfXr 
®2  *  <A2“K2C2)«2  *  K2C1X1  ♦  K2C3X3  K2CrXr 
§3  *  (A3-K3C3)e3  *  K3C1X1  ♦  K3C2X2  ♦  KgCrXr 


(41) 


(42) 


(43) 


Recalling  that  ei  =  ill  ~  Xi  ,  substituting  this  properly 
into  Eq  (40) ,  and  then  combining  Eq  (40)  properly  with  Eqs 
(35),  (36),  and  (37),  the  remaining  components  of  z  are  now 

obtained : 


XI  = 

(Ai>BiGi)xi  + 

BiGiei  + 

BiG2i^2  * 

B1G3X3 

II 

•XI 

( A2-*’B2G2 )  X2  * 

B2G2e2 

B2Gi«i  ^ 

B2G3i^3 

X3  = 

(A3+B3G3)X3  ♦ 

BgGgeg  ♦ 

B3G1JI1  ♦ 

B3G2^2 

xr  = 

ArXf  BrGiJli 

♦  BrG2»2 

♦  BrG3»3 

(44) 


(46) 


(47) 


Equations  (41)  through  (47)  can  now  be  put  into  matrix  form 
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r 

Ai +B1G1 

BlGi 

B1G2 

B1G2 

B1G3 

B1G3 

0" 

0 

Ai+KiCi 

K1C2 

0 

K1C3 

0 

KlCr 

B2G1 

B2G1 

A2*B2G2 

B2G2 

B2G3 

B2G3 

0 

K2C1 

0 

0 

A2-K2C2 

K2C3 

0 

K2Cr 

z 

B3G1 

B3G1 

B3G2 

.  B3G2 

A3+B3G3 

B3G3 

0 

K3C1 

0 

K3C2 

0 

0 

A3-K3C3 

K3Cr 

_BrGi 

BrGl 

BrG2 

BrG2 

BrG3 

BrG3 

ArJ 

(48) 

It  is  now  readily  apparent  how  spillover  can  cause  overall 
instability  in  the  system.  For  example,  notice  how  the 
control  gains  G2  and  G3  make  up  part  of  the  component  state 
equation  for  controller  no.  1  (row  1>.  The  same  type  of 
spillover  from  the  two  other  control  gain  terms  appears  in  the 
component  state  equations  for  controller  no.'s  2  and  3  (rows  3 
and  5) .  Furthermore  all  the  control  gains  appear  in  the 
component  state  equation  for  the  residual  modes,  to  which  no 
control  was  to  be  applied.  If  all  such  terms  could  be  elimin¬ 
ated,  spillover  would  not  occur.  It  is  not  really  necessary, 
however,  to  eliminate  all  such  terms. 

What  is  needed  is  a  block  triangular  form  for  the  ma^or 
matrix.  Since  the  diagonal  terms  of  the  ma^or  matrix  already 
have  "stable"  system  eigenvalues  by  design,  it  is  only  neces¬ 
sary  that  the  matrix  be  block  triangular.  This  is  true 
because  the  eigenvalues  of  a  block  triangular  matrix  are  the 
same  as  those  of  its  diagonal  blocks.  Thus  two  types  of 
triangular ization  can  be  made  to  eliminate  sufficient  block 
terms  and  effect  suppression  of  spillover.  An  upper  block 


triangularlzation  for  suppression  would  require  the  following: 


BtGk  =  KiCk  =0  k  =  (49) 

i  =  k-t-l ,  .  .  .  , N 

A  lower  block  triangularlzation  would  require 

BiGvt  =  KtCk  =  0  i  =  1,...,N-1  (50) 

k  =  L*l,.  .  .  ,N 

Notice  the  residual  terms  are  ignored.  This  means  spillover 
between  the  residual  modes  and  the  controlled  modes  is  not 
suppressed  and  the  system  stability  and/or  control  performance 
may  be  affected.  Table  3-1  shows  which  terms  to  eliminate  for 
upper  or  lower  block  triangularlzation. 

Table  3-1 

Suppressing  Spillover  Among  Three  Controllers 
Upper  Trianoularization  Lower  Triangularization 


B2G1  =  0 

B1G2 

=  0 

BgGi  =  0 

B1G3 

=  0 

B3G2  =  0 

B2G3 

=  0 

K2C1  =  0 

K1C2 

=  0 

K3C1  =  0 

K1C3 

=  0 

K3C2  =  0 

K2C3 

=  0 

n>; 


e 


Before  discussing  the  technique  of  triangularization 
involving  transformation  matrices,  mention  should  be  made  of 
an  alternative  way  to  effect  triangularization.  Aldridge  (3) 
pointed  out  that  it  had  been  shown  previously  that  by  proper 


modes  to  controllers,  the  three-controller  problem  could  be 
reduced  to  a  two-controller  problem.  This  comes  about  through 
the  changes  effected  in  the  control  mapping  matrix  B  (for 
actuation)  and  the  output  mapping  matrix  C  (for  sensing) . 

With  proper  modal  asaignments  one  controller  could  be  made 
orthogonal  to  the  other  two,  thereby  eliminating  spillover 
for  that  controller.  No  attempt  will  be  made  in  thijs  invest¬ 
igation  to  configure  the  structural  model  for  this.  However, 
such  a  discussion  serves  to  introduce  the  next  subsection. 

The  above  discussion  on  sensors/actuators  introduces  the 
idea  that  selection  of  them  impacts  the  control  problem. 
Earlier  it  was  pointed  out  that  an  observer  was  clearly  neces¬ 
sary  due  to  an  insufficient  quantity  of  sensors  to  feed  back 
all  the  states.  The  next  subsection  will  detail  how,  due  to 
the  triangularization,  there  is  Indeed  a  minimum  number  of 
sensors  and  actuators  required  for  complete  observability  and 
controllability,  which  are  themselves  required  for  the  forma¬ 
tion  of  the  observer  and  control  gain  matrices. 


Sensor/ Actuator  Requirements 

Control  or  observation  spillover  among  the  various 
controllers  may  be  eliminated  by  constraining  the  various 
controller  and  observer  gain  matrices.  Using  the  control  term 
BiGi  as  an  example,  notice  that  if  the  columns  of  Gi  are 
constrained  to  lie  in  the  null  space  of  Bi,  then  the  product 
is  zero.  The  product  could  vanish  for  all  values  of  i;<l  if  Gi 
lay  in  the  null  space  of  the  matrix  B2N  defined  as  follows: 
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The  null  space  of  this  matrix  B2N  dimension  P2N  defined  as 


P2N  =  <na  -  r2N> 


(52) 


where 

na  -  number  of  actuators 

r2N  *  rank  of  B2N  <  min(  n2  ♦  ns  ...  ♦  njij,  na> 

Since  the  columns  of  Gi  lie  in  the  null  space  of  B2N  of 
dimension  P2N»  the  number  of  linearly  independent  columns  of 
Gi  is  also  P2N« 

As  Aldridge  pointed  out,  the  number  of  actuators  must 
exceed  the  rank  of  B2N  to  have  B2NG1  =  O  .  Without  this 
condition  being  met,  no  transformation  matrix  to  do  this 
exists.  If  B2N  full  rank,  the  number  of  actuators  needed 
will  be 

N 

na  >  ni  (53) 

i  =  2 

whereas  if  it  is  rank  deficient,  then 

na  >  r2N  (54) 

Similarly,  the  required  number  of  sensors  is,  for  some  C2N 


with  linearly  independent  columns 


N-1 

ns  >  S  ni  (55) 

i  =  l 

and  for  linearly  dependent  columns, 

ns  >  r2N  (here  r2N  the  rank  of  C2N>  (56) 

The  above  are  for  an  upper  block  triangularization .  Noting 
the  conditions  of  Eq  (50) ,  the  required  numbers  of  actuators 
and  sensors  for  a  lower  block  triangularization,  assuming  B2N 
and  C2N  ax'®  full  rank,  are  respectively 

N-1 

na  >  ni  (57) 

i  =  l 

N 

ns  >  ^  ni  (58) 

i  =  2 

Satisfying  Ineq  (51)  through  (58)  provides  complete 
controllability  and  observability  of  the  controlled  modes  and 
ensures  that  spillover  among  them  is  suppressed.  To  give  an 
example,  if  the  finalized  control  model  has  three  controllers 
each  containing  ten  modes,  and  there  are  no  residual  modes,  at 
least  twenty-one  sensors  and  twenty-one  actuators  are  needed 
to  completely  control  and  observe  such  a  system.  This  inves¬ 
tigation  satisfies  the  requirements  on  sensors  and  actuators. 

Based  upon  the  above  conditions  being  met,  the  technique 
of  generating  the  transformation  matrices  follows. 


w. 


IV .  Transformation  Technique 


This  section  describes  how  the  ma^or  matrix  can  be  block 
triangularized,  thereby  eliminating  the  observer  and  control 
spillover  terms  among  controllers.  The  discussion  will 
describe  the  determination  of  transformation  matrices  to 
eliminate  the  control  and  observation  spillover  terms.  These 
matrices  are  referred  to  respectively  as  T  and  L. 

For  a  configuration  using  a  single  controller,  the 
spillover  terms  to  be  eliminated  are  BsG  and  KCs.  The  sub¬ 
script  s  refers  to  those  modes  to  be  suppressed.  The  trivial 
solution  that  G  =  K  ^  O  is  unacceptable  because  it  means  the 
control  and  observer  feedback  loops  have  effectively  been 
severed  by  turning  off  all  gain.  The  system  would  be  without 
feedback  control  inputs,  and  therefore  would  simply  return  to 
the  case  of  a  structure  regulating  itself  through  natural 
damping.  The  objective  here  is  to  find  a  solution  (transfor¬ 
mation)  such  that  BsG  =  0  and  B^G  ^  O  ,  as  well  as  KCs 
=  0  and  KCc  ^  0 

Before  proceeding,  it  is  worth  mentioning  here  that  some 
attempt  could  be  made  to  include  suppression  of  some  or  all  of 
the  residual  modes.  Varhola  (7)  actually  did  this  for  a 
smaller  order  model  with  a  single  controller.  Such 
suppression  will  not  be  addressed  in  this  investigation, 
however,  as  it  will  be  assumed  that  the  multiple  controllers 
already  have  enough  to  handle  in  suppressing  spillover  among 
themselves.  The  residual  modes  are  presumably  beyond  the 
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capacity  for  control  or  suppression,  and  in  fact  may  not  even 
be  known  in  a  real  structure. 

Returning  now  to  discussing  multiple  controllers,  and 
bearing  in  mind  the  conditions  cited  above,  attention  will  now 
be  focused  on  control  spillover  suppression.  The  Bs  matrices 
can  take  the  form  of  Eq  (51)  repeated  here: 


B2N  = 


(51> 


Recall  this  form  was  baaed  upon  the  condition  of  Eq  (49>  for 
upper  block  triangularization : 


Bi^k  “  KiC'n  =  O  k  =  1,...,N~1 

i  =  1.  .  .  .  ,N 


(49) 


saying  that,  for  example,  given  N  controllers,  the  columns  of 
Gi  must  be  made  orthogonal  to  the  (N-1)  rows  of  the  Bi 
matrices  which  were  combined  into  B2N  above.  This  says  that 
the  columns  of  Gi  are  constrained  to  lie  in  the  null  space  of 
B2N*  This  is  ^ust  an  example  of  the  what  the  objective  is. 
The  discussion  will  continue  using  the  generic  term  Bs  to 
represent  any  control  mapping  matrix  (or  set)  that  does  not 
match  in  index  the  control  gain  matrix  of  interest. 


The  control  transformation  matrix  T  is  such  that 


It  will  soon  be  apparent  why,  but  Eq  (59)  is  introduced  here 
to  point  out  again  the  conditions  for  the  existence  of  T. 

From  Eq  (59)  it  is  seen  that  the  row  dimension  of  T  must  match 
the  column  dimension  of  Bg,  which  is  na«  the  number  of 
actuators.  The  column  dimension  of  T  is  found  by  referring 
to  the  conditions  cited  in  Section  III;  namely,  its  column 
dimension  must  be 

P2N  =  <na  -  r2N)  <52) 

where  it  is  seen  that  r2N  must  be  at  least  as  big  the  number 
of  modes  to  be  suppressed:  here  na-  Thus  the  dimensions  of  T 
are  na  by  (na~na)»  given  the  conditions  cited  on  linear  inde¬ 
pendence  of  the  rows  of  83. 

The  application  of  T  will  now  be  illustrated.  Consider  a 
system  containing  controlled  and  suppressed  modes  as  follows: 

Xc  =  AcXc  *  Bcu  (60) 

Xs  =  AsXs  *  BsU  (61) 

where 

u  —  ^  ) 

It  is  necessary  to  eliminate  the  BsU  term  of  Eq  61)  so  that 
the  system  described  by  Eq  61)  no  longer  admits  to  any  of  the 
control  applied  to  the  system  of  Eq  (60).  To  eliminate  this 
term  it  is  necessary  to  find  a  matrix  T  such  that  BsT  =  0 


while  BcT  ^  0  .  Let  this  transformation  matrix  T  be  used  to 

define  a  new  control  v,  such  that 

u  =  Tv  (63) 

Substituting  Eq  (63)  into  Eqa  (60)  and  (61),  and  renaming  the 
term  BcT  as  B*  while  recognizing  that  the  term  BsT  is  now  ^ust 
zero,  results  in  a  transformed  system  in  which  spillover  has 
been  suppressed: 


Xc  =  AcXc  ♦  B*y 

(64) 

Xa  “  AsXs 

(65) 

In  the  same  manner  as  the  construction  was  made  of  the 
feedback  control  vector  u,  the  new  control  vector  is 

y  =  G*«c 

The  determination  of  T  will  now  be  discussed. 

T  is  determined  through  singular  value  decomposition  of 
the  matrix  Bs  (3),  which  has  dimensions  x  n^g .  The  matrix 
can  be  decomposed  as  a  product  of  three  matrices: 

Bs  -  W  2  VT  (67) 

where 

W  ^  an  nmxnm  orthogonal  matrix  of  left  singular  vectors 

V  A  an  naxna  orthogonal  matrix  of  right  singular  vectors 

£  -  an  n^xna  matrix  containing  the  s  singular  values  of 


Bs  in  the  following  form: 


S  :  0 


0  :  O 


where 


.  0 


The  number  of  non-zero  singular  values  s  of  Bs  equals  its 
rank,  and  they  are  all  non-negative.  If  Ba  is  of  full  rank, 
then  a  -  mln(na,  nm) 

The  matrix  U  is  partitioned  into 


W  =  C  Wg  I  Wr  ] 


where 


-  an  n^xs  matrix  of  left  singular  vectors  associated 
with  the  non-zero  singular  values 

-  an  nmxr  matrix  of  left  singular  vectors  associated 
with  the  zero  singular  values 


The  V  matrix  is  partitioned  in  a  like  manner  as 


V  =  C  Va  I  Vp  ] 


where 


m 


Vs  *  an  n^xs  matrix  of  right  singular  vectors 

associated  with  the  non-zero  singular  values 
Vp  -  a  n^xp  matrix  of  right  singular  vectors  associated 
with  the  zero  singular  values 

and 

na  =  3  *  p  <73> 

Since  the  V  matrix  is  an  orthogonal  matrix,  see  that 

VgTvp  =  0  (74) 

allowing  the  following  construction: 

BsT  *  0  =  WsSysTr  (75) 

Examining  Eqs  (73)  and  (74)  for  similarities,  see  that 

T  =  Vp  (76) 

would  satisfy  the  equality  baaed  upon  Eq  (74) .  Necessarily 
Vp  ^  0  ,  else  the  transformation  matrix  is  just  a  matrix  of 

zeroes--an  undesirable  result.  Now  that  the  transformation 
matrix  T  has  been  determined,  it  remains  to  see  how  the  new 
control  gain  matrix  G*  is  found. 

The  G*  matrix  is  found  in  the  same  manner  as  described  in 
Section  III  in  the  subsection  on  Modal  Control.  It  is  the 
optimal  solution 

Gi*  =  - [Rt*] -1 CBt*]TSi  (77) 

in  which  St  is  the  solution  to  the  steady  state  algebraic 


matrix  Riccati  equation: 

SiAi  +  AiTSi  -  SiBi* CRi*] -1  -01=0  <78> 

wherein  the  subscript  i  refers  again  to  a  particular  controller. 
With  the  appropriate  substitutions,  it  is  now  apparent  that 
the  form  of  Eq  (64>  can  be  reiterated  to  show  the  form  of  the 
state  equation  for  any  controller  i,  after  suppression: 

Si  ®  AiXi  ♦  BiTiGi«ili  (79) 

Similarly,  substituting  for  B3,  K'^  for  G,  and  L  for  T  in 

the  above  development,  results  in  L  equalling  Vp.  Again,  the 
number  of  sensors  must  exceed  the  number  of  modes  to  be  sup¬ 
pressed  . 

This  section  completes  the  description  of  equations  used 
in  developing  the  control  model  for  implementing  in  a  computer 
simulation.  A  brief  discussion  on  how  the  time  response  is 
computed  is  included  in  the  following  section,  which  describes 
the  actual  computer  program. 


Computer  Simulat.ion  Implementation 


This  section  will  describe  the  operation  of  a  FORTRAN 
coded  computer  program  which  implements  the  equations  of  Sec¬ 
tion  III.  The  program  is  an  updated  version  of  that  used  by 
Aldridge.  Its  capabilities  have  been  altered  to  deliver  time 
response  of  the  LOS  deviations  mentioned  in  Section  II.  As  is 
often  the  case,  a  lengthy  amount  of  time  in  this  investigation 
was  devoted  to  understanding  what  the  program  does,  how  it 
does  it,  and  what  changes  could  be  made  to  enable  it  to  output 
LOS  time  response  data.  The  following  discussion  will  briefly 
describe  how  time  response  is  generated. 

The  theory  behind  obtaining  any  form  of  time  response  is 
fairly  straightforward.  Hence  it  is  included  here  instead  of 
Section  III.  Of  interest  is  Eq  <33>,  abbreviated  here  as 

z  =  CMAJMlz  C80) 

where  the  term  MAJM  replaces  the  large  block  matrix,  called  the 
major  matrix,  shown  In  Eq  (33).  This  is  a  commonplace 
homogeneous  first  order  matrix  linear  differential  equation. 

The  solution  of  a  scalar  first  order  differential  equation 
z  =  az  is  simply  exp(at),  wherein  t  represents  some  time 
differential.  Analogous  is  the  matrix  solution  to  Eq  <30) : 
exp ( CMAJM] t )  .  To  use  some  common  terms,  the  state  z  can  be 

found  for  any  time  tf  by  propagating  the  state  at  a  previous 
time  ti  forward  (or  background)  through  multiplication  by  the 
fundamental  matrix  (8),  known  also  as  the  state  transition  matri 


here  called  CSTM3 .  where 


CSTM]  =  exp( CMAJM] <tf  -  ti))  (81) 

Thus 

zf  =  CSTMJzi  <82) 

saying  that  a  fixed  time  differential  t  can  be  used  to  propa¬ 
gate  the  state  incrementally  ad  infinitum. 

Once  the  state  vector  can  be  output  at  discrete  time 
points,  it  is  only  a  matter  of  converting  the  state  at  a 
desired  stopping  point  from  modal  coordinates  to  physical 
coordinates.  The  desired  form  of  the  physical  coordinates  in 
this  investigation  is  the  LOS  vector  described  in  Section  II. 
This  version  of  the  program  outputs  a  tabulation  of  LOSX, 

LOSY,  Defocus,  and  Radius  vs.  time  with  the  rest  of  the 
printed  output,  as  well  as  3ust  Radius  and  Defocus  vs.  time  on 
another  output  file  for  plotting.  A  description  of  the  pro¬ 
gram's  capabilities  now  follows. 

The  capabilities  of  the  program  were  tailored  to  the 
investigation.  The  :ser  must  obtain  the  appropriate  products 
of  the  actuator  and  sensor  mapping  matrices  with  the  modal 
matrix  beforehand  (see  Eq  (14)).  Ref  (7)  provides  the  code  to 
read  in  any  modal,  sensor  mapping,  and  actuator  mapping 
matrices.  The  program  accepts  only  modal  coordinates  for 
initial  conditions,  although  once  again,  if  reasonable  physi¬ 
cal  initial  conditions  could  be  specified.  Ref  (7)  contains 
code  for  their  input  and  conversion  to  modal  coordinates. 


The  Q  weighting  matrix  is  assumed  to  be  diagonal  and  therefore 
only  its  diagonal  is  to  be  input.  A  complete  matrix  could  be 
formed  simply  enough,  but  was  not  of  interest  in  this  investi¬ 
gation.  Also  the  control  weighting  matrix,  R,  is  assumed  to  be 
an  identity  matrix.  Extensive  revision  of.  the  code  would  be 
needed  to  include  a  non- identity  R  matrix  if  it  were  of  interest. 
These  limitations  were  of  little  consequence  to  this  investigation. 

The  program  can  operate  with  either  three  or  four  control¬ 
lers  specified.  If  three  controllers  are  specified,  then  a  group 
of  modes  can  be  designated  as  residuals  in  place  of  the  group 
that  a  fourth  controller  would  have  controlled.  This  group  may 
number  as  few  as  zero  (no  residuals).  Thus  three  controllers 
with  or  without  residuals  or  four  controllers  with  no  residuals 
can  be  run.  In  either  case,  the  user  may  also  specify  other 
options:  (1)  time  response  (with  LOS  as  output),  (2)  roots 

(eigenvalues)  of  the  closed-loop  (overall)  system  and  the  open- 
loop  (individual)  syatema--control led ,  error,  and  residual 
states,  or  (3)  both  time  response  and  eigenvalue  analysis. 

Moreover,  the  program  makes  both  an  "unsuppressed"  and  a 
"suppressed"  pass,  and  any  of  these  three  options  can  be  selected 
for  each  pass  as  desired.  A  sensible  choice  in  the  early  part  of 
an  investigation  would  be  to  select  "2"  for  both  the  unsuppressed 
and  suppressed  passes  until  the  user  is  satisfied  with  the  place¬ 
ment  of  the  closed- loop  roots.  Later  the  time  response  option 
could  be  chosen.  This  is  because  the  subroutine  MEXP  (9),  called 
only  when  time  response  is  desired,  requires  a  relatively  large 
amount  of  central  processing  time.  All  options  are  initiated  by 


input  data.  A  guide  illustrating  the  form  and  purpose  of  all  the 
required  input  data  is  included  in  Appendix  B. 

As  the  program  executes  with  its  specified  options,  the 
output  contains  many  statements  analogous  to  prompts  in  an 
interactive  session.  This  simply  allows  the  user  to  follow 
the  execution  and  recollect  input  data  when  reading  the 
output.  Much  of  the  input  data  is  restated  in  the  output. 

Execution  begins  with  the  dimensioning  of  the  numerous 
indexed  variables  the  program  uses  as  well  as  initializing  of 
some  important  parameters  used  by  the  various  subroutines  in 
properly  handling  the  indexed  variables.  Thus  dimensioning 
need  not  be  altered  anywhere  except  in  the  main  program. 

Comment  statements  in  the  code  provide  guidance  on  specifying 
these  parameters  and  dimensions  before  compilation. 

As  a  final  note,  the  program  is  supported  by  several 
facility  routines  from  among  the  International  Mathematical 
and  Statistical  Library  <IMSL>  routines.  Their  code  is  not 
included  in  the  listing  in  Appendix  B. 


As  indicated  in  Section  I,  the  main  objective  of  the 
investigation  is  to  see  how  the  assignment  of  various  groups 
of  modes  can  affect  the  LOS  time  response,  and  see  if  any 
general  conclusions  might  be  warranted.  The  foundation  for 
such  an  investigation  is  baaed  upon  certain  facts  regarding 
the  basic  control  model  setup. 

As  set  up,  this  control  problem  recognizes,  as  is  often 
the  case,  that  there  is  a  limited  number  of  sensors/actuators 
that  can  be  practically  installed,  and  the  placement  of  these 
may  not  provide  all  the  observability  and  controllability 
desired.  For  a  multiple-controller  configuration,  however, 
the  described  block  triangularization  of  the  major  matrix 
results  in  the  various  controllers  commanding  varying  degrees 
of  controllability  and  observability.  To  be  sure,  there  will 
always  be  modes  which  are  not  directly  controlled  or  observed. 
A  large  space  structure  with  many  degrees  of  freedom  must  be 
modelled  with  a  greatly  reduced  number  of  modes  to  accommodate 
the  limited  capacity  and/or  speed  of  an  onboard  computer 
acting  as  the  controller.  Thus  it  is  of  interest  to  see  how 
the  choice  of  modes  to  control  as  well  as  the  assignment  of 
these  modes  in  a  multiple  control  scheme  affects  system  per¬ 
formance  . 

For  this  version  of  the  Draper-2  model  there  are  twenty- 
one  pairs  of  collocated  sensors  and  actuators  with  fixed 
placements.  Three  controllers  were  used  in  the  control  model. 


It  was  assumed  each  had  the  computational  capacity  for  four 
modes.  Thus  each  was  assigned  four  modes.  A  group  of  eight 
modes  was  considered  as  residual  to  represent  the  unmodelled 
modes  of  a  real  structure.  Thus  there  were  twenty  modes  in 
total  representing  a  truth  model  of  the  structure.  With  these 
figures  in  hand,  the  relative  degrees  of  controllability  and 
observability  of  the  individual  controllers  can  be  evaluated. 

Examining  Table  3-1  and  noting  this  version  of  the  simu¬ 
lation  program  performs  a  lower  block  triangular ization ,  a 
pattern  is  seen  as  to  the  amount  of  suppression  of  control  or 
observation  spillover  there  must  be  for  a  particular  control¬ 
ler.  To  suppress  observation  spillover  from  the  other  eight 
modes  assigned  to  Controllers  2  and  3,  Controller  (Ctlr)  1 
must  use  eight  of  its  sensor  inputs,  leaving  it  the  controller 
with  least  observability.  Ctlr  2  uses  only  four  of  its  sensor 
inputs  to  suppress  the  modes  of  Ctrl  3.  It  has  more  observ¬ 
ability.  In  contrast,  Ctlr  3  uses  eight  of  its  actuators  to 
suppress  control  spillover  from  the  other  eight  modes  in  Ctlrs 

1  and  2.  It  has  the  least  controllability.  Once  again,  Ctlr 

2  uses  only  four  actuators  in  suppression,  giving  it  more 
controllability  than  3.  Summarizing,  a  three-controller  con¬ 
figuration  with  lower  block  triangular ization  exhibits  one 
controller  (Ctlr  1)  commanding  the  most  controllability,  one 
(Ctlr  3)  commanding  the  moat  observability,  and  one  (Ctlr  2) 
commanding  a  median  amount  of  both.  Thus  the  model  setup  with 
three  controllers  allows  a  straightforward  study  involving 
assigning  the  various  groups  of  modes  permutati vely  to  the 
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three  controllers.  A  fourth  (residual)  group  provides  more 
realistic  insight  into  the  overall  stability  and  time 
response,  and  was  thus  included.  To  reiterate  the  underpin¬ 
ning  of  this  investigation,  stability  can  be  regarded  a  funda¬ 
mental  requirement  for  any  large  space  structure  control 
system,  but  time  response,  in  particular  LOS  pointing  here, 
may  be  the  driving  specification. 

The  investigation  began  with  a  search  for  candidate  modal 
rankings  which  could  be  divided  into  groups  for  assignment. 

One  way  to  rank  order  is  simply  to  put  the  mode  numbers  in 
ascending  order  (AO) .  This  might  be  the  situation  in  the 
early  stages  of  a  design  study  when  little  is  known  of  the 
relative  effects  of  each  mode  on  LOS  performance.  The  premise 
is  well  known  that  many  of  the  high  frequency  modes  can  be 
expected  to  lie  outside  the  control  bandwidth,  and  may  be 
simply  truncated  (removed  from  the  model).  An  excellent 
alternative  ranking  was  made  by  Varhola  (7) .  It  is  the  so- 
called  modal  cost  ranking  (MC) .  It  ranks  the  relative  effect 
of  a  mode  on  any  index  of  interest  against  the  other  modes. 
Similarly,  another  scheme  called  internal  balancing  (IB)  can  be 
used.  Varhola  generated  an  MC  ranking  for  LOS  performance  and 
an  IB  ranking  for  the  ability  of  the  sensors/actuators  to 
affect  a  particular  mode.  In  other  words,  the  second  ranking 
indicates  the  sensitivity  of  a  mode  to  control  and  control 
spillover,  whether  it  contributes  to  LOS  performance  or  not. 

The  AO  ranking  for  the  investigation  was  made  by  simply 
taking  the  three  rigid  body  rotation  modes  plus  the  next 
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seventeen  in  ascending  order  of  frequency.  The  MC  ranking  was 
made  by  selecting  the  twenty  highest  ranking  modes  among  those 
with  a  mode  number  falling  within  the  first  twenty  so  that 
they  would  match  with  the  AO  ranking.  The  modes  are  con¬ 
sistently  numbered  1  through  20,  with  the  rigid  body  transla¬ 
tion  mode  numbers  having  been  dropped.  The  rank  orderings 
obtained  are  as  follows: 

AO:  1,2, 3, 4  1  5, 6, 7, 8  I  9,10,11,12  I  13.14,15,16,17,18,19,20 

MC:  1,2,3,20  I  4,7,10,6  I  9,19,11,12  I  13,18,5,8,15,14,17,16 

IB:  1,2,3,7,4,8,5,10,20,6,13,9,12,11,19,18,15,14,17,16 

The  dividing  mark  between  the  groups  indicates  how  the  modes 
were  grouped.  Notice  the  first  group  (Group  1)  under  both  AO 
and  MC  analysis  contains  the  three  rigid  body  rotation  modes. 
The  fact  that  they  are  in  the  MC  Group  1  means  they  are  of  top 
importance  to  the  LOS  performance.  Notice  also  no  groupings 
are  presented  for  the  IB  ranking.  This  is  because  this  rank¬ 
ing  was  used  only  in  an  auxiliary  fashion  and  was  not  directly 
examined  in  the  investigation.  However,  to  illustrate  its 
usefulness,  notice  that  it  also  lists  modes  1,  2,  and  3  high¬ 
est.  This  is  good  news  for  effecting  control  as  it  means 
these  are  the  three  modes  most  sensitive  to  control. 

The  AO  and  MC  groups  were  assigned  to  various  controllers 
in  a  permutative  fashion.  For  each  ranking  there  were  six 
possible  combinations:  1-2-3,  1-3-2,  2-1-3,  2-3-1,  3-1-2,  and 
3-2-1.  wherein  the  ordering  implies  the  ordinal  number  of  the 
controllers.  For  example,  the  third  combination,  2-1-3, 


implies  that  the  modes  from  Group  2  were  assigned  to  Ctlr  1 


the  modes  from  Group  1  to  Ctlr  2,  and  the  modes  from  Group  3 
to  Ctlr  3. 

For  each  run  the  initial  conditions  were  such  that  the  LOS 
Radius  was  unity  at  time  zero.  (Note  that  no  units  of  measure 
need  apply  since  modal  coordinates  are  used.)  All  the  runs 
were  also  set  up  to  cover  20.0  seconds  of  propagation  with  a  t 
value  of  O.OS  seconds.  This  value  of  t  would  definitely  meet 
the  requirements  on  sampling  rate  to  reproduce  with  accuracy  a 
time  history  of  LOS  motion.  Some  preliminary  runs  were  made 
to  find  suitable  values  for  the  twenty  elements  of  the  Q  matrix 
diagonal.  (During  this  stage  of  the  investigation  an  error 
was  discovered  in  the  subroutine  FORHQl  algorithm  which  forms 
the  diagonal  Q  matrix.  The  algorithm  is  now  correct.)  The 
objective  was  to  apply  a  fixed  Q  matrix  such  that  all 
combinations  run  would  remain  stable  and  provide  meaningful 
data  on  time  response.  After  some  trial  and  error,  but  with 
much  insight  gained  into  how  each  mode  plays  into  the 
stability  and  movement  of  the  closed-loop  roots,  a  final  set 
of  Q  values  was  determined.  They  are  given  according  to 
ascending  mode  number  as  follows:  200  (x3),  500  (x7),  2500, 
3000,  500  (x6),  15000  (x2).  These  values  provide  for  a  stable 
suppressed  system  for  both  the  AO  and  MC  cases. 

With  the  control  model,  set  up,  six  runs  were  made  for  the 
AO  and  NC  cases.  In  all  twelve  runs  both  time  response  and 
eigenvalue  analysis  ware  made.  It  was  of  interest  to  see  how 
the  c.  ised-loop  suppressed  roots  shifted  from  their  original 
locations  where  they  had  the  given  damping  ratio  of  0.005. 


But  it  was  more  important  to  see  if  any  pattern  might  be 
revealed  by  the  time  response  plots  to  warrant  general  conclu¬ 
sion  about  the  time  response  performance  of  a  large  space 
structure  control -configured  as  described  when  various  modal 
assignments  can  be  made.  The  results  of  the  investigation 
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Figures  7-1  through  7-12  on  the  following  pages  are  plots 
of  the  closed- loop  roots  of  the  system  after  suppression.  For 
the  complex  roots  only  the  root  with  the  positive  imaginary 
part  is  plotted.  The  controlled  modes  have  roots  for  both  the 
controlled  states  and  the  observer  error  states.  Each  resid¬ 
ual  mode,  of  cburse,  has  only  a  single  root.  The  roots  were 
obtained  by  matching  them  with  the  identified  roots  of  the 
individual  (open-loop)  systems  as  generated  by  the  program. 
With  the  small  value  of  damping  ratio  of  0.005  before  closing 
the  loop,  the  mode  can  be  identified  by  recognizing  its  fre¬ 
quency  is  approximately  equal  to  the  imaginary  part  of  its 
root . 

In  all  cases  the  selected  values  for  the  0  matrix 
diagonal  resulted  not  only  in  a  stable  system  regardless  of 
modal  assignment,  but  also  resulted  in  significant  damping 
improvement  of  the  controlled  modes.  This  can  be  seen  in  the 
spread  of  the  roots  away  from  the  constant  damping  ratio  line 
of  0.005.  As  a  further  result  of  selecting  as  residual  the 
eight  lowest  ranking  modes  under  MC  analysis  and  the  eight 
highest  frequency  modes  under  AO  analysis  (frequency  trunca¬ 
tion),  the  residual  modes  remained  quite  steady.  Their 
damping  ratios  hardly  deviated  from  the  original  value  of 
0.005.  In  no  case  did  any  damping  ratio  decrease  signifi¬ 
cantly  below  0.005  whether  for  a  controlled,  error,  or  resid¬ 
ual  state.  To  be  noted  here  for  later  discussion  is  the  fact 
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that  the  root  for  mode  5  (imaginary  part  of  approximately 
1.65),  which  is  a  residual  mode  using  MC  ranking,  actually 
moved  left  under  all  six  combinations. 

Consistent  behavior  of  the  placement  of  the  roots  of  the 
of  the  error  states  relative  to  those  of  the  controlled  states 
for  a  particular  controller  is  exhibited.  The  controller  with 
most  observability,  Ctlr  3,  consistently  achieved  greater 
damping  of  the  error  states  than  the  controlled  states.  On 
the  other  hand,  Ctlr  1  consistently  achieved  greater  damping 
of  the  controlled  states  due  to  its  greater  degree  of 
controllability.  Ctlr  2.  however,  did  neither  in  a  consistent 
fashion.  Its  roots  shifted  so  drastically  after  closing  the 
loop  that  it  was  difficult  to  identify  or  distinguish  between 
the  controlled  and  error  states- -hence  the  notation  on  the 
figures . 

As  previously  mentioned,  under  MC  analysis  mode  5 
(residual)  was  made  more  stable.  In  contrast  notice  that  mode 
11  (controlled,  imaginary  part  approximately  7.69)  did  not 
move  away  from  the  original  damping  ratio  despite  a  relatively 
large  weighting  factor  of  2500.  This  was  true  in  both  MC  and 
AO  analysis.  The  explanation  of  these  two  phenomena  is  found 
in  the  IB  ranking.  Recall  this  ranking  gives  the  relative 
sensitivity  to  control  of  an  individual  mode  regardless  of  its 
importance  to  LOS.  Despite  the  fact  that  mode  5  was  of  little 
importance  to  LOS  and  was  therfore  made  residual,  it  was  still 
quite  sensitive  to  control.  Thus  spillover  moved  its  root 
significantly.  On  the  other  hand,  mode  11  may  have  been 


ranked  slightly  higher  than  mode  5  for  its  contribution  to  LOS 


but  it  was  relatively  insensitive  to  control.  Thus  its  root 
did  not  move  much. 

With  these  general  comparisons  made  among  the  eigenvalue 
plots,  the  telling  results  of  the  overall  system  LOS  per¬ 
formance  can  be  examined. 

Figures  7-13  through  7-24  on  the  following  pages  are  the 
time  response  plots  of  both  the  Radius  and  Defocus.  As  in 
previous  investigations  Defocus  was  found  to  be  of  little 
consequence  to  the  control  problem.  Its  time  dependent  ampli¬ 
tude  has  been  miniscule  compared  to  that  of  the  Radius. 
Therefore  it  will  not  be  discussed  further  here  since,  as 
mentioned,  all  cases  were  set  up  to  be  stable  and  its  time 
response  is  relatively  flat.  The  Radius'  time  response,  how¬ 
ever.  was  found  to  be  highly  dependent  on  modal  assignment. 

It  should  be  noted  that  in  some  cases  the  Radius  was  increas¬ 
ing  significantly  at  the  end  of  the  time  history.  Since  all 
cases  were  previously  shown  to  be  stable,  this  result  simply 
means  that  a  low  frequency  mode  (or  modes)  was  not  well  con¬ 
trolled  in  these  configurations  and  would  take  much  longer  to 
damp  out. 

Comparing  MC  to  AO  analysis,  the  various  combinations 
usually  exhibited  the  same  general  behavior.  When  those  modes 
important  to  control  because  of  their  relatively  high  contri¬ 
bution  to  LOS  were  assigned  to  the  controller  with  most  con¬ 
trollability.  Ctlr  1.  the  amplitudes  were  damped  out  in  a  well 
behaved  fashion.  However,  notice  that  nearly  equivalent  time 
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response  was  achieved  when  the  modes  of  Group  3  were  assigned 
to  Ctlr  2  as  seen  in  Figs  7-13.  14,  19,  and  20.  It  can  be 
seen  that  as  Group  1  is  shifted  to  Ctlr  2  and  then  Ctlr  3,  the 
time  response  deteriorates.  If  the  symbology  of  the  modal 
combinations  is  likened  to  a  cardinal  number,  one  might  be 
tempted  to  conclude  that  the  lower  the  numerical  value  of  the 
combination,  the  better  the  time  response.  Combination  (Comb) 
1-2-3,  analogous  to  the  number  123,  exhibited  very  good  over¬ 
all  time  response,  while  Comb  3-2-1  exhibited  poor  behavior. 
Notice,  however,  that  Comb  3-1-2  (Fig  7-17)  exhibited  better 
behavior  than  combination  2-1-3  (Fig  7-15).  Conversely.  Figs 
7-22  and  24  show  Combs  2-3-1  and  3-2-1  had  nearly  identical 
time  response.  (The  tabulated  values  differed  only  in  the 
third  decimal  place.)  Indeed,  the  beat  time  response  for  MC 


analysis  is  seen  in  Fig  7-14  to  be  Comb  1-3-2  and  not  1-2-3. 

Though  the  time  response  behavior  of  the  various 
combinations  bore  similar  comparison  between  the  MC  and  AO 


analyses,  it  was  not  the  case  that  MC  analysis  resulted  in 
significantly  better  time  response.  When  the  best  time 
response  cases  (Figs  7-13  and  19)  are  compared,  it  is  seen 


simple  frequency  truncation  (AO  analysis)  has  done  as  well  for 
this  model  as  the  effort  of  making  the  Modal  Cost  ranking,  if 
not  slightly  better.  Notice  the  rankings  given  in  Section  VI 
are  quite  similar  and  rank  the  three  rigid  body  modes  highest. 
The  relative  value  of  these  modes  was  so  large  that  Varhola 
could  not  scale  them  onto  his  ranking  charts.  In  summary  of 
the  comparison  of  AO  to  MC  time  response,  the  two  analyses  did 
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not  exhibit  a  clear  pattern  of  improvement  in  time  response  by 


combination.  In  order  of  improving  time  response  the  MC 
analysis  results  are  as  follows:  Combs  2-3-1,  3-2-1,  2-1-3, 
3-1-2,  1-2-3,  1-3-2.  But  for  AO  analysis  the  ordering  is 
different:  3-2-1,  2-3-1,  3-1-2,  2-1-3,  1-3-2,  and  1-2-3. 

With  the  above  results  in  mind  a  discussion  of  conclusions  to 
be  drawn  is  at  hand. 
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A  number  of  conclusions  can  be  drawn  from  the  results  of 
this  investigation. 

From  the  eigenvalue  plots  under  MC  analysis,  it  is  seen 
that  multiple  controllers  suppressing  spillover  can  effec¬ 
tively  control  a  high  frequency  mode  with  little  or  no  adverse 
effect  on  lower  frequency  modes.  For  example,  mode  20  was 
damped  without  driving  other  modes  unstable. 

However,  these  same  plots  show  that  mode  5,  a  residual 
mode,  was  also  affected'-by  spillover.  Though  it  was  damped 
in  this  investigation  under  conditions  of  these  gains,  it 
might  3ust  as  well  have  been  driven  unstable  by  the  spillover. 
It  would  be  better  to  be  able  to  predict  residual  modes  will 
maintain  their  original  damping  ratios  after  suppression  so  as 
to  not  be  of  concern.  Thus  it  should  be  concluded  that  Modal 
Cost  analysis  alone  is  a  risky  approach  to  take.  It  should  be 
accompanied  by  Internal  Balancing  analysis  to  reveal  problem 
areas.  These  arise  when  a  mode  has  a  high  IB  ranking  but  a 
low  enough  MC  ranking  to  leave  it  as  a  residual  mode. 

IB  analysis  can  also  reveal  candidate  modes  for  removal 
from  the  controlled  groups  in  exchange  for  just  such  problem 
modes.  For  example,  mode  11,  exhibiting  little  response  to 
control,  might  well  be  replaced  by  mode  5.  An  advantageously 
greater  range  of  gain  values  (Q  matrix  weighting  v  >lues)  might 
then  be  possible  because  it  could  be  expected  that  the 
residual  modes  would  remain  “steady.*'  It  is  now  obvious  that 


a  better  control  configuration  would  have  the  IB  ranking 
identical  to  the  MC  ranking.  The  designer  would  accomplish 
this  by  moving  sensors  and  actuators  to  increase  or  decrease 
their  effect  on  a  particular  mode,  resulting  in  changing  the 
IB  ranking.  An  even  greater  range  of  gain  could  result 
because  now  the  single  identical  ranking  would  simultaneously 
show  those  modes  most  contributory  to  LOS  also  being  those 
most  sensitive  to  control,  and  equally  important,  those 
contributing  little  to  LOS  being  relatively  insensitive. 

In  addition  to  the  drawbacks  of  applying  MC  analysis 
alone  as  mentioned  above,  another  conclusion  to  be  drawn  is 
that  MC  analysis  alone  does  not  necessarily  offer  significant 
time  response  improvement  over  AO  analysis,  which  is  the 
simplest  approach  to  take.  In  this  investigation,  AO  Comb 
1-2-3,  the  simplest  ordering  and  assignment  possible,  resulted 
in  the  best  time  response.  It  should  be  noted,  however,  that 
the  result  may  have  been  due  to  the  overwhelming  influence 
of  the  three  rigid  body  modes  on  LOS  performance. 

In  summary,  for  any  modal  assignment  approach  taken, 
unsuppressed  residual  modes  will  limit  the  amount  of  gain  that 
can  be  applied  to  improve  time  response.  Though  various 
analytical  rankings  can  be  made  and  used  together  to  effect 
significant  improvement  in  time  response,  this  investigation 
revealed  no  clear  pattern  among  the  analyses  to  predict 
behavior  based  on  modal  assignment.  The  results  and  conclu¬ 
sions  do,  however,  provide  good  ground  for  further  investiga¬ 
tions,  as  described  in  the  next  section. 


IX .  Recommendations 


Several  further  investigations  can  stem  from  this  one. 

As  mentioned,  mode  5  was  found  to  make  little  contribution  to 
LOS  while  being  quite  sensitive  to  control  spillover.  At  the 
same  time  mode  11  had  somewhat  opposite  characteristics. 

These  two  modes  are  good  candidates  to  lead  into  an  investiga¬ 
tion  of  the  effects  on  time  response  and  stability  gain  mar¬ 
gins  for  this  configuration  under  mode-swapping  conditions 
when  the  various  cases  of  modal  assignment  combination  are  run 
There  are  other  modes,  with  less  pronounced  characteristics, 
which  would  also  make  good  candidates  (e.g.,  modes  8  and  12> . 

Some  investigations  might  also  be  carried  out  on  the 
effects  of  new  placements  and  orientations  of  the  sensor/ 
actuator  pairs  based  on  the  modal  rankings  of  Ref  (7).  This 
investigation  could  include  the  addition  of  code  to  the  simu¬ 
lator  for  disturbance  rejection.  Moreover,  such  investiga¬ 
tions  could  include  studying  the  effects  of  ncise  in  the 
sensor  data. 

A  similar  investigation  to  this  one  begs  to  be  run  with 
an  increased  number  of  modes.  This  investigation  used  for 
continuity  of  analysis  a  twenty-mode  truth  model  to  search  for 
general  patterns  of  modal  assignment  effects.  More  revealing 
data  may  be  obtained  if  the  number  of  modes  assigned  to  each 
controller  could  be  increased  and  if  simultaneously  all 
significant  residual  modes  could  be  included  in  the 
suppression.  Doing  this  might  show  that  there  are  some  clear 


patterns  to  the  effects  of  various  modal  assignments  not 
revealed  in  this  investigation.  As  mentioned,  the  over¬ 
whelming  value  of  the  rigid  body  modes  in  the  rankings  may 
have  obscured  such  patterns. 

Along  these  lines  the  code  could  be  revised  to  be  more 
flexible  in  its  dimensioning.  Presently  the  program  must  be 
recompiled  if  the  quantity  of  modes  changes.  This  is  due  to 
the  requirement  that  certain  of  the  routines  (9)  have  exactly 
dimensioned  arrays,  as  opposed  to  the  IMSL  routines  which  do 
not  require  this.  There  was  not  time  in  this  investigation  to 
make  these  revisions,  but  they  would  greatly  enhance  the 
flexibility  of  the  program  and  the  speed  with  which  further 
investigations  could  be  carried  out. 


eigenvalue  analysis  summary  (INVERSE  POWER  METHOD) 


NUMBER  OF  EIGENVALUES  EXTRACTED  .  S4 

NUMBER  OF  STARTING  POINTS  USED .  10 

NUMBER  OF  STARTING  POINT  MOVES  .  0 

NUMBER  OF  TRIANGULAR  DECOMPOSITIONS  ....  61 

TOTAL  NUMBER  OF  VECTOR  ITERATIONS  .  673 

REASON  FOR  TERMINATION  .  7* 

LARGEST  OFF-DIAGONAL  MODAL  MASS  TERM . 20E-13 

...  SI 

MODE  PAIR  . 

...  37 

NUMBER  OF  OFF-DIAGONAL  MODAL  MASS 

TERMS  FAILING  CRITERION  .  0 

(•  1  OR  MORE  ROOT  OUTSIDE  FR. RANGE. 

SEE  NASTRAN  U.M.  PAGE  3.4-12) 
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50  MODES  AM)  frequencies 


real  eigenvalues 


MODE  EXTRACTION  EIGENVALUE 

NO .  ORDER 


RADIAN 
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CYCLIC  generalized  GE^CRALIZED 
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23 

0 

24 

s 

25 

1 

20 

2 

27 

3 

28 

4 

29 

7 

30 

10 

31 

17 

32 

24 

33 

28 

34 

30 

36 

32 

30 

34 

37 

38 

38 

30 

39 

39 

40 

41 

41 

40 

42 

42 

43 

43 

44 

44 

46 

46 

40 

48 

47 

47 

40 

40 

49 

50 

60 

49 

61 

61 

62 

62 

63 

53 

64 

64 

-2.489048E-06 
>l.a67722E>06 
-1.1U428E-06 
>9.82e476E-0e 
-9.6876636-00 
-2.0090916-09 
8.301744E-01 
2.736091E«00 
3.974420£«00 
4.374811£«00 
7.764729E«00 
1.318287E*01 
1.334010£«01 
6.912161E«01 
0.074248E«0t 
7 . 109790E»01 
1 . 109210E«02 
1.3068t3E*O2 
1 . 306862E#03 
1.400966E«02 
2.206301E«02 
3.628249E«02 
3.982044E«03 
4.629673E«02 
l.0618146«0a 
1 .092404€«03 
2.449482E«03 
2.900908E«03 
2.759142E«03 
2.899910E«03 
3.060609E«03 
3.060609E«03 
5.e83329E«03 
S.219047E«03 
9.393348E«03 
7.283877E«03 
7.426014E«03 
7.910140E«03 
9.071943E«03 
1.019817E«O4 
1 .070034E«04 
1.107000E«04 
1 . 161892E*04 
l.264812E*04 
1.4369816«04 
2 . 230987E«04 
2.363083E«04 
2.660074E*04 
3.743684E«04 
2.767676E«04 
3.872906E*04 
3.064867E«O4 
3.076400E«04 
3.213600E«04 


4 . 989037E-03 
3.0847276-03 
3. 3338006-03 
3.1347216-03 
3.0903776-03 
4.6716326-00 
9.1442676-01 
1. 08301 16«00 
1.993690E«00 
2.0916066«00 
2.784731E«00 
3.0a0822E«00 
3.063243E«00 
7.089061E«00 
8.109007E«00 
8.406874€«00 
1.081303E«01 
1 . 142732£«01 
t .  1427436*01 
1 . 1909946*01 
1.4860206«01 
1.878303E«01 
1 . 996606E*01 
2.128279€*01 
3.243108E«01 
3.3061636«0l 
4.949234E*01 
6.099907E«01 
6.262764E«01 
6.3860816*01 
6.637e01E*Ot 
6.6370016*01 
7.1290406*01 
7.2242976*01 
7 . 9968426*01 
8 . 6346036*01 
8.0108626*01 
8 . 8972096*01 
9.834004E«01 
1.0098006*02 
1.038284E«02 
l.062140E«02 
1. 07791 1E*02 
1.120184E«02 
1.190283e*O2 
1.4930166*02 
1.6339706*02 
1.8279066*02 
1.0604076*02 
1.0006966*02 
1 .0949066*02 
1.7470166*02 
1 . 753980E*02 
1 . 7920626*02 


7.941266E-04 
5.8O442OE-04 
6 . 3069206-04 
4.9990836-04 
4.9290396-04 
7.2768206-00 
1.466364E-01 
2.0321226-01 
3.172907E-01 
3.3288926-01 
4.4320386-01 
6.7780336-01 
6.8143176-01 
1.2237616*00 
1.3002336*00 
1.3476466*00 
1.7209476*00 
1.819099E*00 
1.8197326*00 
1.8991806*00 
2.3034926*00 
2.989609E*00 
3.1769406*00 
3.3872826*00 
6.1016026*90 
6.2003166*00 
7.8789346*00 
8.1167666*00 
8.3000186*00 
8.6706226*00 
8.8134016*00 
8.8134026*00 
1.1340236*01 
1.1497836*01 
1.2726786*01 
1.3683186*01 
1.3714166*01 
1.4100446*01 
1 . 6062206*01 
1.0072426*01 
1.0624816*01 
1.0746346*01 
1.7166486*01 
1.7028286*01 
1.9071206*01 
2.3771036*01 
2.4413996*01 
2.5988916*01 
2.0302636*01 
2.0429196*01 
2.0970206*01 
2.7817346*01 
2.7916656*01 
2.0630946*01 


*W 


MeaME.M 

eeveME.M 


-1.367722E.W 
-1.11142SE.06 
-9.S2S47SE-aS 
-S.6t7S63E-0e 
-2.eS9S91E-eS 
S.3S1744E-01 
2.736091E.00 
3.9744266.00 
4.374811E.00 
7.7B4729E.00 
1.3182S7E.01 
1.334019E.01 
S.9121S1E.01 
8.8742486.01 
7 . 1887906.01 
1 . 1892186.02 
1.3088136.02 
1 . 3888826.02 
1 . 4889886.02 
2. 2083016.02 
3.8282496.02 
3.9820446.02 
4.8298736.82 
1.8818146.03 
1.0924046.03 
2.4494826.03 
2.8009086.03 
2.7891426.03 
2.8999106.03 
3.0868696.03 
3. 0068696.03 
8.0823296.03 
8.2190476.03 
6 . 3933486.03 
7.2838776.03 
7 . 4280146.03 
7.9181406.03 
9.6719436.03 
1.8198176.04 
1.0780346.04 
1  .  1070806.04 
1 . 1818926.04 
1 .2848126 '.04 
1.4388816.04 
2.2308876.04 
2.3830836*04 
2.6800746*04 
2 . 7436846*04 
2.7878786*04 
2.8729086*04 
3.8848876*04 
3.0764886.04 
3.2136006*04 


Acossa  MODEL  eiqewata  duplication 
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flD-fll62  977  HOOfll.  ASSIGNMENT  EFFECTS  ON  DECENTRALIZED  CO»ITROLOFfl^2/2 
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AC0SS6  MODEL  EIGEWATA  DUELICATION 


august  13,  IMS  SELEASE  APR.  1M4 


S«  MOOES  AM>  PREOUENCIES 


S  0 

n  r  t 

D  SUL 

K 

DATA 

CARO 

COUMT 

1 

2 

% 

4 

6 

8  .. 

10»- 

CBAR 

118 

4M 

32 

33 

‘.6 

1.8 

104* 

C8AR 

117 

4M 

26 

34 

108- 

CSA0 

118 

4M 

26 

34 

108- 

CBAR 

118 

4M 

36 

34 

10T- 

CBAR 

130 

4M 

36 

36 

108- 

C8AR 

131 

4M 

31 

36 

108- 

C8AR 

133 

4m 

34 

36 

.6 

1.8 

110- 

CBAR 

133 

4m 

32 

36 

111- 

CBAR 

134 

4m 

33 

36 

113- 

CBAR 

138 

4M 

34 

3T 

113- 

CBAR 

138 

4m 

36 

36 

114- 

CBAR 

137 

3M 

26 

3T 

118- 

CBAR 

138 

3m 

26 

36 

118- 

CBAR 

138 

3m 

31 

36 

117- 

CBAR 

130 

3m 

31 

36 

118- 

CBAR 

131 

sm 

46 

46 

.6 

1.0 

118- 

CBAR 

133 

sm 

46 

66 

.6 

1.0 

IM- 

CBAR 

133 

sm 

66 

SI 

.6 

1.0 

121- 

CBAR 

134 

sm 

61 

62 

.6 

1.0 

123- 

CBAR 

138 

sm 

62 

43 

.6 

1.0 

123- 

CBAR 

138 

sm 

46 

53 

.6 

1.0 

124- 

CBAR 

137 

sm 

63 

64 

.6 

1.0 

128- 

CBAR 

138 

sm 

64 

$6 

.6 

1.0 

128- 

CBAR 

130 

sm 

66 

56 

.6 

1.0 

127- 

CBAR 

140 

sm 

66 

57 

.6 

1.0 

128- 

CBAR 

181 

4m 

6 

IS 

128- 

CELAS3 

142 

5.7ms 

4 

2 

42 

1 

130- 

CELAS3 

143 

s.yms 

4 

2 

42 

2 

131- 

CELAS3 

144 

s.ims 

4 

3 

42 

3 

133- 

CaAS3 

148 

s.yms 

3 

1 

46 

1 

133- 

CELAS3 

148 

s.yms 

3 

2 

46 

2 

134- 

CELAS3 

147 

s.yms 

3 

3 

46 

3 

138- 

CELAS3 

148 

s.yms 

6 

1 

47 

1 

138- 

CELAS3 

148 

s.yms 

6 

2 

47 

2 

ECHO 


10 


ri  e 
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ACOSSa  moCL  EIQEWATA  DUfLXCATION  AUGUST  13,  IStS  RELEASE  APR.  1M4 

sa  MOOU  AMD  nffiOUENCICS 

SORTED  aULK  DATA  ECHO 

CARO 


COUNT 

1 

2 

$ 

4 

..  6  .*  8 

7  . .  9  . .  9  .  .  10 

13T- 

CELAS2 

169 

6.79C3 

8 

3  47 

3 

13t> 

CONMJ 

691 

27 

378. 

13a- 

C0MM9 

692 

n 

378. 

14a- 

COMra 

693 

2% 

378. 

141- 

CONMS 

694 

378. 

14J- 

coNMa 

696 

32 

688. 

143- 

COMI2 

696 

33 

888. 

144- 

COMM 

697 

34 

288. 

14S- 

C0NM2 

699 

38 

288. 

14a- 

C0MM2 

699 

8 

688. 

147- 

C0NN2 

619 

18 

888. 

14a- 

C0NM2 

611 

U 

688. 

149- 

coae 

613 

13 

688. 

isa- 

COMME 

629 

14 

8.8 

lai- 

CONNS 

621 

18 

17*8 

tS2- 

CONNS 

632 

18 

17.8 

1S3- 

CONNS 

623 

17 

17.8 

1S4- 

CONNS 

624 

18 

17.8 

isa- 

CONNS 

62S 

18 

8.8 

isa- 

CONNS 

644 

44 

3688. 

»644 

117- 

*644 

2.19E3 

2.18C3 

4.28E3 

ita- 

CONNS 

649 

48 

98.88 

♦848 

ita- 

*649 

379.9 

laa- 

CONNS 

669 

68 

188*8 

♦868 

lai- 

*669 

649.9 

lai- 

CONNS 

662 

62 

98.8 

♦662 

ia3- 

*662 

379.9 

ia4- 

CONNS 

663 

63 

98.8 

♦683 

lea- 

*663 

379.9 

laa- 

CONNS 

666 

88 

188.8 

♦688 

167- 

*666 

649.9 

199- 

CONNS 

667 

67 

98.8 

#667 

199- 

*667 

379.9 

179- 

CRIGOl 

141 

44 

42 

43  48 

48  47 

I 

t 


c 
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Aco$s6  Mooe.  ezccmmta  oun.K*noN 


AUGUST  13,  IMS  RELEASE  ARR.  1M4 


S«  MOOES  AM>  RREQUBICIES 


sorted  SULK  DATA 


CARO 


COIWT 

1 

2 

..  3  ••  4 

6 

8 

171- 

EIOR 

000 

XNV  .# 

i0.0 

SO 

17J- 

•  10 

mass 

173- 

OHIO 

1 

-r.tt 

.0 

.0 

174- 

QRIO 

2 

S.0 

.0 

17t. 

QRtO 

3 

•4.« 

•S.0 

.0 

170- 

GRIP 

4 

.0 

S.0 

.0 

177- 

Goto 

0 

4.« 

S.0 

.0 

170- 

QRtO 

0 

4.0 

•S.0 

.0 

170- 

GRIO 

7 

7.0 

0 

.0 

100- 

GRIO 

0 

•7.0 

.0 

2.0 

101- 

GRIO 

0 

•4.0 

S.0 

2.0 

102- 

GRID 

10 

•4.0 

•S.0 

2.0 

101- 

GRIO 

11 

4.0 

S.0 

2.0 

104- 

GRIO 

12 

4.0 

•S.0 

2.0 

lOt- 

GRIO 

11 

7.0 

.0 

2.0 

IM- 

GRID 

14 

•0.0 

.0 

13. 

107. 

GRID 

IS 

•4.0 

4.0 

12. 

100- 

GRID 

10 

•4.0 

•4.0 

12. 

100- 

GRID 

17 

4.0 

4.0 

12. 

lOO- 

GRID 

IS 

4.0 

•4.0 

12. 

101- 

GRID 

IS 

0.0 

.0 

13.0 

102- 

GRIO 

20 

•S.0 

.0 

23.0 

102- 

GRID 

27 

•4.0 

1.0 

23.0 

104- 

GRID 

20 

•4.0 

•S.0 

23.0 

IM- 

GRIO 

20 

4.0 

S.0 

23.0 

IM- 

GRID 

30 

4.0 

•s.0 

22.0 

107- 

GRIO 

31 

S.0 

.0 

22.0 

IM- 

GRID 

12 

•4.0 

10.0 

23.0 

100- 

GRID 

IS 

4.0 

\0.0 

22.0 

200- 

GRIO 

34 

•4.0 

•10.0 

22.0 

201- 

GRID 

3S 

4.0 

•10.0 

22.0 

303- 

GRID 

30 

•4.0 

S.0 

24.0 

203- 

GRID 

17 

•4.0 

•s.0 

24.0 

204- 

GRID 

30 

4.0 

s.0 

24.0 

ECHO 

I  . .  8  . .  10 

1.0  E-03«10 
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ACQSM  MODEL  EIGEWATA  DUPLICATION 


AUGUST  13,  IMS  SELEASE  APS.  1M4 


$•  MOOES  AW  FREQUENCIES 


SORTED  BULK  DATA  ECHO 


COUNT 

1 

2 

1 

4 

1 

6 

..  T 

MS* 

«» 

12 

4.2 

-1.2 

24.2 

2M* 

GRID 

42 

.2 

.2 

2.6 

MT- 

GRID 

41 

.2 

-1.2 

1.2 

3M- 

GRID 

42 

.2 

1.2 

.2 

2«t- 

GRID 

41 

-2.6 

.2 

.2 

2ia- 

GRID 

44 

.2 

.2 

.2 

211- 

OHIO 

41 

2.2 

.2 

.2 

212- 

GRID 

46 

-4.2 

-1.2 

.2 

211- 

GRID 

47 

4.2 

-1.2 

.2 

214- 

GRID 

46 

-M.6 

.2 

.2 

21S- 

GRID 

42 

-21.22 

.2 

.2 

216- 

GRID 

SR 

-16.6 

.2 

.2 

21T- 

GRID 

fl 

-11.2 

.2 

.2 

216- 

GRID 

12 

-6.2 

.2 

.2 

212- 

GRID 

11 

6.6 

.2 

.2 

226- 

GRID 

14 

11.6 

.2 

.2 

221- 

GRID 

ss 

16.6 

.2 

.2 

222- 

GRID 

M 

21.22 

.2 

.2 

223- 

GRID 

17 

M.6 

.2 

.2 

224- 

MATl 

IM 

1.24C-tt 

.1 

22S- 

OMITl 

SB 

46 

S6 

12 

11 

11  17 

2M- 

OHITl 

4Ba 

2 

16 

11 

12 

14  IS 

22T- 

ONITI 

4H 

XI 

16 

19 

2T 

M  29 

3M- 

OHITl 

4S6 

32 

TM2U 

11 

222- 

OMIT! 

1234SS 

1 

THIU 

• 

236- 

OMITl 

12S4S6 

13 

M 

31 

49 

SI  S4 

231- 

OHITl 

123466 

36 

THiU 

41 

232- 

PARAM 

GRDPNT 

6 

233- 

PBAR 

266 

166 

6 .2B6E-41 . 6BSE-63 .66SE-66 . 1B6E-6 ' 

234. 

PBAR 

366 

166 

3 . 133E-41 .  S6SE-61 .  S6BE-61 . 1 16E-6 

23S- 

PBAR 

422 

166 

3 . 61SE-43 .64SC-63 . 646E..66 . 6fSE-6 

216- 

PBM 

122 

166 

6 . 46TE-41 . 6T4E-S 1 .  tr4E 

.SS.74SE-S 

ENOOATA 

u 


Pr«»aDDed  product  of  ••n«or/ actuator  ■appinq  ■atricea  and  ■odal  natrix. 


(Each  block  balow  is  a  row  in  the  product  matrix.  Tha  nuabar  of  columns 
aquala  tha  nuabar  of  sansora/ actuators.  Tha  number  of  rows  equals  tha  nuabar 
of  modes  in  tha  nodal.  Kara  fifty-four  nodes  are  given.) 


N.B.  THE  FIRST  SIX  BLOCKS  (ROUS)  ARE  FOR  THE  SIX  RIGID  BODY  HODES. 


-.315269E-02 

.683988E-02 

-.466864E-02 

.158803E-02 

-.554972E-03 

.517895E-02 

.274804E-03 

-.142412E-01 

.426474E-03 

-.138475E-02 
- .827424E-02 
-.300683E-02 
-.213272E-02 

-.848535E-03 

-.142413E-01 

.386633E-02 

.272722E-02 

-.111276E-02 

.101214E-01 

.S32660E-02 

-.147629E-02 

.462729E-02 

-.4338S6E-02 

-.7S6388E-02 

-.818312E-02 

-.715295E-02 

-.3SS801E-03 

-.914393E-02 

-.879917E-02 

.4S2030E-04 

-.316205E-02 
. 14771 lE-01 
-.113472E-01 
.271932E-02 

.24300SE-02 

-.879919E-02 

-.131719E-02 

-.473004E-02 

.811893E-03 
-.193392E-02 
-.329975E-02 
. 32001 lE-02 

-.498062E-02 

-.405760E-02 

-.194591E-02 

-.402228E-02 

-.299068E-02 

-.20120SE-02 

-.401340E-02 

.S2470SE-02 

.160SS0E-02 

-.196797E-02 

-.393003E-02 

-.399141E-02 

.991013E-04 

-.471914E-02 

.524703E-02 

-.200323E-02 

-.301276E-02 

.513882E-02 

-.407961E-02 

-.404874E-02 

.946879E-04 

-.474338E-02 

.332965E-02 

.114881E-01 

.S43203E-02 

.806421E-02 

.754460E-02 

.595678E-02 

.S94440E-03 

-.1057S0E-03 

.101718E-01 

.103287E-01 

.727127E-02 

-.470226E-02 

-. 13721 lE-03 
.394416E-03 
.807005E-02 
.673062E-02 

-.466108E-02 

.201362E-02 

.383313E-02 

-.496301E-02 

- . 160182E-02 
.396H1E-03 
-.128230E-02 
.131966E-02 
-.154559E-03 

-.301424E-02 
.155050E-02 
- . 788180E-03 
.111043E-01 

-.183961E-02 

.444494E-02 

.212784E-02 

-.234843E-02 

. 120376E-01 
-.786007E-03 
-.278332E-02 
-.731816E-03 

-.309307E-02 
- . 181398E-03 
.626848E-03 
-.248390E-02 

-.923363E-02 
-.667933E-02 
- . 102700E-02 
-.381S38E-02 
-.295818E-02 

-.639696E-02 

-.309936E-02 

-.315849E-02 

-.268068E-02 

-.281699E-02 

-.618778E-03 

-.130938E-02 

-.977603E-02 

-.2002&8E-02 

-.313849E-02 

-.368097E-02 

-.474816E-02 

-.103184E-01 

-.846932E-02 

-.396334E-02 

-.101340E-01 

-.675074E-04 
-.152481E-02 
. 150731E-02 
-.151452E-02 
.258268E-06 

.132337E-02 
-.151058E-02 
.838672E-04 
- . 18829SE-02 

.131036E-02 

.388829E-02 

-.130637E-02 

-.172713E-03 

-.196338E-02 

.847191E-04 

.132103E-02 

-.332031E-06 

.839319E-04 

-.132630E-02 

-.132188E-02 

-.131139E-03 

A-11 
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461975E-02 

935799E-03 

693351E-03 

806480E-03 

111913E-05 

283244E-03 

S14210E-02 

519414E-02 

517286E-02 

517541E-02 

493235E-02 

783624E-03 

282S48E-02 

103421E-02 

148818E-02 

279027E-03 

779361E-02 

792822E-02 

787447E-02 

787970E-02 

1SS624E-01 

260818E-02 

812426E-02 

312369E-02 

4S4232E-02 

380725E-03 

346788E-02 

360952E-02 

351878E-02 

140472E-04 

606722E-04 

814478E-04 

182670E-03 

138635E-03 

895722E-06 

104122E-02 

S88150E*03 

110052E-02 

303508E-03 

493212E-05 

130713E-02 

S2206SE-03 

205388E>02 

1313S8E-02 

144910E-02 


.935707E-03 

.791320E-03 

.494942E-02 

.209309E-02 


.514215E-02 

.51758SE-02 

.330027E-03 

•.S36401E-07 


.783646E-03 

.148747E-02 

.389S43E-03 

.434250E-07 


.779379E-02 

.788090E-02 

.490472E-03 

.189S4SE-06 


■.262428E-02 

.4S4S19E>02 

•.124810E-02 

.321106E-04 


.345091E-02 
•.351501E-02 
■.600439E-03 
■ , 175764E-01 


.813332E-04 
. 147368E-03 
.S94099E-04 
.569513E-03 


.591641E-03 

.407276E-03 

.124178E-02 

.148712E-02 


.522738E-03 
. 14S579E-02 
.7S4198E-03 
. 129328E-08 


.791995E-03 
■ . 102444E-01 
■.690607E-03 
•.180235E-05 


.S17385E-02 

.334293E-07 

.S19422E-02 

.262004E-03 


.148636E-02 

•.472337E-06 

.282720E-02 

■.493367E-02 


.788088E-02 

.67749SE-07 

.7928S3E-02 

.275233E-03 


.4S1S60E-02 

■.908189E-0S 

.8i6049E-02 

•.155719E-01 


.3S4388E-02 

.152819E-02 

.355857E-02 

.493670E-04 


.147687E-03 

.201893E-03 

.182640E-03 

.223572E-06 


.410S06E-03 

.28S6S2E-02 

■.1106S6E-02 

.945050E-06 


.145S87E-02 

•.151901E-0S 

.20S523E-02 

■.130303E-02 
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■.784817E-02 

.494835E-02 

.916870E-03 

.732430E-06 


•.239311E-05 

.329817E-03 

.515428E-02 

.515777E-02 


.108730E-0S 

.389843E-03 

-.331138E-03 

.351962E-03 


.718050E-05 

.449849E-03 

.782726E-02 

.783367E-02 


.S94745E-04 

.124903E-02 

-.121100E-02 

.960752E-03 


• . 170481E-01 
■.599366E-03 
.347343E-02 
.269304E-05 


.485137E-03 

.623802E-04 

.926131E-04 

.777539E-07 


.290310E-02 

.123112E-02 

■.4&1495E-03 

.427879E-05 


.  126968E-04 
.7S2184E-03 
.685619E-03 
.991604E-03 


.4&1604E-02 
.  1043S8E-02 
.917982E-03 
.646094E-06 


.263218E-03 

.513579E-02 

.S1S426E-02 

.2S9302E-03 


■.49347QE-02 

.1924S0E-02 

■.330840E~03 

•.515918E-02 


.278830E-03 

.778169E-02 

.782722E-02 

.2&8642E-03 


. 155789E-01 
■  .623355E>02 
■.119006E>02 
•.162719E-01 


.479175E-03 

.344473E-02 

.348078E-02 

.522509E-04 


.606673E-04 

.614378E-04 

.936193E-04 

.215776E-05 


.104141E-02 
. 135586E-02 
.467948E-03 
.H3570E-04 


.130948E-02 
■.486307E-04 
.88S142E-03 
. 140079E-02 


• 

-.100571E-01 

.624191E-02 

.397951E-02 

.872781E-03 

-.342586E-04 

-.62S689E-02 

.163093E-03 

-.368026E-02 

-.180601E-02 

-.1422a0E-03 

.954304E-02 

-.402034E-02 

-.139329E-05 

-.1504aiE-01 

-.359604E-02 

-.553053E-02 

-.250137E-04 

. 100722E-01 
.llia30E>01 
.S5a021E-02 
-.716567E-04 

# 

-.403904E-04 

-.823301E-04 

-.671748E~04 

-.754719E-04 

-.701908E-04 

-.838498E-04 

-.706823E-04 

-.191900E-04 

-.374604E-06 

-.712110E-04 

.798046E-06 

-.672465E-04 

-.393685E-04 

-.547136E-06 

-.190752E-04 

-.a362a2E-04 

-.7a3ai6E-04 

-.39113aE-04 

-.922450E-04 

-.a23625E-04 

-.402739E-04 

• 

-.194436E-03 

-.352121E-05 

.13191SE-03 

.665097E-04 

.821806E-04 

-.102196E-04 

.ai8227E-04 

-.111890E-04 

-.16626aE-05 

.800289E-04 

.442a90E-05 

.130320E-03 

-.190aa6E-03 

-.7526iaE-05 

-.112a22E-04 

.a21622E>05 

.3725iaE-04 

- . iaa437E-03 
-.512271E-04 
. 1254a4E-04 
-.197742E-03 

r 

.790560E-08 

.230921E-07 

-.268ailE-07 

.27266aE-07 

.207441E-08 

-.183191E-07 

.262764E-07 

.344111E-08 

-.129133E-07 

-.222229E-07 

.106553E-0a 

.307527E-07 

.112725E-0a 

-.3a7317E-0a 

.3a7022E-oa 

-.192742E-07 

.2197iaE-0a 

-.5a6197E-0a 

.227e97E-07 

.24132aE'07 

.i2a72aE-oa 

# 

-.3S9137E-02 
-.110571E-01 
. 15351 lE-01 
-.133859E-01 
-.337446E-04 

.110446E-01 

-.129291E-01 

-.154157E-02 

.889068E-02 

.129350E-01 

>.5a3540E-03 

-.153706E-01 

.297970E-05 

.207aoaE>02 

-.iaa542E-02 

.115545E-01 

-.560106E-05 

.3ai521E>02 

-.107320E'01 

-.116099E-01 

-.6a4433E-04 

• 

.406412E<03 

-.227008E-03 

-.706913E-03 

-.435593E-03 

-.379977E-03 

-.226147E-03 

-.391974E-03 

-.269481E-03 

-.4d3766E-06 

-.393053E-03 
. 141276E-05 
-.70644aE-03 
.39617aE-03 

-.204132E>04 

-.266115E-03 

-.215413E-03 

-.313a69E-03 

.405306E'03 

.422aaaE'04 

-.216107E~03 

.422012E'03 

.909885E-03 

-.712788E-04 

-.896288E-03 

.325626E-03 

.553053E-05 

.716368E-04 

.354204E-03 

.484372E-04 

-.491253E-03 

-.355725E-03 

-.696177E-04 

.701136E-03 

-.456432E-06 

.6a9764E-03 

.45a267E-04 

.2617a2E-04 

.152662E-05 

-.914633E'03 

-.429031E'03 

-.275099E'04 

.101552E'04 

-.118932E-03 
- . 380480E-03 
.507984E-03 
-.449199E-03 
-.438344E-06 

.379430E-03 
-.391738E-03 
-.420348E-04 
. 102954E-03 

.391920E-03 

.472439E-04 

-.50a923E-03 

.270960E-07 

-.523aa0E-04 

-.44ai42E-04 

.39a334E-03 

-.31164aE-06 

.1193a9E~03 

-.391549E>03 

-.39a301E'03 

-.127353E>05 

e 

-.555878E-10 

.589546E-10 

.120012E-09 

.674301E-10 

.178712E-10 

.5d5491E-10 

.232612E-10 

.391314E-10 

.229694E-12 

.2377aaE-10 

-.aa3301E-12 

.119a66E-09 

-.511310E-10 

.909964E-11 

.374241E-10 

.2a2835E-10 

.4239a3E-10 

-.552615E>10 

-.1966a9E>10 

.285aa3E-10 

-.517914E'10 
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.334063E-04 
.3S5563E-04 
.716573E-04 
.40S964E-04 
.  107336E-04 

.271519E-08 

.285738E-08 

.594557E-08 

.328203E-08 

.874981E-09 

.151482E-02 

.733758E-02 

.180890E-01 

.492848E-02 

.34S359E-04 

.471273E-04 

.193540E-03 

.2440S2E-03 

.814419E-04 

.174910E-04 

,603767E-02 
,111844E-01 
,  143947E-01 
,33675lE-02 
196951E-02 

231762E-04 

375830E-04 

432253E-04 

108560E-04 

402984E-0S 

254078E-05 

S26842E-0S 

123858E-04 

651099E-05 

833832E-06 

556024E-04 

516701E-04 

578764E-05 

dSld84E-0S 

403731E-05 

537188E-03 

311124E-02 

955181E-03 

1S1061E-02 

244244E-02 


.349799E-04 
. 140615E-04 
.235031E-04 
.120033E-06 


-.289700E-08 

-.112267E-08 

-.191492E-08 

-.148776E-10 


>.788804E-02 

.4191S8E-02 

.2a0737E-03 

-.894867E-03 


.18S191E-03 

-.148923E-04 

-.238S94E-04 

-.319109E-08 


-.108939E-01 

-.132992E-02 

.363771E-02 

.431368E-04 


'.361894E-04 

-.394272E-03 

.104030E-04 

-.240884E-06 


'.579623E-03 

-.336410E-06 

.739684E-07 

.217392E-06 


-.480963E-04 
- . 132390E-04 
-.621297E-03 
.291848E-04 


.320830E-03 
. 198102E-02 
.311413E-02 
.704393E-03 


.141842E-04 

-.439932E-06 

.723634E-04 

-.307033E-04 


-.118124E-08 

.393833E-10 

-.379802E-08 

.230274E-08 


-.434479E-02 

.416830E-02 

.179648E-01 

-.119637E-03 


-.163577E-04 

.124739E-05 

.282233E-03 

.382892E-04 


-.130046E-02 

-.190304E-03 

-.131189E-01 

-.640327E-02 


-.400639E-03 

-.676113E-06 

-.437333E-04 

-.268386E-04 


-.6S3088E~06 
. 131493E-06 
-.117262E-04 
.320805E-05 


.12288SE-04 

-.864821E-03 

.734617E-03 

.260441E-06 


.233603E-02 

.114088E-01 

-.723271E-03 

.204818E-04 


.53816SE-0S 

.224606E-04 

.170467E-04 

.234518E-04 


-.461587E-09 

-.183449E-08 

-.137370E-08 

-.207669E-08 


-.400377E-02 

-.332042E-03 

.283214E-02 

-.234803E-03 


.808940E-03 

-.327388E-04 

-.189090E-04 

.783117E-04 


.490010E-04 

.418424E-02 

.260340E-02 

-.538349E-02 


.273338E-06 
. 122041E-04 
.588253E-03 
-.173371E-04 


.887237E-06 
.436827E-06 
26341 lE-03 
-.238316E-05 


-.434313E-04 

-.564735E-05 

-.267803E-03 

-.66S818E-06 


.301840E-03 
. 138503E-02 
.987384E-03 
.843934E-03 


-.331338E-04 
- . 122694E-04 
.171069E-04 
-.310967E-04 


.271067E-08 

.879304E-09 

-.141033E-08 

.233372E-08 


.128023E-02 

-.212311E-01 

-.276292E-02 

.390807E-04 


.474279E-04 

-.273739E-03 

-.230104E-04 

.768187E-04 


-.397131E-02 

.207360E-01 

.285659E-02 

-.779453E-02 


-.236506E-04 

.602995E-04 

.668737E-05 

-.299465E-04 


.298614E-03 

.236420E-0S 

-.274132E-03 

.173710E-05 


.586904E-04 

-.121984E-04 

.130638E-05 

-.357343E-05 


.343037E-03 

.206397E-02 

.877781E-03 

-.237363E-03 


.226417E-03 
.429029E-02 
.522245E-02 
.  106591E-01 
.  187453E-01 

.637010E-04 

.25S410E-03 

.410383E-03 

.733537E-03 

.464027E-04 

.119527E-02 
.161757E-02 
.463689E-02 
.  104078E-01 
.  109491E-02 

.736339E-03 

.111804E-01 

.S40285E-02 

.279214E-02 

.40S806E-02 

.199628E-02 

.961143E<02 

.101007E-01 

.193304E-02 

.829016E-02 

.164774E-02 

.579112E-02 

.199181E-01 

.271237E-03 

.113763E-04 

.177836E-02 
.427994E-02 
.  169200E-01 
.874318E-02 
.738973E-02 

.240188E-02 

.291031E-02 

.897993E-02 

.927309E-02 

.761999E-02 

.  184702E-02 
.430191E-02 
.262082E-01 
.699111E-02 
.982496E-03 


-.632991E-02 

.161342E-01 

-.397380E-02 

-.229962E-03 


-.283002E-03 

-.990240E-03 

.423990E-04 

.280098E-03 


-.297790E-02 

-.114949E-01 

.174024E-02 

.398177E-02 


.909146E-02 

.107226E-01 

.113787E-01 

-.894186E-03 


.109148E-01 

-.110949E-01 

.320126E-02 

-.266018E-03 


.749923E-02 

.101918E-01 

.603894E>03 

.360392E-02 


.128902E-02 

.962841E'02 

-.448668E-03 

-.133946E-03 


-.748290E-02 

.904949E-02 

-.368983E-02 

-.139693E-03 


- . 194987E-02 
. 109187E-01 
'.198301E-02 
-.707292E-03 


.160226E-01 

-.191693E-02 

-.606902E-02 

.981327E-03 


. 107939E-02 
-.436772E-04 
.349771E-03 
.126739E-04 


.141724E-01 

-.249733E-03 

.343649E-02 

.249900E-03 


-.267499E-02 

.262808E-02 

.221183E-02 

.892236E-03 


-.107018E-01 

.193793E-03 

.194420E-01 

-.178092E-06 


-.10g676E-01 
. 248202E-03 
-.195062E-01 
.492917E-03 


.877239E-02 

.676491E-03 

-.776167E-02 

.242971E-02 


.796737E-02 

-.548391E-03 

-.193782E-01 

.242902E-02 


-.119169E-01 

.112967E-03 

-.272186E-01 

-.904793E-03 


.693207E-03 

-.299728E-02 

-.281601E-02 

.852733E-02 


. 142278E-04 
. 48981 3E-04 
-.219148E-03 
.334732E-04 


.111884E-03 

.211879E-02 

-.296028E-02 

.997278E-03 


-.113467E-02 
. 130227E-01 
.188029E-03 
.906370E-03 


.882091E-03 

.199382E-02 

-.297297E-02 

.232023E-02 


.434937E-02 

.248690E-02 

-.120381E-09 

.716848E-03 


.279248E-03 

-.324161E-02 

.107322E-01 

.349981E-02 


.498629E-03 

-.420382E-03 

.241906E-02 

.267994E-02 


-.672602E-03 

-.222609E-02 

-.917448E-02 

.494994E-03 


40331 lE-03 
.195110E-02 
-.296g04E-02 
-.993008E-06 


.798894E-04 

-.219371E-03 

.179096E-03 

-.290639E-04 


. 196016E-02 
-.362017E-02 
.227948E-02 
-.319996E-03 


.217971E-02 

-.196309E-02 

-.939477E-03 

.197179E-03 


-.149119E-02 

-.367990E-02 

-.963129E-03 

.287037E-02 


.170984E-02 

.220379E-02 

.274313E-03 

.966144E-03 


.192959E-02 

-.763846E-02 

.392012E-02 

.164313E-02 


.209999E-02 

-.147144E-01 

.920693E-02 

.930994E-03 


.962223E-03 

-.124927E-01 

.862136E-02 

-.347789E-03 
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552044E-02 

.228056E-01 

-.208466E-02 

-.401412E-03 

-.509181E-02 

186123E-01 

- . 137768E-02 

.678153E-04 

.528961E-02 

-.143599E-02 

290120E-01 

.280988E-02 

-.240259E-01 

.235477E-02 

.105350E-02 

660742E-02 

426498E-03 

-.955036E-04 

-.640794E-02 

.958732E-02 

-.466349E-02 

192S88E-02 

.209248E-01 

-.807090E-02 

-.288086E-02 

.252983E-02 

242720E-01 

.909374E-02 

. 109431E-02 

-.824936E-02 

-.166072E-02 

129665E-01 

-.638088E-02 

.146887E-01 

-.262354E-02 

.313816E-02 

784463E-02 

967373E-03 

-.S31398E-03 

.341206E-03 

-.636178E-03 

.952803E-03 

103436E-02 

-.406277E-02 

-.566725E-02 

-.650086E-03 

- . 190237E-02 

828S80E-02 

-.124772E-01 

-.281084E-03 

-.58843SE-02 

-.455617E-02 

252619E-03 

.781820E-02 

.149419E-02 

.148808E-01 

.3S4281E-01 

550145E-02 

863932E-02 

. 180728E-04 

-.175386E-02 

-.704243E-02 

-.484477E-03 

154460E-02 

.835839E-02 

.149743E-01 

.130651E-02 

.345791E-03 

284185E-02 

-.443142E-02 

-.6H294E-02 

-.426798E-02 

-.137994E-02 

1843S3E-02 

.429431E-02 

.102076E-02 

-.338997E-01 

. 127035E-01 

281911E-02 

504835E-02 

-.5S4804E-04 

.8S81S6E-03 

.349911E-02 

.164146E-03 

142008E-03 

.224142E-02 

.472764E-02 

.889370E-03 

.114589E>03 

378594E-03 

-.323057E-02 

.296670E-01 

.144188E-02 

-.349620E-03 

228201E-02 

.492052E-02 

-.916499E-03 

-.698678E-02 

.513412E-02 

416381E-03 

8228S3E-03 

-.170479E-03 

-.901519E-06 

.912994E-03 

- . 140085E-03 

124063E-02 

. 174579E-02 

.779892E-02 

.701055E-03 

.253186E-02 

318422E-02 

.780943E-02 

. 168804E-02 

-.380284E-02 

-.107403E-02 

127346E-01 

-.498741E-03 

.757852E-02 

.995247E-02 

-.305264E-03 

158982E-01 

648640E-02 

-.383383E-03 

.995313E-03 

.384417E-02 

-.905353E-03 

114790E-07 

. 160498E-06 

.982832E-07 

-.449983E-08 

.632603E-07 

876380E-07 

.203879E-06 

.359133E-07 

- . 121692E-06 

-.303988E-07 

162093E-06 

-.480372E-07 

.233667E-06 

. 172712E-06 

.133042E-07 

252357E-06 

128938E-06 

- . 104950E-07 

.207673E-07 

.690105E-07 

-.114697E-07 

605729E-02 

.219745E-02 

.486531E-02 

-.416147E-02 

-.842852E-02 

529524E-02 

-.124708E-01 

-.499496E-02 

.685893E-02 

.206975E-02 

830728E-02 

-.413295E-02 

-.111182E-01 

.384558E-02 

-.821254E-03 

256678E-01 

188288E-02 

-.281944E-02 

-.802565E-03 

-.213810E-02 

.154408E-02 

133540E-01 

-.118397E-02 

.397182E-02 

.115957E-01 

.188347E-01 

140901E-03 

-.618315E-02 

-.477920E-02 

.338930E-02 

-.469804E-03 

158237E-03 

-.346920E-02 

.129432E-02 

.180813E-02 

.116392E-02 

516107E-02 

209307E-02 

-.143190E-02 

-.290321E-03 

-.281857E-03 

-.202135E-02 

419823E-02 

37a098E-03 

531960E-02 

144181E-01 

333217E-02 


.300807E-02 

.2a0769E-02 

•.895417E-02 

•.144233E-04 


.455697E>02 

•.102dl8E-01 

.537789E-02 

.2126a2E-03 


.272211E-02 
.869922E-02 
.165158E-03 
.  179038E-02 


.306943E-02 

.204967E-02 

.789054E-02 

.361545E-03 


100376E-02 

S48498E-03 

32a094E-02 

767619E-02 

295946E-02 


.308717E-02 
.2225&1E-02 
.S41218E-02 
. 183726E-02 


.784317E-02 

.365427E-02 

•.358549E-02 

.191140E-03 


.956176E-03 
.376257E-02 
.769268E-03 
.  146357E-02 


.688409E-03 

.979819E-03 

.389560E-02 

.137414E-03 
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Card 

or 

Cards  <  C  ]  ) 

Range 

of 

Values 

Remarks 

DEC 

3,4 

Number  of  Controllers 

NC1,NC2,NC3,NR 

MACT,NSEN,ZETA 

N/A 

Numbers  of  modes  asgn'd. 

Number  sensors,  actuators,  zeta 

CPHIA] 

N/A 

Premapped  modal -sensor  matrices 

CPHIS] 

N/A 

Premapped  modal-actuator  mat.* 

CUCI)  ] 

N/A 

Modal  frequencies*  (asc.  order) 

ICKI) 

N/A 

Modes  assigned  to  Ctlr  1 

IC2(I) 

N/A 

Modes  assigned  to  Ctlr  2 

IC3<I) 

N/A 

Modes  assigned  to  Ctlr  3 

IR<I) 

N/A 

Modes  assigned  to  4th  Group 

Q 

0,1 

O^do  not  print  A,B,C,Q  matrices 

BB<I) 

N/A 

Q  matrix  diagonal  elements 

SKIP 

1,2,3 

1*TR*EA,  2aEA,  3*TR  ** 

DUMMY 

1.2,3 

Conditional  SKIP*** 

CINITCI, J)3 

N/A 

Initial  conditions  state  vector 

DT,TMAX,PDT 

N/A 

TR;  t,  tmax,  output  counter 

CPHIL< J,I)] 

N/A 

*5^03  transposed* 

DUMMY 

1,2,3 

Conditional  SKIP*** 

*  See  the  code 

for  how 

the  data  should  be  formatted. 

**  TR  stands  for  time  response;  EA  for  eigenvalue  analysis. 
***  Covers  all  possible  cases  (see  code  for  implied  meanings) . 


Case  no. : 

Cl) 

(2) 

(3) 

(4) 

1st  pass: 

EA 

TR(TR*EA) 

TR(TR+EA) 

TR 

2nd  pass: 

TR<TR^-EA) 

TR(TR+EA) 

EA 

EA 
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PROGRAM  AC0SS2  < IMPUT , OUTPUT , TAPE6 , TAPE6=0UTPUT , TAPE? ) 


C 

C  pfn:  CSDLLOWERTIME 

C  THIS  PROGRAM  GENERATES  A  LOWER  TRIANGULAR  TRANSFORMATION 
C  THE  THREE  CONTROLLER  SOLUTION  WILL  INCLUDE  RESIDUALS. 

C 

REAL  A1(21,21),A2(21,21),A3(21,21),A4(21,21) 

REAL  B1 (21, 21), B2(21, 21), B3(21,21>, 84(21,21) 

REAL  C1(21,21),C2(21,21),C3(21,21),C4(21,21) 

REAL  CTCCl (21,21), CTCC2 (21,21), CTCC3 (21,21), CTCC4 (21,21) 
REAL  SAT ( 21 , 21 ) , SAT2 (21,21), SAT3 (21,21), SAT4 (21,21) 

REAL  AKC(21,21),ACT(21,21),BCG(21,21),KCC(21,21) 

REAL  P(21,21),S(21,21) 

REAL  QAl(21,21),aA2(21,21),QA3(21.21).QA4(21,21) 

REAL  ACG1(21,21),ACG2(21,21),ACG3(21,21),ACG4(21,21) 

REAL  ABGl (21 , 21 ) , ABG2 (21,21), ABG3(21 , 21 ) , ABG4 (21,21) 

REAL  GAIN1(21,21),GAIN2(21,21),GAIN3(21,21),GAIN4(21,21) 
REAL  KT1(21,21),KT2(21,21),KT3(21,21),KT4(21,21) 

REAL  K0B1(21,21),K0B2(21,21),K0B3(21,21),K0B4(21,21> 

REAL  GANMAl (21,21), GAMNA2 (21,21), GANMA3 (21,21) 

REAL  T2(21,21),T3(21,21),T4(21,21),BB(21) 

REAL  TRT (21 , 21 ) , TEN(21 , 21 ) , CT (21 ,21 ) , V (21 , 21 ) 

REAL  RK(21,21),RK1(21,21),RK2(21,21),RK3(21,21) 

REAL  RG2(21,21),RG3(21,21),RG4(21,21> 

REAL  D(21),W(21),T0L,DT,ABT,TMAX 
REAL  CV(3,21),XL(3) 

REAL  ZETA,CL0S(3,21),SING(21),XTR(21,21),ST0R(21,21) 

REAL  PHIA(21,21),PHIS(21,21),N0DE(2,21),IMIT(4,21) 

REAL  PHIL(3,21),RPHIL(3,21) 

REAL  NA JN  ( 64 , 64 ) ,  STM  ( 64 ,  G4) ,  WORK  (64,64),  XHAT  (3,64) 

REAL  X0(64) ,X1 (64) ,XC(64) ,PETA(21,64) ,CPHIL(3,64) 

INTEGER  N , N2 , NCI , NC2 , NC3 , NC12 , NC22 , NC32 , NR , NR2 , NNODE 
INTEGER  IC1(21),IC2(21),IC3(21),I,J,K,L,M,KK,LL,MM,PDT 
INTEGER  DEC ,  Q ,  NACT ,  NSEN ,  IR(21 ) ,  lER , SKIP ,  NCOL ,  NCOLl 
INTEGER  NDA , NDIM , NDAl , NDINl , ZZ, E2, E3 , E4 , PI , P2 , P3 , AA 
COMPLEX  Z(64),W1(21) 

C0MM0N/MAIN1/NDIM,NDIM1,TEN,X(4096) 

COMMON/MAINA/NDA,NDAl 
CONMON/MAINB/NCOL, NCOLl 
C0MM0N/NAIN2/ST0R 
CONNON/NAIN3/XTR 
C0MM0N/SAVE/T(100) ,TS(100) 

COMMON/ INOUT/KOUT , TAPE 

COMNON/NUN/ ICl , IC2 , ICS , IR , NCI , NC2, NC3 , NR 


INITIALIZATIONS  AND  COMMENTS  ON  DIMENSIONING 

NDIM  IS  THE  ACTUAL  NUMBER  OF  MOOES  (ALSO  REP'D  BY  N>  FOR  THIS  VERSION. 
NDIM  »  20 
NDINl  >  NDIM  «  1 

C  NCOL  IS  THE  MAX  NUMBER  OF  MODES  ENVISIONED  AS  POSSIBLE  TO  RUN 
C  IN  ANY  CONFIGURATION  OF  THE  4  CONTROLLERS.  FOR  THIS  VERSION  ASSIGN 
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C  NOT  LESS  THAN  1  NOR  MORE  THAN  10  MODES  TO  ANY  OF  THE  1ST  3  CTLRS. 

C  THE  A,B,C  MATRICES  NEED  DIM  >>  2  •  «  NODES  ASGN'D  TO  THE  INDIV  CTLR. 

C  SPECIFYING  THIS  PARAN  ALLOWS  DIMENSIONING  OF  THE  INDEXED  VAR'S 
C  TO  A  STD  NUMBER  AND  AIDS  USING  THE  INSL  ROUTINES. 

NCOL  -  21 
NCOLl  >  NCOL  1 

C  NDA  IS  THE  MAX  NUMBER  OF  STATES,  >  4  >  #  NODES  (ETAS)  (-  N  ETAS  « 

C  N  ETA  DOTS  «  N  ERR  «  N  ERR  DOTS)  NDA  IS  BASED  ON  NDIN,  NOT  NCOL, 

C  BECAUSE  IT  REPRESENTS  4  •  THE  ACTUAL  NO.  OF  STATES,  AND  IS 
C  MEANT  TO  SET  UP  SPACE  FOR  'NN',  THE  ACTUAL  NO.  OF  STATES.  MM'S 
C  VALUE  DEPENDS  ON  WHETHER  THERE  ARE  RESIDUALS  (IN  PLACE  OF  A  4TH 
C  CONTROLLER)  OR  NOT  (3  OR  4  ACTIVE  CONTROLLERS).  IF  THE  NO. 

C  OF  RESIDUALS  *  0,  THEN  MM  -  NDA,  ELSE  NN  *  4  ■  «  CTLD  NODES  > 

C  2  •  «  RES  MODES.  NOTE:  REDUNDANTLY,  THIS  VERSION  CONTAINS 
C  THE  EQUIVALENT  TERNS  'MM'  AND  'N'. 

C  WARNING:  UNTIL  SOMEONE  DEVELOPS  A  SUBSTITUTE  FOR  SUBROUTINES  'NEXP, ' 
C  'NNUL,'  AND  'NNULl,'  ETC. 

C  IT  WILL  BE  NECESSARY  TO  MAKE  NDA  AND  ITS  ASSOCIATED  DIMENSIONS  IN  THE 
C  MAIN  DIMENSION  STATEMENTS  EXACTLY  THE  SANE  AS  NN, 

C  AND  RECOMPILE  WHENEVER  THE  VALUE  OF  'MM'  CHANCES  ACCORDING  TO  NEW 
C  CONTROLLER  ASSIGNMENTS.  NAJM  IS  INPUT  TO  NEXP,  A  KLEINNAN  ROUTINE, 

C  WHICH  REQUIRES  THE  MATRIX  TO  HAVE  EXACT  DIMENSIONS,  IN  CONTRAST  TO 
C  THE  IMSL  ROUTINES  WHICH  ALLOW  THE  USER  TO  TELL  THEM  HOW  A  MATRIX  IS 
C  ACTUALLY  DIMENSIONED  RELATIVE  TO  THE  SPACE  IT  WAS  GIVEN  IN  THE 
C  DIMENSION  STATEMENT  IN  THE  MAIN. 

CTEMP«**«  NDA  MUST  EQUAL  NN;  THIS  VERSION  EXPECTS  20  MODES  WITH  8  OF 
CTENP«»*«  THEM  ASSIGNED  TO  CTLR  «  4  WITH  'DEC'  »  3.  THEREFORE,  NN-64. 
NDA  «  64 
NDAl  »  NDA  «  1 
KOUT  =  6 
TAPE  »  9 
Q  -  0 
lER  a  0 

C  ZZ  IS  THE  FLAG  DRIVING  THE  SUPPRESSED  AND  UNSUPPRESSED  PASSES 
ZZ  =  0 

C  AA  IS  THE  FLAG  DRIVING  THE  SINGLE  READ- IN  OF  INITIAL  CONDITIONS 
AA  -  0 


PRINT'(////)' 


PRINT*,' 
PRINT*,' 
PRINT*, ' 

3  0R4C0NTR0LLER 

•  •••• 

PRINT*,' 

»«««« 

(OBSERVER) 

•  •••• 

PRINT*,' 

•  ••«« 

LOWER-RESIDUAL 

PRINT*,' 

BLOCK 

«•••• 

PRINT*,' 
PRINT*,' 
PRINT* , ' 

••••••< 

C  S  D  L  II 

«•••• 

!•••••••• 

PRINT' (//)' 

PRINT*,'  THIS  PROGRAM  GENERATES  A  SOLUTION', 

'  USING  A  LOWER  TRIANGULAR  TRANSFORMATION  ' 
PRINT' (////)' 


on  oooooo  ooo 


C 

C  THIS  CODE  ALLOWS  ASSIGNMENT  OF  THE  NODES  TO  EITHER 
C  ANY  OF  3  CTLRS  PLUS  1  RESID  MODE  SET  OR 

C  ANY  OF  4  CTLRS 
C 
C 

C  INITIAL  SELECTION  FOR  THREE  OR  FOUR  CONTROLLERS 
C 

PRINT*,'  FOR  A  THREE  CONTROLLER  RUN,  ENTER  3,  OR,  ' 
PRINT*,'  FOR  A  FOUR  CONTROLLER  RUN,  ENTER  4  >' 

READ(8,*)  DEC 

DEC  DEFAULT  SWITCH 

IF  (DEC.NE.4)  DEC  =  3 
PRINT*,'  ' 

PRINT*,'  THIS  IS  A  ',DEC,'  CONTROLLER  RUN  ' 


PHI  MATRICES  AND  CONTROLLER  ENTRIES 
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i:  . 


PRINT' <///)' 

IF  (DEC.EQ.S)  THEN 

PRINT*,'  ENTER  NC1,NC2,NC3,NR,NACT,NSEN,ZETA  >' 
ELSE 

PRINT*,'  ENTER  NC1,NC2,MC3,NC4,NACT,NSEN,ZETA  >' 
ENDIF 

READ ( d , * }  NCI , NC2 , NC3 , NR , NACT , NSEN , ZETA 
PRINT* , NCI , NC2 , NC3 , NR , NACT , NSEN , ZETA 
PRINT*,'  ' 


NOW  N  IS  CALCULATED  TO  EQUAL  THE  NUMBER  OF  MODES  IN  THE  MODEL  (’^^NDIM) 
N  =  NCI  ♦  NC2  ♦  NC3  ♦  MR 


PRINT*,'  ENTER  THE  ',NACT,'  ELEMENTS  FOR  EACH  PHIA  ' 
PRINT*,'  ' 

DO  1  1=1, N 

PRINT*, 'ENTER  PHIA  >' 

READ(8,*)  (PHIA(I,J) ,J=1,NACT> 

PRINT*,'  ',(PHIA<I,J),J=1,NACT) 

CONTINUE 
PRINT' (//)' 

PRINT*,'  ENTER  THE  ',NSEN,'  ELEMENTS  FOR  EACH  PHIS  ' 
PRINT*,'  ' 

DO  2  1=1, N 

PRINT*, 'ENTER  PHIS  >' 

READ(8,*)  (PHIS(I,J),J=1,NSEN) 

PRINT*,'  ',(PHIS<I,J),J=1,NSEN) 

CONTINUE 
PRINT' <//)' 
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C  OMEGAS 

C 

C 

PRIMT»,  '  ENTER  THE  VALUE  FOR  EACH  OMEGA  ' 

PRINT*,'  ' 

DO  4  1^1, N 

PRINT*, 'ENTER  OMEGA  >' 

READ(8,*)  W(I) 

PRINT*,'  ',«<!> 

D(I)  =  -2.  *  ZETA  *  W(I) 

4  CONTINUE 
PRINT' (/)' 

C 

C  NODE  ASSIGNMENT  TO  CONTROLLERS 

C 

C 

PRINT' <///>' 

PRINT*,'  THE  FOLLOWING  NODES  ARE  ENTERED  ACCORDING  TO  THE  ' 
PRINT*,'  ORDER  IN  WHICH  THEY  ARE  ENTERED  IN  THE  DATA  FILE  ' 
PRINT*,'  AND  NOT  ACCORDING  TO  THEIR  ACTUAL  MODE  NUMBER.  ' 
PRINT' (//)' 

PRINT*,'  ENTER  THE  ',NC1,'  CONTROLLER  1  NODES  >' 


C 

C 


READCa,*)  <IC1(I),I*1,NC1) 

PRINT*,'  ',(IC1<I),I*1,NC1> 

PRINT*,'  ' 

PRINT*,'  ENTER  THE  ',NC2,'  CONTROLLER  2  MODES  >' 
READ(8,*)  <IC2(I),I»1,NC2) 

PRINT*,'  ',<IC2<I),I=1,NC2> 

PRINT*,'  ' 

PRINT*,'  ENTER  THE  ',NC3,'  CONTROLLER  3  MODES  >' 
READ(8,*)  <IC3<I),I*1,NC3) 

PRINT*,'  ',<IC3(I),I=1,NC3> 

PRINT*,'  ' 

IF  (DEC.E0.3)  THEM 

PRINT*,'  ENTER  THE  ',NR,'  RESIDUAL  MODES  >' 

ELSE 

PRINT*,'  ENTER  THE  ',NR,'  CONTROLLER  4  MODES  >' 
ENDIF 

READ(8,*)  (IR(I),I=1,NR) 

PRINT*,'  ',(IR(I),I=1,NR) 

PRINT*,'  ' 


NC12  -  2  •  NCI 
MC22  =  2  •  MC2 
NC32  =  2  •  NC3 
N2  a  2  *  N 
NR2  >  2  ■  NR 

N  IS  THE  ACTUAL  NUMBER  OF  STATES  (VALUE  DEPENDS  ON  HAVING  RESID  MODES) 
IF  (DEC.EQ.3)  THEN 

M  «  2  *  MC12  ♦  2  *  MC22  ♦  2  *  NC32  ♦  MR2 
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CALL  PRNT(C1,NSEN,NC12) 

PRINT*,'  THE  CONTROLLER  1  WEIGHTING  MATRIX  IS  ' 

CALL  PRNT(QA1,NC12.NC12} 

PRINT*,'  THE  CONTROLLER  2  A  MATRIX  IS  ' 

CALL  PRNT(A2,NC22,NC22) 

PRINT*,'  THE  CONTROLLER  2  B  MATRIX  IS  ' 

CALL  PRNT(B2,NC22,NACT) 

PRINT*,'  THE  CONTROLLER  2  C  MATRIX  IS  ' 

CALL  PRNT(C2,NSEN,NC22) 

PRINT*,'  THE  CONTROLLER  2  WEIGHTING  MATRIX  IS  ' 

CALL  PRNT(QA2,NC22,NC22) 

PRINT*,'  THE  CONTROLLER  3  A  MATRIX  IS  ' 

CALL  PRNT(A3.NC32,NC32) 

PRINT*,'  THE  CONTROLLER  3  B  MATRIX  IS  ' 

CALL  PRNT(B3,NC32,NACT) 

PRINT*,'  THE  CONTROLLER  3  C  MATRIX  IS  ' 

CALL  PRNT(C3,NSEN,NC32) 

PRINT*,'  THE  CONTROLLER  3  WEIGHTING  MATRIX  IS  ' 

CALL  PRNT(aA3,NC32,NC32) 

C 

IF  (NR.EQ.O)  THEN 

PRINT*,'  NO  RESIDUAL  TERMS  ' 

GOTO  115 
ENDIF 
C 

IF  (DEC.EQ.3)  THEN 

PRINT*,'  THE  A  RESIDUAL  MATRIX  IS  ' 

CALL  PRNT(A4,NR2,NR2) 

PRINT*,'  THE  B  RESIDUAL  MATRIX  IS  ' 

CALL  PRNT(B4,NR2,NACT) 

PRINT*,'  THE  C  RESIDUAL  MATRIX  IS  ' 

CALL  PRNT(C4,NSEN,NR2} 

ELSE 

PRINT*,'  THE  CONTROLLER  4  A  MATRIX  IS  ' 

CALL  PRNT(A4,NR2,NR2) 

PRINT*,'  THE  CONTROLLER  4  B  MATRIX  IS  ' 

CALL  PRNT(B4,NR2,NACT) 

PRINT*,'  THE  CONTROLLER  4  C  MATRIX  IS  ' 

CALL  PRNT(C4,NSEN,NR2) 

PRINT*,'  THE  CONTROLLER  4  WEIGHTING  MATRIX  IS  ' 

CALL  PRNT(QA4,NR2,NR2} 

ENDIF 

C 

ENDIF 

C 

115  CONTINUE 
C 
C 

C  THIS  SECTION  GENERATES  THE  RICCATI  SOLUTIONS 
C  AND  THE  GAIN  MATRICES  OF  EACH  CONTROLLER 
C 
C 

CALL  VNULFP ( B1 , B1 , NC12 , NACT , NC12 , NCOL , NCOL , SAT , NCOL , lER ) 


IF  (ZZ.EO.O)  THEN 

CALL  VHULFM (Cl , Cl , NSEN , MC12, MC12 , MCOL , NCOL , CTCCl , NCOL , lER ) 
ENDIF 

120  CONTINUE 
lER  >  0 
TOL  =  0.001 
PRINT' <//>' 

CALL  NRIC(NC12,A1,SAT,QA1,S,ABG1,T0L,IER) 

IF  (ZZ.EO.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  «  BCG  #1  IS  ' 

PRINT*,'  lER  =  ',IER 
CALL  PRNT(S,NC12,NC12) 

ENDIF 

CALL  VNULFN ( B1 , S , NC12 , N ACT , HC12 , NCOL , NCOL , GA IN 1 , NCOL , lER ) 
PRINT*,'  THE  G1  GAIN  MATRIX  IS  ' 

CALL  PRNT(GAIN1,NACT,NC12) 
lER  =  0 
TOL  =  0.001 

CALL  TFR(ACT,A1,NC12,NC12,1,2) 

PRINT' (//)' 

CALL  MRIC ( NC12 , ACT , CTCCl , QAl , P , ACG 1 , TOL , lER ) 

IF  (ZZ.EO.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  -  KCC  #1  IS  ' 

CALL  PRNT(P,NC12,NC12) 

ENDIF 

CALL  NMUL(C1,P,NSEN,NC12,NC12,KT1) 

IF  (ZZ.EO.l)  THEN 

CALL  VMULFP ( RKl , GAHHAl , PI , PI , NSEN , MCOL , NCOL, STOR , NCOL , lER ) 
CALL  NMUL ( STOR , KTl , PI , NSEN , NC12 , KCC ) 

CALL  NHUL ( GAHHAl , KCC , NSEN , PI , NC12 , KTl ) 

PRINT*,'  THE  Kl*  GAIN  MATRIX  IS  ' 

ELSE 

PRINT*,'  THE  Kl  GAIN  MATRIX  IS  ' 

ENDIF 

CALL  TFR(K0B1, KTl, NSEN, NC12, 1,2) 

CALL  PRNT(K0B1,NC12,NSEN) 

125  CONTINUE 

IF  (ZZ.EO.O)  THEN 

CALL  VMULFP ( B2 , B2 , NC22 , N ACT , NC22 , NCOL , NCOL , S AT2 , NCOL , I ER ) 
ENDIF 

IF  (ZZ.EO.O)  THEN 

CALL  VMULFM ( C2 , C2 , NSEN , NC22 , NC22 , NCOL , NCOL , CTCC2 , NCOL , lER ) 
ENDIF 

140  CONTINUE 
lER  =  0 
TOL  =  0.001 
PRINT' (//)' 

CALL  MRIC ( NC22 , A2 , SAT2 , 0A2 , S , ABG2 , TOL , lER ) 

IF  (ZZ.EO.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  #2  IS  ' 

PRINT*,'  lER  =  ',IER 
CALL  PRNT(S.NC22,NC22) 

ENDIF 


CALL  VNULFN ( B2 , S , NC22 , N ACT , NC22 , NCOL , NCOL , G AI N2 . NCOL , lER > 
IF  (ZZ.EQ.l)  THEN 

CALL  VNULFN ( T2 , GAIN2 , NACT , E2 . NC22 , NCOL , NCOL , STOR , NCOL , I ER ) 
CALL  NNUL(RG2, STOR , E2 . E2 . NC22 , TEN ) 

CALL  NNUL ( T2 . TEN , N ACT , E2 , NC22 , GA I N2 ) 

PRINT*,'  THE  G2*  GAIN  NATRIX  IS  ' 

ELSE 

PRINT*,'  THE  G2  GAIN  NATRIX  IS  ' 

ENDIF 

CALL  PRNT(GAIN2,NACT,NC22) 
lER  =  0 
TOL  =  0.001 

CALL  TFR(ACT,A2,NC22,NC22,1,2) 

PRINT' (//)' 

CALL  NRIC ( NC22 , ACT , CTCC2 , QA2 , P , ACG2 , TOL , lER ) 

IF  (ZZ.EQ.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  -  KCC  #2  IS  ' 

CALL  PRNT(P,NC22,NC22) 

ENDIF 

CALL  NNUL  <  C2 . P , NSEN , NC22 , NC22 , KT2 ) 

IF  (ZZ.EQ.l)  THEN 

CALL  VNULFP ( RK2 , GAHHA2 , P2 , P2 , NSEN , NCOL , NCOL , STOR , NCOL , lER ) 
CALL  NNUL ( STOR , KT2 , P2 , NSEN , NC22 , KCC > 

CALL  NNUL ( GANNA2 , KCC , NSEN , P2 , NC22 , KT2  > 

PRINT*,'  THE  K2*  GAIN  NATRIX  IS  ' 

ELSE 

PRINT*,'  THE  K2  GAIN  NATRIX  IS  ' 

ENDIF 

CALL  TFR (K0B2,KT2, NSEN, NC22, 1,2) 

CALL  PRNT(K0B2,NC22,NSEN) 

145  CONTINUE 

IF  (ZZ.EQ.O)  THEN 

CALL  VNULFP ( B3 , B3 , NC32 , NACT , NC32 , NCOL , NCOL , S AT3 , NCOL , I ER ) 
ENDIF 

IF  (ZZ.EQ.0.0R.DEC.EQ.3)  THEN 

CALL  VNULFN ( C3 , C3 , NSEN , NC32 , NC32 , NCOL , NCOL , CTCC3 , NCOL , I ER ) 
ENDIF 

150  CONTINUE 
lER  =  0 
TOL  =  0.001 
PRINT' (//)' 

CALL  NRIC  ( NC32 ,  A3 ,  SAT3 ,  aA3 ,  S ,  ABG3 ,  TOL ,  lER ) 

IF  (ZZ.EQ.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  N3  IS  ' 

PRINT*,'  lER  =  ',IER 
CALL  PRNT(S,NC32,NC32) 

ENDIF 

CALL  VNULFN ( B3 , S , NC32 , N ACT , NC32 , NCOL , NCOL , GA I N3 , NCOL , I ER ) 
IF  (ZZ.EQ.l)  THEN 

CALL  VNULFN  ( T3 ,  G  A IN3 ,  N  ACT ,  E3 ,  NC32 ,  NCOL ,  NCOL ,  STOR ,  NCOL ,  I  ER ) 
CALL  NNUL ( RG3 , STOR , E3 , E3 , NC32 , TEN ) 

CALL  NNUL (T3 , TEN , NACT , E3 , NC32 , GAIN3 ) 

PRINT*,'  THE  G3*  GAIN  NATRIX  IS  ' 


ELSE 

PRINT*,'  THE  G3  GAIN  MATRIX  IS  ' 

ENDIF 

CALL  PRNT( GAINS, NACT,NC32> 
lER  =  0 
TOL  =  0.001 

CALL  TFR(ACT,A3,NC32,NC32,1,2) 

PRINT' <//)' 

CALL  NRIC ( NC32 , ACT , CTCC3 , QA3 , P , ACG3 , TOL , lER ) 

IF  (ZZ.EQ.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  -  KCC  #3  IS  ' 

CALL  PRNT(P,NC32,NC32) 

ENOIF 

CALL  NNUL ( C3 , P , NSEN , NC32 . NC32 , XT3 ) 

IF  (ZZ.EQ.1.AND.DEC.EQ.4)  THEN 

CALL  VMULFP ( RK3 , GAMMAS , P3 , P3 , NSEN , NCOL , NCOL , STOR , NCOL , lER ) 
CALL  MMUL  <  STOR , KT3 , P3 , NSEN , NC32 , KCC  > 

CALL  MMUL ( GAMMA3 . KCC , NSEN , P3 , NC32 , KT3 ) 

PRINT*,'  THE  K3*  GAIN  MATRIX  IS  ' 

ELSE 

PRINT*,'  THE  K3  GAIN  MATRIX  IS  ' 

ENDIF 

CALL  TFR(X0B3,XT3,NSEN,NC32,1,2> 

CALL  PRNT(K0B3,NC32,NSEN> 

155  CONTINUE 

IF  (DEC.EQ.4)  THEN 
IF  (ZZ.EQ.O)  THEN 

CALL  VMULFP  <  B4 , B4 , MR2 , NACT , NR2 , NCOL , NCOL , SAT4 , NCOL , lER  > 
ENDIF 

CALL  VMULFM  <  C4 , C4 , NSEN , NR2 , NR2 , NCOL , NCOL , CTCC4 , NCOL , I ER  > 
160  CONTINUE 
lER  =  0 
TOL  =  0.001 

CALL  NRIC ( NR2 , A4 , SAT4 ,QA4 , S , ABG4 , TOL , lER > 

IF  (ZZ.EQ.O)  THEN 

PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  ♦  BCG  #4  IS  ' 

PRINT*,'  lER  =  ',IER 
CALL  PRNT(S,NR2,NR2) 

ENDIF 

CALL  VMULFM ( B4 , S , NR2 , N ACT , NR2 , NCOL , NCOL , G A IN4 , NCOL , I ER ) 

IF  (ZZ.EQ.l)  THEN 

CALL  VMULFM ( T4 , GA I N4 , N ACT , E4 , NR2 , NCOL , NCOL , STOR , NCOL , I ER ) 
CALL  NNUL (RG4, STOR, E4,E4,NR2, TEN) 

CALL  NNUL  ( T4 ,  TEN ,  N ACT ,  E4 ,  NR2 , GAIN4 ) 

PRINT*,'  THE  G4*  GAIN  MATRIX  IS  ' 

ELSE 

PRINT*,'  THE  G4  GAIN  MATRIX  IS  ' 

ENDIF 

CALL  PRNT(GAIN4,NACT,NR2) 
lER  *  0 
TOL  »  0.001 

CALL  TFR(ACT,A4,NR2,NR2,1,2) 

CALL  NRIC ( NR2 , ACT , CTCC4 , QA4 , P , ACG4 , TOL , lER) 
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CALL  MMUL(C4,P,NSEN,NR2,NR2,KT4) 

CALL  TFR(K0B4,KT4,NSEN,NR2,1,2> 

•  IF  (ZZ.EQ.O)  THEN 

PRINT*,'  THE  K4  GAIN  MATRIX  IS  ' 

ELSE 

PRINT*,'  THE  K4*  GAIN  MATRIX  IS  ' 

ENDIF 

CALL  PRNT(K0B4,NR2,NSEN) 

^  165  CONTINUE 

ENDIF 

C 

C 

C  THIS  SECTION  GENERATES  THE  BLOCK  SEGMENTS 
C  OF  NAJM  AND  PUTS  THEM  INTO  THE  MAJH  MATRIX 

•  c 

C  THE  THREE  CONTROLLER  MATRIX  WILL  CONTAIN 
C  RESIDUAL  TERMS  (SEE  DIAGRAM  BELOW). 

C  THE  FOUR  CONTROLLER  MATRIX  DOES  NOT  IN- 
C  CLUDE  RESIDUALS  (YET). 

C 

^  C  THE  THREE  CONTROLLER  MATRIX  (MAJM)  WITH 

C  RESIDUAL  TERMS  WILL  LOOK  LIKE: 

C 


c 

A1*BG1 

BlGl 

B1G2 

B1G2 

B1G3 

B1G3 

0 

• 

c 

c 

0 

Al-KCl 

K1C2 

0 

K1C3 

0 

KICR 

c 

c 

B2G1 

B2G1 

A2*BG2 

B2G2 

B2G3 

B2G3 

0 

c 

c 

K2C1 

0 

0 

A2-KC2 

K2C3 

0 

K2CR 

A 

c 

c 

B3G1 

B3G1 

B3G2 

B3G2 

A3*BG3 

B3G3 

0 

c 

c 

K3C1 

0 

K3C2 

0 

0 

A3-KC3 

K3CR 

c 

c 

BRGl 

BRGl 

BRG2 

BRG2 

BRG3 

BRG3 

AR 

K  »  2  •  NC12 

KK  =  K  ♦  MC22 

L  *  2  *  NC22  ♦  K 

LL  =  L  ♦  MC32 

PI  »  2  *  MC32  ♦  L 

IF  (DEC.EQ.3)  THEM 

HN  «  2*NC12  ♦  2*MC22  ♦  2*NC32  ♦  NR2 

ELSE 

MM  »  2*MC12  ♦  2*NC22  ♦  2*NC32  ♦  2*NR2 

P2  »  PI  ♦  NR2 

ENDIF 
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DO  200 

DO  200  J=1,MM 
NAJH(I,J)  =  0.0 
DO  201  I=1,NC12 
DO  201  J»1,NC12 
NAJM(I.J)  s  ABG1<Z,J) 

DO  202  I»1,NC22 
DO  202  J>1,NC22 
MAJM(1^K,J«K>  -  ABG2(I,J) 

DO  203  I»1,NC32 
DO  203  J=1,MC32 
MAJM<I-»L.J«L)  =  ABG3<I,J) 

CALL  TFR(AKC,ACG1,NC12,NC12.1,2) 

DO  204  I>1,NC12 
DO  204  J>1,NC12 
NAJM<I*MC12,J*(IC12)  =  AKC(I,J) 

CALL  TFR(AKC,ACG2,NC22,NC22,1,2) 

DO  205  I=1,NC22 
DO  205  Jsl,NC22 
MAJMCI  +  KK.J^^KK)  =  AKC(I,J> 

CALL  TFR<AKC,ACG3,NC32,NC32,1.2) 

DO  206  I>1,NC32 
DO  206  J=1,NC32 
MAJM(I*LL,J>LL)  =  AKC(I,J) 

CALL  MMUL<B1,GAIN1,NC12,NACT,NC12,BCG) 
DO  207  I»1,MC12 
00  207  J»1,MC12 
HAJN(I.J4’NC12)  s  BCG(I,J) 

CALL  NMUL ( B1 , GAIN2 , NC12 , N ACT , NC22 , BCG ) 
DO  208  I-1,NC12 
DO  208  J=1,NC22 
I1AJM<I,J*K)  =  BCG(I,J) 
t!AJII<I,J*KK)  =  BCG<I,J) 

CALL  HMUL ( B1 , GAIN3 , NC12 , N ACT , NC32 , BCG ) 
DO  209  I=1,NC12 
DO  209  J=1,NC32 
HAJH<I,J*L)  =  BCG(I,J) 

«AJH(I,J*LL)  =  BCG<I,J) 

CALL  HHUL ( B2 , GAINl , NC22 , N ACT , NC12 , BCG  > 
DO  210  I=1,NC22 
DO  210  J=1,NC12 
NAJNd^K.J)  =  BCG(I,J) 

HAJM(I*K,J*NC12)  =  BCG<I,J) 

CALL  HHUL ( B2 , GAIN2 , NC22 , NACT , NC22 , BCG > 
DO  211  I»1,NC22 
DO  211  J»1,NC22 
NAJN(I«K,J4^KK)  »  BCG(I,J) 

CALL  HHUL ( B2 , GA I N3 , NC22 , N ACT , NC32 , BCG  > 
DO  212  I*1,MC22 
DO  212  J»1,NC32 
NAJN(I*K,J«L)  >  BCG(1,J) 

MAJH(I*K,J*LL)  »  BCGd.J) 
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CALL  MMUL ( B3 , GAIN 1 , NC32 , N ACT , NC12 , BCG> 
DO  213  I=1,MC32 
DO  213  J=1,NC12 
MAJM<I*L,J)  =  BCG<I,J> 

«AJM<I*L,J+NC12)  =  BCG<I,J) 

CALL  NHUL ( B3 , GAIN2 , NC32 , N ACT , NC22 , BCG ) 
DO  214  I-1,NC32 
DO  214  J-1,NC22 
NAJH(l4^L,J4^K)  -  BCG(I,J) 

MAJM(I*L,J*KK)  =  BCG<I,J) 

CALL  NNUL ( B3, GAIN3, NC32 , NACT, NC32, BCG) 
DO  215  I»1,NC32 
DO  215  J»1,MC32 
MAJM(I*L,J*LL)  =  BCG<I,J) 

CALL  MMUL ( KOBl , C2 . NC12 . NSEN , NC22 , KCC ) 
DO  216  I=1,NC12 
DO  216  J-1,NC22 
HAJM(I*NC12,J*K)  =  KCC<I,J) 

CALL  NNUL ( KOBl . C3 , NCI 2 , NSEN » NC32 , KCC ) 
DO  217  1=  1,NC12 
DO  217  J=  1,NC32 
NAJK<I^NC12,J*L)  =  KCCd.J) 

CALL  NNUL ( K0B2 , Cl , NC22 . NSEN , NC12, KCC) 
DO  218  I-1,NC22 
DO  218  J»1,NC12 
NAJN<I*KK,J)  s  KCC(I,J) 

CALL  NNUL ( KOB2 . C3 , NC22 , NSEN , NC32 , KCC) 
DO  219  I»1,NC22 
DO  219  J»1,NC32 
NAJN(I*KK,J>L)  =  KCCd.J) 

CALL  NNUL ( K0B3 . Cl , NC32 . NSEN . NCI 2 . KCC ) 
DO  220  I=1,NC32 
DO  220  J=1,NC12 
NAJKd*LL,J)  =  KCCd.J) 

CALL  NNUL ( K0B3 . C2 . NC32 , NSEN . NC22 . KCC ) 
DO  221  I=1.NC32 
DO  221  J=1.NC22 
NAJNd+LL.J+K)  =  KCCd.J) 

CALL  NNUL ( B4 . GAINl . NR2 . NACT . NC12 . BCG ) 
DO  222  1=1. NR2 
DO  222  J=1.NC12 
HAJNd*Pl.J)  =  BCGd.J) 

HAJNd  +  Pl.J*NC12)  =  BCGd.J) 

CALL  NNUL ( B4 . GA IN2 . NR2 . NACT . NC22 , BCG ) 
DO  223  1=1. NR2 
DO  223  J=1.NC22 
NAJNd*Pl,J*K)  =  BCGd.J) 
NAJNd^Pl.J^KK)  =  BCGd.J) 

CALL  NNUL  ( B4 .  G A IN3 .  NR2 .  N ACT .  NC32 .  BCG ) 
DO  224  1=1. NR2 
DO  224  J=1.NC32 
NAJNd+Pl.J+L)  =  BCGd.J) 
NAJNd*Pl.J*LL)  »  BCGd.J) 
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CALL  NMUL  ( KOBl . C4 ,  NC12 ,  NSEN ,  NR2 .  KCC ) 
DO  225  I=1,NC12 
DO  225  J=1,NR2 

225  MAJII(I*MC12,J*P1)  =  KCC(I,J) 

CALL  NHUL ( KOB2 , C4 . NC22 , NSEN , NR2 , KCC > 
DO  226  I-1,NC22 
DO  226  J-1,NR2 

226  MAJH(I*KK,J*P1)  =  KCC(I,J) 

CALL  NMUL ( K0B3 . C4 . NC32 , NSEN . NR2 , KCC  > 
DO  227  I=1,MC32 
DO  227  J=1,NR2 

227  MAJM(I*LL,J*P1)  =  KCCd.J) 


IF  (DEC.EQ.4}  THEN 


DO  230  1=1, NR2 
DO  230  J=1,NR2 

230  MAJM(I*P1,J*P1)  =  ABG4<I,J) 

CALL  TFR(AKC,ACG4,NR2,NR2,1,2) 

DO  231  1=1, NR2 

DO  231  J=1,NR2 

231  MAJM(I*P2,J*P2)  =  AKC(I,J) 

CALL  NMUL ( B1 , GAIN4, NC12 , NACT , NR2 , BCG) 
DO  232  I»1,NC12 
DO  232  J=1,NR2 
MAJM<I,J*P1)  =  BCG<I,J) 

232  MAJM<I,J*P2)  =  BCG(I,J) 

CALL  NHUL ( B2 , GAIN4 , NC22 , NACT , NR2 , BCG > 
DO  233  I=1,NC22 
DO  233  J=1,MR2 
NAJM(I*K,J+P1)  =  BCGCI.J) 

233  NAJN(I*K,J*P2)  =  BCG(I,J) 

CALL  NNUL ( B3 , GAIN4 , NC32 , N ACT , NR2 , BCG ) 
DO  234  I=1,MC32 
DO  234  J=1,NR2 
MAJM<I»L,J*P1)  =  BCG(I,J) 

234  MAJM(I*L,J+P2)  =  BCG(I,J) 

CALL  NNUL ( B4 , GA IN4 , NR2 , NACT , NR2 , BCG  > 
DO  238  1  =  1, NR2 
DO  238  J=1,NR2 

238  MAJM(I*P1,J*P2)  =  BCG(I,J) 

CALL  NNUL ( K0B4 , Cl , NR2 , NSEN , NC12 , KCC ) 
DO  242  1=1, NR2 
DO  242  J=1,NC12 

242  NAJM(I*P2,J)  =  KCC(I,J) 

CALL  NHUL ( K0B4 , C2 , NR2 , NSEN , NC22 , KCC ) 
DO  243  1=1, NR2 
DO  243  J=1,NC22 

243  HAJM(I*P2,J*K>  =  KCC(I,J) 

CALL  NMUL  <  K0B4 , C3 , NR2 , NSEN , NC32 , KCC ) 
DO  244  1=1, NR2 


o  o  o  o  o 


DO  244  J*1,MC32 
244  MAJM(I+P2,J*L)  =  KCC(I,J) 

C 

ELSE 

C 

CALL  F0RMA(A4,D,U,NR.NR2,IR} 
DO  250  I»1,NR2 
DO  250  Jsl,NR2 

250  MAJM(I^P1,J*P1)  =  A4(I,J) 
ENDIF 


CONSTRUCTIOM  OF  THE  MAJOR  MATRIX  'MAJM'  IS  MOW  COMPLETE 


IF  (DEC.EQ.4)  THEN 

PRINT*,'  THE  FOUR  CONTROLLER  MAJM  IS  ' 

ELSE 

PRINT*,'  THE  THREE  CONTROLLER  MAJM  W/RESIDUALS  IS  ' 

ENDIF 

CALL  PRNTXL(MAJN,NM,NM} 

C 

C 

C  NEXT,  THE  TIME  RESPONSE  SECTION  FOLLOWED  BY  THE  E' VALUE  ANALYSIS  SECTION 
C  THE  USER  HAS  THE  OPTION  OF  GETTING  BOTH  TIME  RESP  AND  E'VALUE  ANALYSIS 
C  OR  'SKIP' PING  EITHER  OF  THEM.  IF  TIME  RESP  IS  SKIPPED,  HO  IC'S  NEEDED. 

C 

C 

IF  (ZZ.EQ.l)  THEN 

PRINT*,'  THE  SUPPRESSED  ANALYSIS  IS  CALCULATED  ' 

PRINT'(//)' 

ENDIF 

PRINT*,'  FOR  THE  TIME  RESPONSE  AND  THE  EIGENVALUE' 

*  ,'  ANALYSIS  ENTER  1  ' 

PRINT*,'  FOR  ONLY  THE  EIGENVALUE  ANALYSIS  ENTER  2  >  ' 

PRINT*,'  FOR  ONLY  THE  TIME  RESPONSE  ENTER  3  >  ' 

READ (8,*)  SKIP 
IF  (SKIP.EQ.l)  THEN 

PRINT*,'  YOU  ASKED  FOR  BOTH  TIME  RESP  AND  EIGENVALUE  ANAL' 

PRINT*,'  ' 

PRINT*,'  THE  TIME  RESPONSE  IS  FIRST...  ' 

ENDIF 

IF  (SKIP.EQ.2)  THEN 
PRINT*,'  ' 

PRINT*,'  YOU  CHOSE  "2"  SO  WE  WILL  FORGET  THE  TIME  RESPONSE  ' 

GOTO  300 
ENDIF 

IF(SKIP.EQ.3)  THEN 

PRINT*,'  YOU  ASKED  FOR  THE  TIME  RESPONSE  ONLY' 

PRINT*,'  ' 

ENDIF 

C  INITIAL  CONDITIONS 

C  THESE  ARE  READ  IN  ONLY  ONCE  FOR  EACH  JOB  (WHEN  AA*0.) 


C  REGARDLESS  OF  THE  NUMBER  OF  CONSECUTIVE  RUNS 
C 

C  REMEMBER  --  EVEN  THOUGH  RESIDUAL  MODES  DO  NOT  HAVE  E  AND 
C  E  DOT  TERMS,  THIS  READ  STATEMENT  WILL  BE  LOOKING  FOR  THEM, 

C  SO  INPUT  E  AND  E  DOT  TERMS  FOR  THE  RESIDUALS  ALSO. 

C  DON'T  WORRY  —  SUBROUTINE  FORMXO  WILL  FILTER  THEM  OUT  LATER. 

C 

C 

IF(AA.EQ.O)  THEN 

C  READ  IN  A  DUMMY  DATA  CARD.  ALLOWS  STRUCTURING  OF  INP  DATA  FILE  THE 
C  SAME  WAY  FOR  ALL  4  CASES  OF  DESIRED  OUTPUT:  TIME  RESP/EIG  ANAL  VS 
C  ONE  OR  THE  OTHER  ACCORDING  TO  SUPPRESSED  OR  UNSUPPRESSED  PASS. 

C 


c 

n 

CASE  1 

CASE  2 

CASE  3 

CASE 

4 

c 

EA 

TR<TR*EA) 

TR(TR*EA> 

EA 

ZZ  =  0 

c 

c 

TR(TR+EA) 

TR<TR*EA) 

EA 

EA 

ZZ  =  1 

c 

c 

IMPLIED  MEANING 

1  OF  THE  DATA  CARDS  IN  THAT 

AREA  OF  THE 

INPUT  DATA 

c 

PASS  WITH 

I  ZZ  =  0 

O 

c 

PASS  WITH 

[  ZZ  =  1 

c 

c 

SKIP 

SKIP 

SKIP 

SKIP 

c 

SKIP 

DUMMY 

DUMMY 

SKIP 

c 

c 

DATA  BLKS 

DATA  BLKS 

DATA  BLKS 

DATA 

BLKS(EXTRAN) 

c 

DUMMY (EXTR AN) 

SKIP 

SKIP 

DUMMY 

(EXTRANEOUS) 

IF(ZZ.EQ.O) 

READ(8,*) 

DUMMY 

PRINT*,'  ENTER  THE  INITIAL  CONDITIONS  FOR  ',N,'  MODES  ' 

PRINT*,'  ROWl=ETA,  R0W2=ETA  DOT,  R0W3=E,  ROW4=E  DOT  ' 

DO  90  1=1,4 

90  READ<S,*)  (INIT(I,J),J=1,N) 

PRINT' (//)' 

PRINT*,'  THE  INITIAL  CONDITIONS  ARE  ' 

DO  95  1=1,4 

95  PRINT*, <INIT<I,J),J=1,N) 

PRINT' (//)' 

PRINT*,'  THE  INITIAL  STATE  VECTOR.  Z  (XO)  IS  ' 

CALL  FORHXO<XO,INIT,DEC) 

CALL  PRNT(X0,M,1) 

C  READ  IN  THE  TIME  PARAMETERS 
READ<8,*)  DT,THAX,PDT 

C  READ  IN  THE  PHI -LOS  MATRIX,  SIMULTANEOUSLY  FORMING  IT  TO  (3XN) 

DO  9B  1=1, N 

98  READ(8,*)  (PHIL(J,I) ,J=1,3) 

C  REFORM  THE  PHI -LOS  MATRIX  I AW  THE  REORDERING  OF  NODES 
C  SIMULTANEOUSLY  FORMING  THE  PRODUCT  CPHIL  =  RPHIL*PETA  ONCE  AND  FOR  ALL 
CALL  RFMPHIL  <  PHIL , RPHIL , NCI , NC2 , NC3 , NR , ICl , IC2 , IC3 , IR , 


•PETA,CPHIL,MM> 

AA  =  1 
ENDIF 

PRINT*,'  THE  TIME  INCREMENT  IS  ',DT,'  SECONDS  ' 

TOL  =  0.001 

CALL  NEXP(NN.NAJN,DT.STN) 

PRINT' (//)' 

PRINT*,'  THE  SOLUTION  IS  THE  STATE  TRANS  MATRIX,  STM  ' 
CALL  PRNTXL(STN,MM,MM) 

C 

C 

C  STM  IS  NOW  THE  SOLUTION  TO  ZDOT  »  MAJM  *  Z 
C  WE  NOW  PROPAGATE  THE  STATE  IN  OT  STEPS 
C 

C  FIRST  GENERATE  THE  [C1:C2:C31  MATRIX 

C 

C 


CS 

DO  260  Isl,NSEN 

cs 

DO  260  J=1,NC1 

CS 

260 

V(I,J)=C1(I,J) 

cs 

DO  270  I»1,NSEN 

cs 

DO  270  J=1,NC2 

cs 

270 

V(I,J*NC1)=C2(I,J> 

cs 

DO  280  I*1,NSEN 

cs 

DO  280  J-1,NC3 

cs 

280 

V(I,J*NC1*NC2)*C3<I,J) 

cs 

NN0DE>NC1*NC2*NC3 

cs 

IF  (NR.GT.O)  THEN 

cs 

DO  290  I«1,NSEN 

cs 

DO  290  J»1,NR 

cs 

290 

V(I,J*NM0DE)=C4(I,J) 

cs 

NNODE-NNODE^NR 

cs 

n 

ENDIF 

cs 

PRINT*,'  THE  C  PARTITION  MATRIX  IS  ' 

cs 

n 

CALL  PRNT(V,NSEN,NNODE) 

c 

cs 

CALL  YHAT ( CLOS , NSEN , NCI , NC2 , NC3 , NR , NN , XH AT , PETA , CV , V ) 

C 

CALL  TIMEX(STM,MM,DT,X0,PDT,TMAX,X1,XC,CPHIL,XL) 
IF<SKIP.EQ.3)  GOTO  410 


C 

C 

C  END  OF  TINE  RESPONSE  SECTION 

C 

C 

300  CONTINUE 
C 
C 

C  EIGENVALUE  ANALYSIS  SECTION 

C 

C 
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PRIMT' <///)' 

PRINT*,'  OVERALL  SYSTEM  EIGENVALUES  ' 

NOTICE  EIGRF  INPUT  HERE  IS  MAJM  WHICH  HAS  ACTUAL  DIMENSION  'MM' 
SEE  COMMENTS  ON  DIMENSIONING  NEAR  BEGINNING  OF  MAIN 
CALL  EIGRF ( MAJM . MM . NDA , 0 , 2 , TEN . NCOL , WORK , lER  > 

PRINT*,'  lER  «  ',IER 
DO  400  1*1, MM 

400  PRINT*,'  ',2(1) 

PRINT' <//)' 

C 

PRINT*,'  EIGENVALUES  OF  AC  ♦  BCG  SYSTEM  1  ' 

CALL  EIGRF ( ABGl , NC12, NCOL . 0 , W1 , TEN. NCOL , STOR , lER) 

PRINT*,'  lER  =  ',IER 
DO  401  I=1,HC12 

401  PRINT*,'  ',W1(I) 

PRINT' <//)' 

C 

PRINT*,'  EIGENVALUES  OF  AC  -  KCC  SYSTEM  1  ' 

CALL  TFR(AKC,ACG1,NC12,NC12,1,2) 

CALL  EIGRF ( AKC , MC12 , NCOL , 0, W1 , TEN , NCOL , STOR , lER ) 

PRINT*,'  lER  =  ',IER 
DO  402  I*1,NC12 

402  PRINT*,'  ',W1(I) 

PRINT' (//)' 

C 

PRINT*,'  EIGENVALUES  OF  AC  ♦  BCG  SYSTEM  2  ' 

CALL  EIGRF  <  ABG2 , MC22 , NCOL , 0 , Wl , TEN , NCOL , STOR , lER  > 

PRINT*,'  lER  »  ',IER 
DO  403  I=1,NC22 

403  PRINT*,'  ',W1(I) 

PRINT' (//)' 

C 

print*,'  eigenvalues  of  AC  -  KCC  SYSTEM  2  ' 

CALL  TFR<AKC,ACG2,NC22,NC22,1,2> 

CALL  EIGRF ( AKC , MC22 , NCOL , 0 , Wl , TEN , NCOL , STOR , lER ) 

PRINT*,'  lER  =  ',IER 
DO  404  I=1,NC22 

404  PRINT*,'  ',W1<I) 

PRINT' (//)' 

C 

print*,'  EIGENVALUES  OF  AC  »  BCG  SYSTEM  3  ' 

CALL  El GRF  <  ABG3 , NC32 , NCOL , 0 , Wl , TEN , NCOL , STOR , I ER ) 

PRINT*,'  lER  =  ',IER 
DO  405  I=1,NC32 

405  PRINT*,'  ',W1<I) 

PRINT' (//)' 

C 

print*,'  eigenvalues  of  AC  -  KCC  SYSTEM  3  ' 

CALL  TFR(AKC,ACG3,NC32,NC32,1,2) 

CALL  EIGRF ( AKC , NC32 , NCOL , 0 , Wl , TEN , NCOL , STOR , lER ) 

PRINT*,'  lER  *  ',IER 
DO  406  I*1,NC32 

406  PRINT* , ' 


WKI) 


id 


PRINT' (//)' 

IF  (NR.EQ.O)  THEN 

PRINT*,'  NO  RESIDUAL  TERM  EIGENVALUES  ' 

GOTO  410 
ENDIF 

IF  (DEC.EQ.4)  THEN 

PRINT*,'  EIGENVALUES  OF  AC  ♦  BCG  SYSTEM  4  ' 

CALL  EIGRF  <  ABG4 , NR2 , NCOL , 0 , U1 , TEN , NCOL  r STOR , lER) 
PRINT*,'  lER  »  ',IER 
DO  407  1-1, NR2 

407  PRINT*,'  •  ',W1<I) 

PRINT' (//)' 

PRINT*,'  EIGENVALUES  OF  AC  -  KCC  SYSTEM  4  ' 

CALL  TFR(AKC,ACG4,NR2,NR2.1,2) 

CALL  EIGRF ( AKC , NR2 , NCOL , 0 , W1 , TEN , NCOL , STOR , I ER ) 
PRINT*,'  lER  *  '.lER 
DO  406  1*1, NR2 

408  PRINT*,'  ',W1<I) 


PRINT*,'  EIGENVALUES  OF  THE  A  RESIDUAL  MATRIX 

CALL  EIGRF ( A4 , NR2 , NCOL , 0 , W1 , TEN , NCOL , STOR , lER > 

PRINT*,'  lER  *  ',IER 
DO  409  1*1, NR2 

409  PRINT*,'  ',W1<I) 

ENDIF 

PRINT' (///)' 

C  END  OF  EIGENVALUE  ANALSIS  SECTION 

410  CONTINUE 

IF  (ZZ.EO.l)  GOTO  600 


C  THIS  SECTION  FORMS  THE  TRANSFORMATION  MATRICES. 
C  TO  GET  MAJM  IN  LOWER  TRIANGULAR  FORM,  IT  IS 
C  NECESSARY  TO  DRIVE  THE  B1G2,  B1G3,  B2G3,  K1C2, 

C  K1C3,  AND  K2C3  TERMS  TO  ZERO  (THREE  CTLRS) . 

C 

C  WHEN  FOUR  CONTROLLERS  ARE  USED,  THERE  WILL 
C  BE  ADDITIONAL  G4  AND  K4  TERMS.  IT  WILL 
C  THEN  BE  NECESSARY  TO  DRIVE  THE  B(I)G4,  AND 
C  K(I)C4  TERMS  TO  ZERO  ALSO. 


C  AFTER  THE  TRANSFORMATION  IS  COMPLETE, 

C  THE  THREE  CONTROLLER  MAJM  (WITH  RESIDUALS) 
C  WILL  LOOK  LIKE: 

C 

C  *••• 

c  • 
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(, 


ri 

)5 


Al^BGl 

BlGl 

0 

0 

0 

0 

0 

Al-KCl 

0 

0 

0 

0 

B2G1 

B2G1 

A2«BG2 

B2G2 

0 

0 

K2C1 

0 

0 

A2-KC2 

0 

0 

B3G1 

B3G1 

B3G2 

b3G2 

A34^BG3 

B3G3 

K3C1 

0 

K3C2 

0 

0 

A3-KC3 

BRGl 

BRGl 

BRG2 

BRG2 

BRG3 

BRG3 

:re  the 

NON-ZERO 

TERMS 

INCLUDE  THE 

C  TRANSFORMATION  MATRICES. 

C  (THE  RESIDUALS  ARE  NOT  ZEROED  SINCE  THEY 
C  ARE  NOT  ADDRESSED  IN  THIS  PROGRAM) 


C  ON  WITH  THE  TRANSFORMATION  MATRICES! 

C 

C  FIRST  THE  OBSERVER  GAIN  MATRIX,  K1 

C 

c 

CALL  TFR(CT,C2,NSEN,NC22,1,2) 

DO  500  I«1,NC2 
DO  500  J»1,NSEN 

500  V(I,J)  =  CTd.J) 

CALL  TFR<CT,C3,NSEN,NC32,1,2) 

DO  501  1=1, NC3 
DO  501  J=1,NSEN 

501  V(I*NC2,J)  =  CT(I,J) 

IF  (DEC.EQ.4}  THEN 

CALL  TFR(CT,C4,NSEN,NR2,1,2) 

DO  502  1=1, NR 
DO  502  J=1,NSEN 

502  V(I+NC2*NC3,J)  =  CTd.J) 

MRV  =  NC2  ♦  NC3  ♦  NR 
PRINT*,'  V  (C2/C3/C4)  IS  ' 

ELSE 

PRINT*,'  V  (C2/C3)  IS  ' 

NRV  =  NC2  ♦  NC3 
ENDIF 

CALL  PRNT(V,NRV,NSEN) 

CALL  LSVDF (V.NCOL, NRV, NSEN, TEN, NCOL,-!, SING, STOR.IER) 
PRINT*,'  ' 

PRINT*,'  LSVDF  1  lER  »  ',IER 
PRINT' (//)' 
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PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT(V,NSEN,NSEN) 

PI  =  NSEN  -  NRV 
IF  (Pl.LT.l)  THEN 
DO  503  I -1,  NSEN 

503  GANMAKI.l)  =  V(I,NSEN) 

PI  »  1 

ELSE 

DO  504  I-1,NSEN 
DO  504  J»1,P1 

504  GANMAKI.J)  -  V(I,J*NRV) 

ENDIF 

PRINT*,'  TRANSFORHATION  MATRIX  GAMMAl  ' 

CALL  PRNT(GAMMA1,NSEN,P1) 

C 

C  CHECK  TO  SEE  THAT  GAMMAl  IS  ORTHOGONAL  TO  BOTH  C2  AND  C3 
C 

C  NOTE:  AKC  IN  THIS  SECTION  IS  JUST  A  WORK  AREA  TO  TEST 

C  THE  ORTHOGONALITY  OF  CT  *  TR.  IN  ALL  CASES  IT 

C  SHOULD  BE  A  BLOCK  ZERO  MATRIX. 

C 

CALL  TFR(CT,C2, NSEN, NC22, 1,2) 

CALL  MMUL ( CT , GAMMAl , NC22 , NSEN , PI , AKC > 

PRINT*,'  C2T  *  GAMMAl  ' 

CALL  PRNT(AKC,NC22,P1) 

CALL  TFR<CT, C3, NSEN, NC32, 1,2) 

CALL  MMUL (CT , GAMMAl , NC32 , NSEN , PI . AKC) 

PRINT*,'  C3T  *  GAMMAl  ' 

CALL  PRNT<AKC,NC32,P1) 

IF  (DEC. £0.4)  THEN 

CALL  TFR (CT,C4, NSEN, NR2, 1,2) 

CALL  MMUL  ( CT ,  GAMMAl ,  NR2 ,  NSEN ,  P 1 ,  AKC ) 

PRINT*,'  C4T  *  GAMMAl  ' 

CALL  PRNT(AKC,NR2,P1) 

ENDIF 

C 

IF  (DEC.EQ.4)  THEN 

PRINT*,'  C234  SINGULAR  VALUES  ' 

ELSE 

PRINT*,'  C23  SINGULAR  VALUES  ' 

ENDIF 

CALL  PRNT(SING,NRV,1) 

CALL  TFR(TRT, GAMMAl, NSEN, PI, 1,2) 

CALL  MMUL  ( TRT , GAMMAl , PI ,  NSEN , PI ,  RK ) 

CALL  GMINV(P1,P1,RK,RK1,J,TAPE) 

CALL  TFR(CT,C1,NSEN,MC12,1,2) 

CALL  MMUL  ( THT , Cl , PI , NSEN , NC12 , AKC) 

CALL  MMUL ( CT , GAMMAl , NCI 2 , NSEN , PI , KOB 1 ) 

CALL  MNUL(K0B1,RK1,NC12,P1,P1,ST0R) 

CALL  MMUL  ( STOR ,  AKC ,  NCI  2 ,  PI ,  NC12 ,  CTCCl ) 

C 

C  THIS  CTCCl  WILL  BE  SUBSTITUTED  BACK  INTO  MRIC 
C  SYSTEM  1  TO  GET  A  NEW  Kl. 


I 


i 

*r 


r 

i 


§1 

I 

t 

*al 


I 


I 


NOW  THE  OBSERVER  GAIN  MATRIX,  K2 


CALL  TFR(CT,C3,NSEN,NC32,1.2) 

DO  505  Isl,NC3 
DO  505  J»1,NSEN 

505  V(I,J)  =  CT(I,J) 

IF  (DEC.EQ.4)  THEN 

CALL  TFR(CT,C4,NSEN,NR2,1,2} 

DO  506  I»1,NR 
DO  506  J»1.NSEN 

506  V(I*NC3,J)  =  CTd.J) 

NRV  =  NC3  ♦  NR 
PRINT*,'  V  (C3/C4)  IS  ' 

ELSE 

NRV  >  NC3 

PRINT*,'  V  <C3)  IS  ' 

ENDIF 

CALL  PRNT(V,NRV,NSEN> 

CALL  LSVDF  < V , NCOL , NRV , NSEN , TEN , NCOL , - 1 , SI NG , STOR , lER ) 
PRINT*,'  ' 

PRINT*,'  LSVDF  2  lER  «  ',IER 
PRINT' <//)' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT(V, NSEN. NSEN) 

P2  »  NSEN  -  NRV 
IF  {P2.LT.1)  THEN 
DO  507  1*1, NSEN 

507  GAMMA2<I,1)  =  V(I,NSEN) 

P2  =  1 

ELSE 

DO  508  1*1, NSEN 
DO  508  J=1,P2 

508  GAHHA2(I,J>  *  V(I,J*NRV) 

ENDIF 

PRINT*,'  TRANSFORMATION  MATRIX  GAMMA2  ' 

CALL  PRNT(GAMMA2,NSEN,P2) 

C 

C  CHECK  TO  SEE  THAT  GAHMA2  IS  ORTHOGONAL  TO  C3 
C 

CALL  TFR(CT,C3, NSEN, NC32, 1,2) 

CALL  MMUL ( CT , GAMMA2 , NC32 , NSEN , P2 , AKC  > 

PRINT*,'  C3T  *  GAHHA2  ' 

CALL  PRNT(AKC,NC32,P2) 

IF  (DEC.EQ.4)  THEN 

CALL  TFR(CT, C4, NSEN, NR2, 1,2) 

CALL  MMUL  ( CT ,  G AHMA2 ,  NR2 ,  NSEN ,  P2 ,  AKC  > 

PRINT*,'  C4T  *  GAMMA2  ' 

CALL  PRNT(AKC,NR2,P2) 

ENDIF 
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' 


o  o  o 


IF  (DEC.EQ.4)  THEN 

PRINT*,'  C34  SINGULAR  VALUES  ' 

ELSE 

PRINT*,'  C3  SINGULAR  VALUES  ' 

ENDIF 

CALL  PRNT(SING.NRV,1) 

C 

CALL  TFR(TRT,GANNA2,NSEN,P2.1,2) 

CALL  NNUL ( TRT , G ANMA2 , P2 , NSEN , P2 . RK  > 

CALL  GMINV(P2,P2.RK.RK2,J,TAPE) 

CALL  TFR(CT,C2,NSEN,NC22,1,2> 

CALL  NHUL  <  TRT , C2 , P2 , NSEN . NC22 , AKC  > 

CALL  NNUL ( CT , GANNA2 , NC22 , NSEN , P2 , K0B2} 

CALL  NNUL ( K0B2 , RK2 , NC22 , P2 . P2 , STOR ) 

CALL  NNUL ( STOR . AKC , NC22 , P2, NC22 , CTCC2  > 

C 

C 

C  CTCC2  WILL  BE  SUBSTITUTED  BACK  INTO 
C  NRIC-SYSTEN  2  FOR  A  NEW  K2. 

C 

C 

C  NOW  THE  OBSERVER  GAIN  NATRIX,  K3.  WHEN 
C  USING  FOUR  CONTROLLERS 
C 
C 

IF  (DEC.EQ.4)  THEN 

CALL  TFR <CT,C4, NSEN, NR2, 1,2) 

DO  509  1*1, NR 
DO  509  J«1,NSEN 

509  V(I,J)  »  CT<I,J) 

PRINT*,'  V  <C4)  IS  ' 

CALL  PRNT<V,NR,NSEN) 

CALL  LSVDF ( V , NCOL , NR , NSEN , TEN , NCOL , - 1 , SING , STOR , lER ) 
PRINT*,'  ' 

PRINT*,'  LSVDF  K3  lER  =  ',IER 
PRINT' (//)' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT(V, NSEN, NSEN) 

P3  =  NSEN  -  NR 
IF  (P3.LT.1)  THEN 
DO  510  1=1, NSEN 

510  GAMMA3(I,1)  =  V(I,NSEN) 

P3  =  1 

ELSE 

DO  511  1=1, NSEN 
DO  511  J=1,P3 

511  GANHA3(I,J)  =  V(I,J*NR) 

ENDIF 

PRINT*,'  TRANSFORMATION  NATRIX  GAMHA3  ' 

CALL  PRNT (GANNAS, NSEN, P3) 

CHECK  TO  SEE  THAT  GANNAS  IS  ORTHOGONAL  TO  C4 


o  o  n  o  o 


CALL  MMUL ( CT , G AMNA3 , NR2 , NSEN , P3 , AKC ) 

PRINT*,'  C4T  •  GAHMA3  ' 

CALL  PRNT(AKC,NR2,P3> 

C 

PRINT*,'  C4  SINGULAR  VALUES  ' 

CALL  PRNT(SING,NR,1) 

CALL  TFR(TRT,GAHNA3,NSEN,P3,1,2) 

CALL  NNUL ( TRT , GANHA3 , P3 . NSEN . P3 , RK> 

CALL  GMINV(P3,P3,RK,RK3,J,TAPE) 

CALL  TFR(CT,C3, NSEN, NC32, 1.2) 

CALL  NHUL ( TRT , C3 . P3 , NSEN . NC32 , AKC ) 

CALL  NNUL ( CT , GANNA3 , NC32 , NSEN , P3 , K0B3 ) 

CALL  NNUL ( K0B3 , RK3 , NC32 , P3 . P3 , STOR ) 

CALL  NNUL  ( STOR ,  AKC ,  NC32 ,  P3 ,  NC32 ,  CTCC3  > 

ENDIF 

C 

C 

C  CTCC3  WILL  BE  SUBSTITUTED  BACK  INTO 
C  NRIC-SYSTEN  3  FOR  A  NEW  K2  (WHEN 
C  USING  FOUR  CONTROLLERS) 

C 

C 

C  NOW  THE  CONTROLLER  GAIN  NATRIX,  G2 

C 

C 

DO  512  1*1, NCI 
DO  512  J»1,NACT 

512  V<I,J)  «  B1<I*NC1,J) 

PRINT*,'  V  (Bl)  IS  ' 

CALL  PRNT<V,NC1,NACT) 

CALL  LSVDF ( V , NCOL , NCI , NACT , TEN , NCOL , - 1 , SING , STOR , lER ) 
PRINT*,'  ' 

PRINT*,'  LSVDF  CONTROL  2  lER  =  ',IER 
PRINT' (//)' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT(V, NACT, NACT) 

E2  =  NACT  -  NCI 
IF  (E2,LT.l)  THEN 
DO  513  1=1, NACT 

513  T2<I,1)  =  Vd.NACT) 

E2  =  1 

ELSE 

DO  514  1=1, NACT 
DO  514  J=1,E2 

514  T2(I,J)  =  V(I,J*NC1) 

ENDIF 

PRINT*,'  TRANSFORMATION  MATRIX  T2  ' 

CALL  PRMT(T2,MACT,E2) 

CHECK  TO  SEE  THAT  T2  IS  ORTHOGONAL  TO  Bl 

NOTE:  IN  THIS  SECTION,  BCG  IS  THE  WORK  AREA 

FOR  B  *  T.  IN  ALL  CASES  THESE  SHOULD 


aoo  noooooocio  o  o  ooo 


BE  BLOCK  ZERO  MATRICES. 


CALL  MMUL ( B1 , T2 , NC12 , NACT , E2 , BCG  > 

PRINT*,'  B1  •  T2  ' 

CALL  PRNT(BCG,NC12,E2) 

PRINT*,'  B1  SINGULAR  VALUES  ' 

CALL  PRNT (SING, NCI, 1) 

CALL  VHULFN ( T2 , T2 , NACT , E2 , E2 , NACT , NACT , RK , N ACT , I ER ) 

CALL  GHINV(E2,E2,RK,RG2,J,TAPE> 

CALL  MMUL ( B2 , T2 , NC22 , N ACT , E2 , K0B2  > 

CALL  MMUL ( K0B2 , RG2 , NC22 , E2 , E2 , SAT2> 

CALL  VMULFP ( SAT2 , T2 , NC22 , E2 , NACT , NCOL , NACT , K0B2 , NCOL , lER ) 
CALL  VMULFP ( K0B2 , 82 , NC22 , NACT , NC22, NCOL , NCOL , SAT2 , NCOL , lER  > 


THIS  SAT2  WILL  BE  SUBSTITUTED  BACK  INTO  MRIC 
SYSTEM  2  FOR  A  NEW  G2. 


NOW  THE  CONTROLLER  GAIN  MATRIX,  G3 


DO  515  1=1, NCI 
DO  515  J»1,NACT 

515  V(I,J)  »  B1<I*NC1,J) 

DO  516  1=1, NC2 

DO  516  J=1,NACT 

516  V<I*NC1,J)  =  B2(I*NC2,J) 

V  IS  NOW  AN  NC1*NC2  BY  NACT  MATRIX 

PRINT*,'  V  <B1/B2)  IS  ' 

NRV  =  NCI  ♦  NC2 
CALL  PRNT <V, NRV, NACT) 

CALL  LSVDF ( V , NCOL , NRV , NACT , TEN , NCOL, - 1 , SING , STOR , lER ) 
PRINT*,'  ' 

PRINT*,'  LSVDF  CONTROL  3  lER  =  ',IER 
PRINT' <//)' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT (V, NACT, NACT) 

E3  =  NACT  -  NRV 
IF  (E3.LT.1)  THEN 
DO  517  I  =1,NACT 

517  T3(I,1)  =  V(I,NACT) 

E3  =  1 

ELSE 

DO  518  1=1, NACT 
DO  518  J=1,E3 

518  T3(I,J)  =  V(I,J*NRV) 

ENDIF 
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PRINT*,'  TRANSFORMATION  MATRIX  T3  ' 

CALL  PRNT{T3,NACT,E3) 

C 

C  CHECK  TO  SEE  THAT  T3  IS  ORTHOGONAL  TO  B1  AND  B2 
C 

CALL  MMUL ( B1 . T3 . NC12 , NACT , E3 , BCG ) 

PRINT*,'  B1  *  T3  ' 

CALL  PRNT(BCG,NC12,E3) 

CALL  MMUL  ( B2 ,  T3 ,  NC22 ,  N ACT ,  E3 ,  BCG) 

PRINT*,'  B2  *  T3  ' 

CALL  PRNT(BCG.NC22,E3) 

C 

PRINT*,'  B12  SINGULAR  VALUES  ' 

CALL  PRNT(SING,NRV,1) 

C 

CALL  VMULFM ( T3 , T3 , N ACT , E3 , E3 , NACT , NACT , RK , N ACT , I ER ) 

CALL  GNINV(E3,E3,RK,RG3,J,TAPE> 

CALL  MMUL<B3,T3,NC32, NACT, E3,K0B3) 

CALL  MMUL ( X0B3 , RG3 , NC32 , E3 , E3 , SAT3) 

CALL  VMULFP ( S AT3 , T3 , NC32 , E3 , N ACT , NCOL , N  ACT , K0B3 , NCOL , I ER ) 
CALL  VMULFP ( K0B3 , B3 , MC32 , NACT , NC32 , NCOL , NCOL , SAT3 , NCOL , lER ) 
C 
C 

C  SAT3  WILL  BE  SUBSTITUTED  BACK  INTO 
C  MRIC-SYSTEM  3  FOR  A  MEW  G3. 

C 

C 

C  NOW  THE  CONTROLLER  GAIN  MATRIX,  G4, 

C  TO  BE  FOUND  WHEN  USING  FOUR  CONTROLLERS 

C 

C 

IF  (DEC. EG. 4)  THEN 
DO  519  1=1, NCI 
DO  519  J=1,NACT 

519  V(I,J)  =  B1(I*NC1,J) 

DO  520  1=1, NC2 

DO  520  J=1,NACT 

520  V(I*NC1,J)  =  B2<I*NC2,J) 

DO  521  1=1, NC3 

DO  521  J=1,NACT 

521  V(I*MC1*NC2,J)  =  B3(I*NC3,J) 

C 

C  V  IS  NOW  A  <NC1»MC2*NC3)  BY  NACT  MATRIX 
C 

PRINT*,'  V(B1/B2/B3)  IS  ' 

NRV  =  NCI  ♦  MC2  ♦  NC3 
CALL  PRNT<V, NRV, NACT) 

CALL  LSVDF < V, NCOL , NRV , NACT , TEN , NCOL, - 1 , SING , STOR , lER > 
PRINT*,'  ' 

PRINT*,'  LSVDF  CONTROL  lER  =  ',IER 
PRINT' <//)' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNTCV, NACT, NACT) 


ooooo  ooooooo  o  o  ooo 


Ml'*  ■ "  1«M  ■■  I  ng 


DO  522  I=1,NACT 

522  T4(I,1)  =  Vd.NACT) 

E4  1 

ELSE 

DO  523  I=1,NACT 
DO  523  J-1.E4 

523  T4<I,J)  =  V<I,J*NRV> 

ENDIF 

PRINT*,'  TRANSFORMATION  MATRIX  T4  ' 

CALL  PRMT(T4,MACT,E4) 

CHECK  TO  SEE  THAT  T4  IS  ORTHOGONAL  TO  B1,B2,B3 

CALL  NMUL ( B1 , T4 , NC12 , NACT , E4 , BCG ) 

PRINT*,'  B1  *  T4  ' 

CALL  PRNT(BCG,NC12,E4) 

CALL  MMUL ( B2 , T4 , NC22 , N ACT , E4 , BCG ) 

PRINT*,'  B2  *  T4  ' 

CALL  PRNT(BCG,NC22,E4) 

CALL  MMUL ( B3 , T4 , NC32 , N ACT , E4 , BCG ) 

PRINT*,'  B3  *  T4  ' 

CALL  PRNT(BCG,NC32,E4) 

PRINT*,'  B123  SINGULAR  VALUES  ' 

CALL  PRNT(SING,NRV,1> 

CALL  VMULFM  <  T4 , T4 , NACT , E4 , E4 , NACT , NACT , RK , NACT , lER ) 

CALL  GHINV(E4,E4,RK,RG4,J,TAPE> 

CALL  NMUL ( B4 , T4 , NR2 , NACT , E4 , K0B4  > 

CALL  MMUL ( K0B4 , RG4 , NR2 , E4 , E4 , SAT4 > 

CALL  VMULFP  <SAT4 , T4 ,  NR2 ,  E4 , NACT ,  NCOL ,  NACT ,  KOB4 ,  NCOL ,  lER ) 
CALL  VMULFP ( KOB4 , B4 , NR2 , NACT , NR2 , NCOL , NCOL , SAT4 , NCOL , lER ) 
ENDIF 


SAT4  WILL  BE  SUBSTITUTED  BACK  INTO  MRIC- 
SYSTEM  4  FOR  A  NEW  G4  WHEN  USING  FOUR 
CONTROLLERS. 


Z2  =  1 
GOTO  H5 
600  CONTINUE 


THE  PROBLEM  IS  NOW  COMPLETE 


PRINT'<///)' 

PRINT*,'  THIS  RUN  HAS  BEEN  COMPLETED  ' 
PRINT' (///)' 


END 

SUBROUTINE  FORMXO(XO, INIT.DEC) 

CONMON/NAINA/NDA 

COHMON/MAINB/NCOL 

C 

C 

C  FORMXO  IS  THE  SAME  FOR  THREE  OR  FOUR 
C  CONTROLLERS  SINCE  THE  RESIDUALS  BECOME 
C  THE  FOURTH  CONTROLLER. 

C 

C0MM0N/NUM/IC1<21),IC2(21),IC3<21),IR(21),NC1,MC2,NC3,NR 
REAL  X0(N0A),INIT(4.NC0L> 

INTEGER  N,I,J,K,L 
DO  1  1  =  1,  MCI 
N  =  ICKI) 

XO<I)  =  1NIT<1,M) 

XO(I*NCl)  =  INIT(2,M) 

XO<I*NCl*2)  =  IMIT(3,M> 

1  XO<I*NCl*3)  =  IMIT<4,M) 

J  =  NCl-4 

DO  2  1=1, NC2 
N  =  IC2<I) 

XO<I*J)  =  INIT(1,M) 

X0<I*J+NC2)  =  INIT(2,M) 

X0a*J*NC2*2)  =  INIT(3,N) 

2  XO(I*J-MC2»3)  =  IMIT(4,M> 

K  =  J  ♦  NC2«4 

DO  3  1=1, NC3 
M  »  IC3(I) 

XO<I+K)  =  INIT<1,M) 

X0(I*K*NC3)  =  INIT(2,N) 

X0<I*K*NC3«2)  =  INIT(3,M> 

3  X0<I*K*NC3»3)  =  INIT(4,M) 

IF  (NR.ME.O)  THEN 

L  =  K  ♦  NC3*4 
DO  4  1=1, NR 
H  =  IR(I) 

XO<I*L)  =  INIT<1,M) 

XO<I+L*NR)  =  INIT(2,M) 

IF  (DEC.EQ.4)  THEN 
XO{I*L«-NR»2)  =  INIT(3,M) 

X0<I»L+NR»3)  =  INIT(4,M) 

ENDIF 

4  CONTINUE 
ENDIF 
RETURN 
END 

SUBROUTINE  FACTOR(N, A,S,MR) 

C  A=S'S 

DIMENSION  A(1),S(1} 

COHMON/MAINB/  NCOL.NCOLl 

COMMON/INOUT/KOUT 

T0L=l.E-6 
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T0L1=0. 

DO  1  I-l,NN,NCOLl 
R-ABS(A(I)) 

1  IF  (R.GT.TOLl)  TOLl=R 

T0L1»T0L1»1.E-12 
11*1 

DO  SO  1*1, N 
IM1*I-1 

DO  5  JJ>I,NN,NCOL 
5  S(JJ)*0. 

ID=II^IM1 

R*A<ID)-DOT(II!l,S<II),S{II)) 

IF  (ABS(R).LT.<TOL»A<ID)*TOLl))  GO  TO  SO 
IF  (R)  IS, SO, 20 
IS  MR*-1 

WRITE (KOUT, 1000) 

1000  FORMAT (37H0TRIED  TO  FACTOR  AN  INDEFINITE  MATRIX 
RETURN 

20  S(ID)*SQRT(R) 

MR=MR>1 

IF  (I.EQ.N)  RETURN 

L*II>NCOL 

DO  25  JJ*L,NN,NCOL 

IJ*JJ^IN1 

25  S<IJ)»(A(IJ>-D0T<IM1,S(II),S(JJ)>)/S(ID> 

50  llall^NCOL 
RETURN 
END 

SUBROUTINE  F0RNA(A,0,W,N,N2,IC> 

COMMON/MAINB/NCOL 

REAL  A(NCOL,NCOL) ,W(NCOL) ,D<NCOL} 

INTEGER  IC(N),I,J,N,M 
DO  1  1*1, N2 
DO  1  Jal,N2 
A(I,J)=0.0 
L  CONTINUE 

DO  2  1*1, N 
M=  IC(I) 

A<<I*N>,(I*N))*D<M) 

A(I,<I»N))  =  1.0 
A(<I*N),I)  =  -(W(H)*«2) 

>  CONTINUE 

RETURN 
END 

SUBROUTINE  FORMB ( B , PHI , N , N2 , N ACT , IC ) 

COMMON/MAINB/NCOL 

REAL  B(NCOL,NCOL)  ,PHI(NCOL,NCOL> 

INTEGER  IC(N),NACT,N,M,I,J 
DO  1  1*1, N2 
DO  1  J»1,NACT 
B(I,J)  *  0.0 
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1  CONTINUE 
DO  2  1=1, N 
N  =  IC(I) 

DO  2  J=1,NACT 
B(<N^I),J)  =  PHI<M,J> 

2  CONTINUE 
RETURN 
END 

SUBROUTINE  FORMC ( C , PHIS , N . N2 , NSEN , IC ) 

CONNON/MAINB/NCOL 

REAL  C(NCOL,NCOL) ,PHIS(NCOL,NCOL> 

INTEGER  IC(N),N,NSEN,N,N2.I,J 

DO  1  1=1, NSEN 

DO  1  J=1,N2 

C(I,J)  =  0.0 

1  CONTINUE 

DO  2  1=1, NSEN 
DO  2  J=1,N 
M  =  IC(J) 

C(I,J)  =  PHIS(H,I> 

2  CONTINUE 
RETURN 
END 

SUBROUTINE  F0RNQ1(Q,A,N,IC) 

COMHON/MAINB/NCOL 

REAL  A<NCOL),Q(NCOL,NCOL) 

INTEGER  I,J,K,N,N,N2,IC(NC0L) 

N2  =  N  •  2 
DO  1  1=1, N2 
DO  1  J=1,N2 
Q(I,J)  =  0.0 

1  CONTINUE 
DO  2  1=1, N 
N  =  IC(I) 

0(1,1)  =  A(H) 

Q(I*N,I*N)  =  0(1, I) 

2  CONTINUE 
RETURN 
END 

SUBROUTINE  GH INV ( NR , NC , A , U , NR , NT ) 
DINENSION  A(1),U(1) 

COHNON/NAINl/  NDIN,NDIN1,S<1) 

C0NH0N/NAINB/NC0L,NC0L1 

CONNON/INOUT/KOUT 

T0L=1.E-12 

NR=NC 

NRN1=NR-1 

T0Ll=l.E-20 

JJ=1 

DO  100  J=1,NC 
FAC=DOT(NR,A(JJ) , A<JJ) ) 

JN1=J-1 

JRM=JJ*NRM1 
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DO  20  I=JJ,JCM 
U(I)»0. 

U<JC1!)  =  1.0 

IF  (J.EQ.l)  GO  TO  54 
KK>1 

DO  30  K>1,JN1 

IF  (S(K).EQ.l.O)  GO  TO  30 

TEHP>-DOT(NR,A(JJ),A(KK)) 

CALL  VAOD(K,TEMP,U(JJ),U(KK)> 

KK«KK«NCOL 

DO  50  L>1.2 

KK>1 

DO  50 

IF  (S<K).EQ.&.)  GO  TO  50 
TEIIP»-DOT<IIR,A<JJ),A(KK)) 

CALL  VADD(NR,TEMP,A(JJ),A(KK)) 

CALL  VADO(K,TEMP.U<JJ),U(KK)) 

KK^KK^NCOL 

T0L1»T0L»FAC 

FAC>DOT<NR,A<JJ).A(JJ)) 

IF  ( FAC. GT. TOLD  GO  TO  70 

DO  55  I»JJ,JRM 

A<I)>0. 

S(J)«0. 

KK-1 
DO  65 

IF  (S(K).EQ.O.)  GO  TO  65 
TEMP»-DOT(K,U<KK) ,U<JJ)) 

CALL  VADD(MR,TEHP.A(JJ>,A<KK)> 
KK^KK^NCOL 

FAC»DOT(J,U<JJ),U<JJ)) 

HR=I!R-1 
GO  TO  75 
S<J)sl.O 
KK*1 

DO  72  K»1,JK1 
IF  (S(K).EO.l.)  GO  TO  72 
TEHP=-OOT<MR,A(JJ) ,A<KK) ) 

CALL  VADD<K,TEKP,U<JJ),U(KK)) 

KK-KK»NCOL 

FAC=1. /SORT (FAC) 

DO  80  I=JJ,JRN 
A(I)=A<I)*FAC 
DO  85  I=JJ,JCM 
U(I)=U<I)*FAC 
JJsJJ*NCOL 

IF  (NR.EQ.NR.OR.HR.EQ.NC)  GO  TO  120 
IF  (MT.ME.O)  WRITE<K0UT,110)MR.MC,MR 
F0RNAT(I3,1HX,I2,8H  H  RANK, 12) 
NEND>NC>NCOL 
JJ=1 

DO  135  J«1,NC 
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DO  125  1*1, NR 

II«I-J 

S<I>*0. 

DO  125  KK*JJ,NEND,NCOL 
125  S(I)=S<I)^A(II*KK)*U(KK) 

II*J 

DO  130  1*1, NR 
U(II>*S(I) 

130  II*II«NCOL 
135  JJ»JJ*NC0L1 
RETURN 
END 

SUBROUTINE  INTEG(N,A,C,S,T)  ' 

:  S* INTEGRAL  EA»C>EA  FROM  0  TO  T 

:  C  IS  DESTROYED 

DIMENSION  A(1),C(1),S(1> 
COMHON/MAINl/  NDIM,NDIH1,  X(l) 
C0NM0N/NAINB/NC0L,NC0L1 
C0HN0N/NAIN2/C0EF(100) 

NN*N>NCOL 

NH1*N-1 

IN0*0 

ANORM*XNORM(N,A) 

DT*T 

5  IF  <ANORM«ABS(DT).LE.O.S)  GO  TO  10 

DT=DT/2. 

IND»IND*1 
GO  TO  5 

10  DO  15  I*1,MN,NC0L 
J»I^NM1 
DO  15  JJ=I,J 
15  S(JJ)=DT*C<JJ) 

Tl=DT*»2/2. 

DO  25  IT* 3, 15 

CALL  NMUL(A,C,N,N,N,X) 

DO  20  1=1, N 
II=(I-1)*NC0L 
DO  20  JJ*I,NN,NCOL 
11*11*1 

C<JJ)=<X(JJ)*X(II))»T1 
20  S(JJ)*S<JJ)*C(JJ) 

25  T1=DT/FL0AT<IT) 

IF  (IND.EQ.O)  GO  TO  100 
C0EF(11)=1.0 
DO  30  1*1,10 
11=11-1 

30  C0EF(II)=DT*C0EF(II*1)/FL0AT<I) 
11*1 

DO  40  I*1,NN,NC0L 
J*I*NM1 
DO  35  JJ*I,J 
35  X<JJ)*A(JJ)«C0EF(1) 

X(II)*X(II)*C0EF(2) 
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40  I1=II*MC0L1 
DO  55  L=3,ll 
CALL  MMUL(A.X,N,N.N,C) 

Tl=COEFCL) 

DO  55  I»1,NN,NC0L 

DO  50  JJsI.J 
50  XCJJ)»CCJJ> 

X(II)=X<II)*T1 
55  II»II*MC0L1 
C  X=EXP(A*DT) 

L»0 

60  L»L*1 

CALL  HMUL(X,S,N,N,N,C) 

11=1 

DO  90  1=1, N 
J=II 

IF  (I.EQ.l)  GO  TO  75 
DO  70  JJ=I,II,NCOL 
S(JJ)=S(J) 

70  J«J*1 

75  DO  85  JJ=I,M 
KK=JJ 

DO  80  X>I,NN,NCOL 
S<J)»S<J)*C(K)*X(KK) 

80  KK=KK«NCOL 
85  J»J*KCOL 

DO  87  JJ»I,NN,NCOL 
87  C<JJ)=X(JJ) 

90  II  =  II*IICOL 

IF  (L.EQ.IHD)  GO  TO  100 
CALL  NI1UL(C,C,N,N.N,X) 

GO  TO  60 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  NEXP(N,A,T,EA) 

DIMENSION  Ad)  ,EA(1)  ,C(84>  ,0(85)  ,E(84> 
COnMON/NAINl/NDIH,NDINl,TEN,X<l) 

NN  =  N«*2 
MMl  =  N-1 
NPl  =  N+1 

IF  (N.GT.l)  GO  TO  5 
EA(1)=EXP(T*A(1)) 

RETURN 
5  W=1.0 

DO  10  I=1,NN.N 
IL=I«NH1 
DO  10  J=I,IL 
10  EA<J)*A<J> 

C1=XN0RH1(N,A} 

IND=0 


T1=T 

15  IF  (ABS<T1»C1).LE.3.0)  GO  TO  20 
Tl=Tl/2. 

IND>IND«1 
GO  TO  IS 
20  C2>0. 

DO  25  I«1,NN,NP1 
25  C2»C2-EA(I> 

C2»C2/FL0AT(L) 

C(L>»C2 

D{L>1)»0. 

EdD-U 

11*1 

DO  35  I*1,NN.N 
IL»I*MM1 
DO  30  J=I,IL 
30  X(J)=EA(J) 

X(II)*X<II>*C2 
35  1I»H*MP1 

IF  (L.EQ.N)  GO  TO  40 
CALL  HMUL1(X,A,N,N,N,EA> 

W»W*T1 /FLOAT <L) 

L=L*1 
GO  TO  20 
40  COMTIMUE 

****  CAN  CHECK  X  0  FOR  ACCURACY  (63  IS  OVERKILL  FOR  CONVERGENCE) 
J»N^63 
DO  50  L*N,J 
DO  45  K*1,N 

D<K)s(D<K*l)-W*C(K))*Tl/FLOAT<L) 

45  E(K)»E<K)*D<K) 

50  W*D(1) 

11=1 

DO  60  I=1,NN,N 
IL=I*NM1 
DO  55  J=I,IL 
55  EA(J)=E(1)«A<J) 

EA(II)=EA<II)«E<2} 

60  II=I1*NP1 

IF  (N.Ea.2)  GO  TO  85 
DO  80  L=3,N 

CALL  NMUL1(EA,A,N,N,N,X) 

11  =  1 
C2=E(L) 

DO  75  I=1,NN,N 
IL=I«NN1 
DO  70  J»I,IL 
70  EA(J>»X(J> 

EA<II)»EA(II)*C2 
75  II«II*NP1 
80  CONTINUE 
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dS  IF  (IND.EQ.O)  RETURN 
DO  100  L>1,IND 
DO  90  I<1.NN,N 
IL»I*NII1 
DO  90  JsI.lL 
90  X(J)»EA(J) 

100  CALL  NHUL1(X,X,N,N,N,EA) 

RETURN 
END 

SUBROUTINE  NLINEQ ( N , A , C , X , TOL , lER ) 
SOLVES  A'X^XA*C»0 
A  AND  X  CAN  BE  IN  SANE  LOCATION 
ANSWER  RETURNED  IN  C  AND  X 
DIHENSION  A(1).C(1),X(1) 
COHNON/HAINB/  NCOL,  NCOLl 
C0HH0N/HAIN3/F(1) 

ADV=T0L*l.E-6 
DT=.5 
DT1»0. 

NN-N»NCOL 

DO  5  II-1,NN»NC0L1 
5  DT1»DT1-A(II) 

DT1*DT1/N 

IF  <DT1.GT.4.0)  DT=DT»4.0/DT1 
II=»1 

DO  20  lal.N 
DO  IS  JJ*I,NN,NCOL 
15  X(JJ)»DT»A<JJ) 

X(II)*X<Il)-.5 
20  II»II ♦NCOLl 

CALL  GMINV(N,N,X,F,NR,0) 

IER»4 

IF  (NR.NE.N)  RETURN 
CALL  HNUL(C,F,N,N,N,X) 

C  INITIALIZATION  OF  X,F 
1  =  1 

DO  40  II=1,NN,NC0L 
J=II 


G 

25 

30 

^  35 

40 

G 


IF  (I.EQ.l)  GO  TO  30 
DO  25  JJ=I,II,NCOL 
C(J)=C(JJ) 

J»J*1 


ID=J 

DO  35  JJ=II,NN,NCOL 

C(J)=DT*OOT(N,F<II),X<JJ)) 

J=J*1 

F<ID)=F(ID)^1.0 

I=I+1 

DO  90  IT=1,20 
NEZ>0 

CALL  MNUL(C,F,N,N,N,X) 

1=1 

II»1 
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J=1 

GO  TO  70 
60  JslI 

DO  65  JJsI,II,NCOL 
C(J)=C(JJ) 

65  J»J*1 

70  IDsJ 

DT1=C(J) 

DO  75  JJ»II,NN,NCOL 
C(J)*C(J)*DOT(M,F(II),X<JJ)) 

75  J=J^1 

J*J-1 
DO  80 

80  X(JJ)»F(JJ) 

IF  (ABS(C(ID)).GT.1.E150)  GO  TO  95 

IF  (ABS(C(ID)-DT1).LT.<ADV*T0L»ABS(C(ID>))>  ME2=MEZ+1 
1=1*1 

II=II*MCOL 

IF  (I.LE.N)  GO  TO  60 
IF  (NEZ.EQ.N)  GO  TO  150 
CALL  NMUL(X,X,N,N,N,F) 

90  CONTINUE 
95  IER=1 
RETURN 

150  CONTINUE 
NN1=N-1 

DO  155  I=1,NN,NC0L 
rr»i*NNi 
DO  155  JJ=I,II 
155  X<JJ)»C<JJ) 

IER=0 

RETURN 

END 

SUBROUTINE  NNUL(X, Y,N1 ,N2,N3,Z) 

CONHON  /NAINB/  NCOL 

DIMENSION  X ( NCOL , 1 ) , Y ( NCOL , 1 ) , 2 ( NCOL , 1 ) 

DO  3  J=1,N3 
DO  2  1=1, N1 
S=0. 

DO  1  K=1,N2 

1  S=S*X(I,K)*Y<K,J) 

2  Z(I,J)=S 

3  CONTINUE 
END 

SUBROUTINE  HHUL1(X,Y,N1,N2,N3.Z) 

REAL  X(1),Y(1),Z(1) 

COHNON/NAINA/NDA 

NEND3=NDA>N3 

NEND2=NDA»N2 

DO  1  1  =  1, N1 

DO  1  J=1,NEND3,NDA 

TM=0. 

K=I 


B-37 


o  o  o 


5 


KK>KK«1 

TM=TM*X<K)«Y(KK) 

K>K«NDA 

IF  (K.LE.NEND2)  GOTO  S 
1  Z(J)»TI1 
RETURN 
END 

SUBROUTINE  NRIC(N.A,S,Q.X,Z,TOL. lER) 
DIMENSION  A(l),S<l),a(l).X<l),Z(l> 
CONNON/HAINl/NDIN,  NDIMl,  F(l> 
CONHON/NAINB/NCOL.NCOLl 
C0NH0N/NAIN2/TR(1) 

COMMON/ INOUT/KOUT 

ADV=T0L*l.E-6 

NN«N«NCOL 

NMlsN-1 

IND>1 

C0UNT=0. 

IF  (lER.EO.l)  C0UNT»99. 

IF  (lER.EQ.l)  NR>N 
IF  (lER.EQ.l)  GO  TO  100 
Tl»-1. 

300  CONTINUE 
IER«0 

C0UNT»C0UNT>1. 

DO  IS  I>1,N 
DO  IS  J>I,NN,NCOL 

15  X(J)»-S(J) 

CALL  INTEG(N,A,X,Z,T1) 

CALL  FACTOR (N,Z,X, NR) 

IER=1 

IF  (NR.LT.O)  GO  TO  200 
IER=0 

CALL  GNINV(N,N,X,Z,NR,0) 

CALL  TFR(TR,Z,N,N,1,2) 

CALL  NMUL(Z,TR,N,N,N,X) 

DO  IS  II-1,NN,NC0L1 
I  =  II 

DO  17  J=II,MN,NCOL 
X(J)=<X(J)*X(I))/2. 

X(I)=X(J) 

17  I=I>1 

18  CONTINUE 
100  CONTINUE 

DO  16  I-1,N 

16  TR<I)=-1.0 
A«SX  IS  STABLE 

POSSIBLE  UNCONTROLLABILITY  IF  NR.NE.N 
JIN  DILLOW'S  REPUTATION  PRECEDES  HIM 
T0L1»T0L/10. 

HAXIT»40 
DO  40  IT«1,MAXIT 
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IF  (lER.EQ.l)  GO  TO  101 
CALL  NMUL(S,X.N,N,N.F) 

CALL  HMUL(X,F,N.N.N,Z) 

DO  20  I»1,NN,NC0L 
II«I*NM1 
DO  20  J=I,II 
X(J)«A(J)-F<J) 

20  Z<J)*2(J)^0(J) 

101  CONTINUE 
IER>0 

CALL  HLINE0<N,X.Z,X,T0L1.IER) 

IF  (lER.NE.O)  GO  TO  200 

L=0 

C1»0.0 

II»1 

DO  25  I»1,N 

IF  <ABS<X<II)-TR(I)).LT.(ADV*TOL»X(II)))  L=L»1 
TR(I>»X(II) 

II*II*NC0L1 
25  C1=C1*TR<I) 

IF  <ABS(C1).GT.1.E20)  GO  TO  50 
IF  (L.NE.N)  GO  TO  40 
CALL  GMINV(N,N.2.F,NR,0) 

CALL  MNUL(S,X.N.N,N,Z) 

DO  30  I>1,NN,NC0L 
II=I*NM1 
DO  30  J»I,II 
30  2(J)»A<J)-Z<J) 

IF  (NR.NE.N)  WRITE(K0UT,3S)NR 
35  FORMAT (27H0RICCATI  SOLN  IS  PSD--RAMK  ,I3> 

GO  TO  65 
40  CONTINUE 

WRITE <KOUT, 45)  MAXIT 

45  FORMAT (27H0RICCATI  NON -CONVERGENT  IN  ,I2,11H  ITERATIONS) 

GO  TO  60 

50  WRITE<K0UT,55)IT,T1 

55  FORMAT OOHORICCATI  BLOW-UP  AT  ITERATION  ,I2,12H  INITIAL  T=  ,F10.5) 
60  IER=1 
65  RETURN 

200  IF  <IND.EQ.2)  GO  TO  250 
IF  ( COUNT. GE. 10.)  RETURN 
Tl=Tl/<2. ••COUNT) 

IND>2 
GO  TO  300 

250  T1=T1^<2.^^C0UNT) 

IND-1 
GO  TO  300 
END 

SUBROUTINE  PRNT(MAT,N,N) 

CONNON/MAINB/NCOL 
REAL  NAT(NCOL.NCOL) 

INTEGER  N,I,J,K,M 
PRINTS,'  ' 


▼ 


w 


\  # 


t* 


c 
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IF  (M.GT.12)  GOTO  2 
DO  1  1=1, N 

PRINT' (1X,12F10.4)' ,(MAT(I,J),J=1,M) 

1  CONTINUE 
GOTO  10 

2  CONTINUE 

IF  (M.GT.24)  THEN 
CALL  PRNTXL(NAT,N,N) 

RETURN 
ENDIF 
DO  3  1*1, N 

PRINT' (1X,12F10.4> ' , <MAT<I,J) ,J=1,12) 

3  CONTINUE 
PRINT' <//)' 

DO  4  1*1, N 

PRINT' <1X,12F10.4) ' , <HAT(I, J) ,J=13,M) 

4  CONTINUE 

10  PRINT' <///)' 

RETURN 

END 

SUBROUTINE  PRNTXL(HAT,N,M) 
COHHON/HAINB/NCOL 
COMHON/NAINA/NDA 
REAL  NAT(NDA,NOA) 

INTEGER  I,J,K,L,N,N 
PRINT*,'  ' 

DO  1  L*1,N,12 
K  »  L  ♦  11 

IF  (H-L.LT.ll)  K  =  N 
DO  2  1*1, N 

PRINT' (1X,12F10.5)' , (HAT(I,J) ,J=L,K) 

2  CONTINUE 

PRINT' (//)' 

1  CONTINUE 
PRINT' (///)' 

RETURN 

END 

SUBROUTINE  PRNTSH<NAT,N,H) 

CONHON/NAINB/NCOL 

REAL  NAT(3,NC0L) 

INTEGER  I,J,K,L,M,N 
PRINT*,'  ' 

DO  1  L*1,H,12 
K  =  L  ♦  11 

IF  (H-L.LT.ll)  K  =  M 
DO  2  1*1, N 

PRINT' (1X,12F10.5) ' , (HAT(I, J) ,J=L.K) 

2  CONTINUE 
PRINT' <//)' 

1  CONTINUE 

PRINT' (///)' 

RETURN 

END 
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SUBROUTINE  TFR<X,A,N,M,K,I> 

I  -  1  GIVES  X  =  A 

2  GIVES  X  =  A' 

3  GIVES  X  =  A  AS  A  VECTOR 

4  GIVES  A  =  X  WHERE  X  WAS  A  VECTOR 

DIMENSION  X(1),A(1} 

COHMON/NAINB/NCOL 

JS=<K-l)*NCOL«M 

JEND»H>NCOL 

GO  TO  (10, 30, 50, 70), I 

DO  20  11=1, N 

DO  20  JJ=II,JEND.NCOL 

X(JJ)=A(JJ*JS) 

RETURN 

DO  40  11  =  1, N 
KK=(II-1)»NC0L 
DO  40  JJ=1,M 
LL=<JJ-1)*NC0L*II 
X<KK*JJ)=A(LL*JS) 

RETURN 

KK=0 

DO  60  II=1,JEND,NC0L 

LL=II>M-1 

DO  60  JJ=I1,LL 

KK=KK^1 

X<KK)-A(JJ*JS) 

RETURN 

KK=M»N*1 

DO  80  11=1, M 

LL=<«-II)»NC0L*1 

DO  80  IJ=1,N 

KK=KK-1 

JJ=LL*N-IJ 

A<JJ*JS)=X(KK) 

RETURN 

END 

FUNCTION  DOT (NR, A, B) 

DIMENSION  A(1),B(1) 

D0T=0. 

DO  1  1=1, NR 
DOT=DOT*A(I)*B(I) 

RETURN 

END 

SUBROUTINE  VADD(N,C1 ,A,B) 

DIMENSION  A(1).B(1) 

DO  1  1  =  1, N 
A(I)=A(I)*C1*B(I) 

RETURN 

END 

FUNCTION  XNORMKN.A) 

COMPUTES  AN  APPROXIMATION  TO  NORM  OF  A--  MOT  A  BOUND 


o  o  o 


DIMENSION  A(N*N) 

NN=N*N 
NPl  =  N*1 
C1=0. 

TR=A<1) 

IF  (N.EQ.l)  GO  TO  20 
1=2 

00  10  II=NP1,NN.N 
J=II 

DO  5  JJ=I,II,N 
C1=C1*ABS(A(J)«A(JJ)) 

5  J=J+1 

TR=TR*A(J) 

10  1=1*1 

TR=TR/FLOAT(N) 

DO  15  11=1, NN, NPl 
15  C1=C1*<A<II)-TR)**2 

20  XN0RM1=ABS(TR)*SQRT(C1) 

RETURN 

END 

FUNCTION  XNORH(N,A) 

C  COMPUTES  AN  APPROXIMATION  TO  NORM  OF  A  --  NOT  A  BOUND 
DIMENSION  A(l) 

COHMON/MAINB/NCOL, NCOLl 

NN=N>NCOL 

C1=0, 

TR»A<1) 

IF  (W.EO.l)  GO  TO  20 
1=2 

DO  10  II»NC0L1,NN,NC0L 
J=II 

DO  5  JJ=I.II,NCOL 
C1=C1*ABS(A(J)*A(JJ)) 

5  J=J*1 

TR=TR*A(J) 

10  1=1*1 

TR=TR/FLOAT(N) 

DO  15  11=1, NN, NCOLl 
15  C1=C1*(A(II)-TR)**2 

20  XN0RM=ABS(TR)*SQRT(C1) 

RETURN 
END 

SUBROUTINE  RFMPHIL ( PHIL , RPHIL , NCI , NC2 , NC3 , NR , ICl . IC2 . ICO , IR , 
*PETA,CPHIL,MM) 

COMMON/MAINA/NDA 
COMHON/MAINB/NCOL 

REAL  PHILO, NCOL) ,RPHIL<3,NC0L) ,CPHIL(3,NDA) ,PETA{NCOL,NDA) 
INTEGER  I,J,K,KK,L,M,MM,NC1,NC2,NC3,NR,NM0DE 
INTEGER  ICKNCOL)  ,IC2<NC0L)  ,IC3(NC0L)  ,IR(NCOL) 
NM0DE=NC1*NC2*NC3*NR 

FIRST  REFORM  PHILOS  lAW  MODE  NUMBER  REORDERING 


o  o  o 


L=0 

DO  1  1=1, NCI 
M  =  ICKI) 

DO  1  J=l,3 
RPHIL(J,I*L) 

1  CONTINUE 

L  =  L  ♦  NCI 
DO  2  1=1, NC2 
N  =  IC2(I) 

DO  2  J=l,3 
RPHIL<J,I*L) 

2  CONTINUE 

L  =  L  ♦  NC2 
DO  3  1=1, NC3 
N  =  IC3<I> 

DO  3  J=l,3 
RPHIL(J,I*L) 

3  CONTINUE 

IF  (NR.GT.O) 
L  =  L  ♦  NC3 
DO  4  1=1, NR 
M  =  IR(I) 

DO  4  J=1.3 
RPHIL(J,I>L) 

4  CONTINUE 
ENDIF 


=  PHIL(J,M) 


=  PHIL(J,N) 


=  PHIL(J,M) 
THEN 


=  PHIL(J.N) 


NOW  FORM  P(ETA)  TO  PICK  THE  ETAS  OUT  OF  THE  STATE  VECTOR 


DO  60  I=1,NM0DE 
DO  60  J=1,MM 

60  PETA(I,J)=0.0 
DO  61  1=1, NCI 
DO  61  J=1,NC1 

IF  (I.EQ.J)  THEN 

PETA{I,J)=1.0 

ENDIF 

61  CONTINUE 
K=NC1 
KK=4»NC1 

DO  62  1=1, NC2 
DO  62  J=1,NC2 
IF  (I.EQ.J)  THEN 
PETA(I*K,J*KK)=1.0 
ENDIF 

62  CONTINUE 
K=K*NC2 
KK=KK»4»NC2 
DO  63  1=1, NC3 
DO  63  J=1,NC3 

IF  (I.EQ.J)  THEN 

PETA(I*K,J*KK)=1.0 

ENDIF 


63  CONTINUE 

IF  (NR.EQ.O)  GOTO  65 

K=K*NC3 

KK=KK*4»NC3 

DO  64  I»1,NR 

DO  64  J=1,NR 

IF  (I.EQ.J)  THEN 

PETA{I*K,J+KK)=1.0 

ENDIF 

64  CONTINUE 

65  CONTINUE 

C  PRINT*,'  THE  MATRIX  P<ETA)  IS  ' 

C  PRINT*,'  ****TRANSPOSED****  ' 

C  PRINT' <//)' 

C  DO  66  J-1,HM 

C  66  PRINT' <1X,20F4.0)',(PETA<I,J),  1=1,20) 

C 

C  NOW  GET  CPHIL  =  RPHIL*PETA 
C 

CALLVHULFF(RPHIL,PETA,3,NH0DE.MH.3,NC0L,CPHIL,3,IER) 

RETURN 

END 

SUBROUTINE  TIMEX  < STM , MM , DT, XO , PDT , TMAX , XI , XC , CPHIL , XL) 
COHMON/MAINA/NDA 

REAL  DT , TM AX , ABT , STM ( NO A , ND A ) , XO  <  ND A ) , XI ( NDA ) , XC ( ND A ) 

REAL  CPHIL(3,NDA),XL(3) 

INTEGER  I, J, MM, PDT 

C 

C  THIS  ROUTINE  PROPAGATES  THE  STATE  VECTOR  IN  TIME  AND  OUTPUTS  LOS  DATA 
C  TO  BOTH  THE  MAIN  PRINTED  OUTPUT  FILE  AND  TO  A  PLOT  FILE.  IT  GIVES 
C  LOS  RADIUS  AND  OEFOCUS  VS  TIME  TO  THE  PLOT  FILE,  BUT  X,Y,2,R  AND  T 
C  TO  THE  PRINT  FILE. 

C 

ABT*0. 

C  MAKE  A  COPY  OF  THE  I.C.  VECTOR  AND  USE  THE  COPY. 

DO  10  1=1, MM 
XC(I)  =  XO(I) 

10  CONTINUE 

C  GET  THE  LINE-OF-SIGHT  DATA  OUTPUT  BY 
C  WRITING  IT  TO  THE  PRINTED  OUTPUT, 

C  AND  TO  TAPE?. 

WRITE(6, '(//)') 

WRITE(6,*)'  TIME  LOSX  LOSY  DEFOCUS  RADIUS 

WRITE(6,' (/)') 

C  CHECK  THAT  MAX  TIME,  'TMAX'.  IS  NOT  REACHED  AND  CONTINUE  TO 
C  PROPAGATE  THE  STATE  VECTOR  AND  GET  LOS  EVERY  'PDT'  TIME  INCREMENTS 
C  FIRST  OUTPUT  THE  INITIAL  CONDITION 

CALL  VMULFF  < CPHIL , XC , 3 , MM , 1 , 3 , NDA , XL , 3 , lER ) 

RAD  =  ((XL(1)**2>  ♦  <XL{2)**2))**0.5 
WRITE(6,1000)ABT,XL(1) ,XL(2) ,XL(3) ,RAD 
WRITE (7 , 1001 ) ABT , RAD , XL ( 3 ) 

C  NOW  BEGIN  PROPAGATION 
20  DO  40  1  =  1, PDT 


CALL  VNULFF ( STM , XC , HM , MM . 1 , NDA , NDA , XI . NDA , lER) 

DO  50  J=1,MM 
XC(J>  =  XKJ) 

50  CONTINUE 
40  CONTINUE 

C  INCREMENT  THE  ABSOLUTE  TIME  AND  CALCULATE  LOS  FOR  OUTPUT 
ABT=ABT+DT»PDT 

CALL  VNULFF(CPHIL.XC,3.MM,1,3.NDA,XL,3,IER} 

RAD  =  (<XL(1)*«2)  ♦  <XL(2)*«2))*»0.5 
WRITE(6,1000)ABT,XL(1).XL(2>,XL(3),RAD 

1000  F0RNAT<1X,1F7.3,4E12.4) 

WRITE(7, 1001) ABT, RAO. XL(3) 

1001  F0RNAT<2X,1F7.3,2E14.6) 

IF(ABT.LE.TMAX)  GOTO  20 

60  CONTINUE 
RETURN 
END 
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deterministic  observers,  is  developed  and  implemented  in  a  computer  simulation. 
The  technique  of  spillover  suppression  and  the  conditions  assuring  stability  of 
the  control  system  are  developed  and  implemented  as  wall. 

The  simulation  is  tailored  to  address  the  control  of  the  Draper-2  large 
flexible  space  structure  model.  The  model  has  been  used  previously  for  optical 
pointing  (line-of-slght  <LOS)>  studies.  Here,  position  sensors  and  point  force 
actuators  are  used  to  effect  feedback  control  (regulation)  of  the  damped  un¬ 
forced  structural  vibrations.  The  simulation  can  output  both  the  unsuppreased 
and  suppressed  case  open-loop  and  closed-loop  eigenvalues  and  the  LOS  tine 
response  for  stability  and  performance  analysis. 

With  the  control  problem  formulated  for  nodal  control,  an  Investigation  is 
made  into  the  effects  on  tine  response  of  assigning  three  groups  of  four  nodes 
in  a  permutatlve  fashion  to  the  three  controllers.  A  fourth  residual  set  of 
eight  modes  is  carried  without  spillover  suppression  to  represent  the  un- 
nodelled  nodes  of  a  real  structure.  The  groups  are  assembled  based  upon  a 
previous  investigation's  results  from  applying  nodal  cost  analysis  for  LOS 
performance.  Simple  high  frequency  truncation  la  also  used  for  comparison. 

Controllers  based  on  model  coat  analysis  alone  are  found  to  yield  marginal 
stability  and  mediocre  LOS  performance  due  to  little  insight  into  the  sensitiv¬ 
ity  of  the  residual  modes  to  spillover.  However,  specific  problem  modes  are 
readily  identified  by  examining  the  results  of  an  internal  balancing  analysis 
for  nodal  sensitivity.  Simple  frequency  truncation  is  found  to  give  the  best 
time  response  here  when  the  nodes  most  contributory  to  LOS  are  assigned  to  the 
controllers  with  more  controllability.  However,  the  relatively  small  quantity 
of  nodes  and  the  overwhelmingly  large  relative  contribution  of  the  three  rigid 
body  nodes  included  nay  obscure  some  conclusions.  Results  indicate  more  re¬ 
vealing  results  night  be  obtained  if  more  modes  are  added  to  the  model  and/or 
if  some  of  the  residual  nodes  are  suppressed. 
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